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Abstract 

The  report  covers  the  period  June  1  2001  to  April  30th  2004  on  the  Office  of  Naval 
Research  Grant  N00014  01-1-870  on  the  topic  “Thermodynamic  Phenomenology  for 
Perovskite  Structure  Ferroelectric  Solid  Solutions  with  Morphotropic  Phase  Boundaries”. 
Topics  for  study  which  will  be  reported  include:  Lower  temperature  (down  to  12K) 
studies  of  the  dielectric  properties  of  single  crystal  lead  zinc  niobate:  lead  titanate. 
Studies  of  a  new  two  order  parameter  model  for  lead  zirconate  titanate  which  leads 
naturally  to  the  occurrence  of  the  lower  temperature  monoclinic  phase  interleaving  the 
tetragonal:  rhombohedral  classical  MPB.  Application  of  the  flexoelectric  measurement 
technique  developed  on  the  ONR  US;  Japan  Global  studies  program  to  explore  strain 
gradient  drive  of  ferroelastic  domain  walls  and  the  development  of  a  new  thermodynamic 
function  for  BaTiCb  to  better  describe  the  high  temperature  phase  transitions  evidenced 
under  large  epitaxial  strains  in  thin  BaTiCb  films  on  non  matching  substrates. 

1.  INTRODUCTION 

The  program  on  the  “Thermodynamic  Phenomenology  for  Perovskite  Structure 
Ferroelectric  Solid  Solutions  with  Morphotropic  Phase  Boundaries”  was  originally 
funded  under  N00014-01-1-870  for  the  two  years  01  June  2001  to  31  May  2003  and 
to  help  complete  the  proposed  studies  a  supplementary  period  from  01  June  2003  to 
April  30,  2004  was  funded.  This  final  report  will  cover  the  whole  period,  but  since 
much  of  the  earlier  work  was  described  in  the  proposal  for  supplementary  funding, 
this  focus  will  be  upon  the  completed  studies.  Following  earlier  precedent,  work 
already  published  will  be  provided  as  technical  appendices  from  which  the  important 
salient  features  will  be  discussed  in  this  text. 

2.  NEW  TWO-ORDER  PARAMETER  MODEL  FOR  LEAD  ZIRCONATE 
TITANATE  (PZT) 

The  proposed  two  order  parameter  thermodynamic  function  for  PZT  is  discussed 
in  detail  in  appendices  (1,  2,  3).  The  free  energy  in  this  model  comprises  three  sets  of 
terms,  two  expansions  of  polarization  representing  contributions  to  the  total 
polarization  of  each  end  member  in  the  solid  solution  weighted  linearly  with  respect 
to  composition,  plus  terms  representing  the  coupling  between  the  two  contributors. 
Employing  the  coefficients  for  lead  titanate  (PbTiOa)  and  extrapolated  values  for  a 
hypothetical  lead  zirconate  (PbZrOs)  taken  from  the  earlier  phenomenology  of  Haun 
et  al  (1,2, 3, 4, 5).  The  function  predicts  the  composition  of  the  Morphotropic 
ferroelectric  tetragonal:  rhombohedral  phase  boundary  to  lie  at  0.55  mole  fraction 
PbZrC>3  independent  of  the  strength  of  coupling  between  the  two  order  parameters. 
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The  model  also  predicts  a  monoclinic  interleaving  phase  at  lower  temperature  where 
composition  range  does  depend  on  the  strength  of  the  coupling  terms  between  the 
order  parameters.  The  model  with  a  suitable  choice  of  coupling  gives  a  good  fit  to 
lattice  parameters  of  the  Morphotropic  Pb(Zro.52Tio.4s)03  composition  over  a  wide 
temperature  range.  Details  are  given  in  the  appendices  1,  2,  3. 

3.  FLEXOELECTRIC  STUDIES  OF  FERROELASTIC:  FERROELECTRIC 
DOMAIN  WALL  MOTION  IN  PZT. 

During  studies  of  flexoelectric  effects  in  cubic  perovskite  paraelectrics  such  as 
(BaxSri.x)Ti03,(BST)  Pb(Mgi/3Nbi/3)Ti03  (PMN)  etc.,  which  are  in  point  group  m3m 
and  cannot  show  piezoelectricity  under  uniform  elastic  stress  it  was  pointed  out  by 
Dr.  Wenhui  Ma  that  unpoled  PZT  which  is  in  the  Curie  group  oo  oo  m  also  cannot 
show  a  piezoelectric  effect  in  a  uniform  stress  field  so  that  it  could  be  interesting  to 
explore  flexoelectricity  i.e.  charge  separation  induced  by  controlled  elastic  strain 
gradient. 

Several  questions  which  might  be  answered  by  simple  transverse  gradients  of 
strain  which  could  be  generated  in  a  flexing  bar  of  PZT  included: 

a.  Would  the  measured  values  scale  with  weak  field  dielectric  permittivity  as  in 

PMN  and  BST. 

b.  Would  the  value  of  the  scaling  parameter  y  be  similar  to  PMN,  indicating  the 

importance  of  the  lead  chemistry  or  similar  to  BST  indicating  the  importance 
of  soft  mode  behavior  as  distinct  from  relaxor  ferroelectric  response. 

c.  Since  stress  levels  could  be  easily  achieved  which  lead  to  a  remanent  curvature 

of  the  beam  indicating  ferroelastic  poling  of  the  ferroelastic:  ferroelectric 
domain  structure  would  there  be  a  domain  wall  contribution  to  the 
flexoelectric  y . 

d.  Could  the  higher  strain  gradient  levels  lead  to  a  sufficient  poling  of  the 

ferroelectric  component  in  the  domain  structure  to  give  a  measurable  residual 
piezoelectric  response. 

Detailed  answers  to  these  questions  are  given  in  appendices  4  and  5.  Briefly: 
Measured  values  do  scale  with  permittivity  but  the  scaling  factor  y  changes  with  drive 
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level.  Static  measurements  show  a  clear  onset  of  domain  wall  motion  in  the  soft  PZT- 
5H  which  coincides  with  a  change  in  elastic  stiffness.  For  low  drive  levels  the  scale 
factor  y  is  similar  to  that  in  PMN  not  BST  showing  relation  to  the  chemistry  and  not 
the  relaxor  character.  Even  for  very  high  static  stress  levels  which  clearly  evidenced 
remanent  curvature  no  measurable  piezoelectricity  was  induced. 

4.  STUDIES  RELATED  TO  PZN:PT  SINGLE  CRYSTALS 

4. 1  Orthorhombic  Phase  Stability 

In  appendix  6,  the  question  of  the  orthorhombic  phase  which  can  be  stabilized  in 
both  PZN:PT  and  PMN:PT  single  crystals  at  composition  close  to  the  MPB  is 
discussed.  It  is  pointed  out  that  even  in  the  PZT  phenomenology  the  orthorhombic 
phase  comes  very  close  in  free  energy  to  the  tetragonal  and  rhombohedral  at  the  MPB 
composition. 

4.2  Dielectric  Relaxation  in  PZN:PT  Crystals  at  Composition  0.955  PZN:  0.45PT 
Dielectric  studies  covering  the  temperature  range  from  12  K  to  300  K  are  reported 

in  appendix  7.  Two  additional  dielectric  anomalies  occurring  near  250  K  and  near  100 
K  are  reported.  The  lower  temperature  anomaly  has  Arrhenius  character  with  T  ~  10' 
115  and  activation  energy  E  -  0.24  eV,  the  upper  peak  is  not  sensitive  to  DC  bias  and 
appears  more  intrinsic  in  character.  Unipolar  high  field  studies  show  that  high  strain 
is  retaining  for  the  001  poled  crystal  into  the  cryogenic  region. 

4.3  Dielectric  Switching  in  PZN;PT  Crystals 

Anisotropy  of  fatigue  behavior  between  001  and  111  poled  PT:4.5%  PZN  is 
striking.  In  appendix  8  polarization  decay  after  1/  Ec  fields  was  explored  in  virgin 

and  105  cycled  crystals.  Relaxation  was  shown  to  lengthen  markedly  in  111  poled 
crystals  but  not  under  001  poling.  It  may  be  noted  that  in  both  orientations  switching 
is  primarily  by  domain  rotation  and  not  by  180°  switching.  For  001  poling  of  the 
rhombohedral  phase  this  is  clearly  likely,  that  is  also  occurs  for  111  poling  is  clearly 
evidence  by  the  study  described  in  appendix  9  where  the  elastic  strain  associated  with 
reversal  clearly  indicates  that  it  is  not  by  180°  wall  motion. 

4.4  Hydrostatic  dn  in  PZT,  PZN:PT  and  PYN:PT 

It  has  long  been  realized  that  ferroelastic  ferroelectric  domains  can  not  be  driven 
by  uniform  hydrostatic  stress  so  that  for  ferroelectric  materials  where  domains  can 
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contribute  strongly  to  direct  and  converse  piezoelectric  d33  and  d3i  in  uniaxial  phases, 
dn  as  measured  under  hydrostatic  stress  may  not  be  equivalent  to  (d33  +2d3i)  as 
measured  by  resonance  methods.  It  must  be  remembered  also  however  that  in  relaxor 
ferroelectric  crystals  the  nano-polar  regions  which  can  persist  in  both  paraelectric  and 
ferroelectric  phases  do  change  orientation  under  electric  field  and  can  contribute 
markably  to  dielectric  response.  It  is  very  clear  from  these  true  hydrostatic 
measurements  in  Appendix  10  that  in  all  the  systems  studied  there  is  a  discrepancy 
between  dn  and  (d33  +  2  d3i).  In  soft  PZT  this  is  most  probably  a  true  domains  wall 
contribution.  PZN:PT  and  PYN:PT  are  however,  both  embracing  relaxor  end 
members.  Now  there  is  much  more  massive  disagreement  but  since  nanopolar  regions 
which  also  will  not  be  driven  by  hydrostatic  stress  are  strongly  associated  it  is  not 
clear  whether  this  should  be  classed  as  an  extrinsic  effect  or  not. 

5.  BARIUM  TITANATE  PHENOMENOLOGY 

To  uncover  the  reason  for  the  disagreement  between  predicted  and  observed 
phase  transitions  driven  by  electric  field  in  1 1 1  oriented  BaTi03  single  crystals,  it  was 
hoped  to  make  a  detailed  study  of  the  Fujikura  BaTi03  crystals  used  in  the  111  high 
field  studies  at  TIT  in  Japan.  Unfortunately  due  to  personnel  changes  on  the  contract 
this  became  impossible.  It  was  also  clear  however,  that  Bell/Cross  free  energy 
function  used  in  the  analysis,  which  was  designed  primarily  to  describe  dielectric  and 
electric  field  behavior  near  the  10°  C  and  -90°  C  orthorhombic  and  rhombohedral 
transitions  was  becoming  unusable  at  higher  temperatures.  The  need  for  a  modified 
form  is  highlighted  by  recent  thin  film  studies  where  induced  strain  associated  with 
epitaxial  constraint  on  a  SiTiC>3  substrate  is  sufficient  to  raise  Tc  to  more  than  950  K. 
(Appendix  11)  To  explore  this  higher  temperature  region  a  new  energy  function 
embracing  up  to  eighth  order  teams  has  been  developed  (Appendix  12).  In  spite  of  the 
fact  that  only  the  second  order  terms  are  linearly  temperature  dependent  the  function 
nicely  describes  the  entropy  changes  at  al  three  ferroelectric  transitions  and  enables 
prediction  of  higher  temperatures  high  strain  properties.  Unfortunately  the  spread  of 
current  values  in  the  measured  electrostrictive  Q  constants  still  leaves  a  wide  spread 
in  predicted  changes  of  Tc  (6)  and  do  clearly  highlight  the  need  for  more  precise 
measurements  which  are  now  well  within  the  capability  of  modern  ultra  dilatometers. 
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6.  POLARIZATION  SWITCHING  IN  PbCYbmNb^Oa-PbTiOa  SINGLE 
CRYSTALS 

Evidence  of  the  importance  of  disorder  in  the  Yb/Nb  arrangement  in  the  high 
temperature  grown  Pb(Ybi/2Nbi.2)03-PbTi03  crystals  is  addressed  from  the 
ferroelectric  switching  behavior  discussed  in  Appendix  13.  This  important  paper 
highlights  the  role  of  random  fields  and  nano-polar  direction  in  the  switching  process, 
and  the  strong  role  which  they  play  in  the  single  crystal  systems  which  exhibit 
massive  piezoelectric  properties. 
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A  two-parameter  thermodynamic  model  is  proposed  for  the  ferroelectric  solid  solution  system 
Pb(Zri_xTi*)C>3.  The  free  energy  comprises  three  sets  of  terms:  two  expansions  of  polariza¬ 
tion,  one  representing  each  of  the  end-members,  weighted  linearly  with  respect  to  composition, 
and  terms  representing  the  coupling  between  the  two  polarization  contributions.  Employing 
previously  published  coefficients  for  PbTi03  and  values  extrapolated  from  Zr-rich  PZT  com¬ 
positions  for  PbZrC>3,  calculations  predict  the  morphotropic  phase  boundary  lies  at  x  =  0.45 
and  is  independent  of  the  strength  of  coupling  between  the  two  order  parameters.  A  mono- 
clinic  phase  is  predicted,  the  range  of  existence  of  which  depends  upon  the  strength  of  the 
parameter  coupling,  but  in  general,  is  broadest  at  low  temperature.  It  is  proposed  that  the  two 
contributions  to  the  polarization  have  their  origins  in  the  two  types  of  cation  shifts  identified 
in  PZT:  the  major,  homogenous  shifts  parallel  to  the  polar  direction  and  randomly-oriented,  or 
heterogeneous,  minor  shifts  away  from  the  polar  axis.  The  field-induced  ordering  of  the  minor 
shifts  is  identified  as  the  origin  of  the  large  intrinsic  piezoelectric  coefficient  in  PZT. 

Keywords:  Thermodynamic  model;  PZT 


INTRODUCTION 

Lead  zirconate-titanate,  Pb(Zri„xTix)03  or  PZT  [1],  has  been  widely  stud¬ 
ied  due  to  the  outstanding  piezoelectric  properties  of  compositions  close  to 
x  =  0.5.  The  peak  in  the  piezoelectric  coefficient  for  these  compositions 
has  often  been  associated  with  the  presence,  at  x  =  0.48,  of  a  morphotropic 
phase  boundary  between  the  zirconium-rich  rhombohedral  perovskite  and 
the  titanium-rich  tetragonal  perovskite  phases.  Haun  [2-6]  successfully  mod¬ 
elled  the  PZT  phase  diagram,  as  it  had  been  known  for  30  years,  using  the 
thermodynamic  approach  of  Devonshire  [7],  in  which  the  free  energy  associ¬ 
ated  with  the  polar  order  parameter  is  expressed  in  terms  of  a  power  series  of 
the  pseudo-cubic  components  of  the  polarization.  A  sufficiently  accurate  fit 
to  experimental  data  has  usually  been  achieved  with  expressions  expanded 
up  to  the  6th  power  of  polarization. 
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figure  1.  Revised  phase  diagram  of  PZT,  by  combining  those  of  Jaffe  [1]  and 
Noheda  |8|. 

The  recent  discovery  of  a  monoclinic  phase  at  the  morphotropic  phase 
boundary  in  PZT  [8]  (Fig.  1)  and  the  confirmation  that  other  systems  with 
morphotropic  phase  boundaries  that  exhibit  very  large  piezoelectric  effects 
191  also  exhibit  other  phases  close  to  the  phase  boundary  [10,  1 1],  has  led  to 
renewed  interest  in  how  best  to  represent  ferroelec tries  with  morphotropic 
phase  boundaries  in  the  context  of  the  Landau-Devonshire  model.  The  mon- 
oclinie  phase  in  PZT  was  the  first  known  phase  in  a  ferroelectric  perovskite 
in  which  die  non-zero  components  of  the  polarization  with  respect  to  the 
cubic  axes  are  not  equal  under  conditions  of  zero  applied  field  and  stress. 
Such  phases  are  not  solutions  of  the  conventional  6th  power  Devonshire 
model,  however  Vanderbilt  [12]  has  shown  that  the  inclusion  of  8th  power 
terms  in  the  perovskite  free  energy  expansion  can  lead  to  the  possibility  of 
monoclinic  phases,  whilst  the  inclusion  of  12th  power  terms  allows  triclinic 
solutions.  The  suggestion  is  that  due  to  the  implicit  disorder  in  solid  solu¬ 
tions  such  as  PZT,  the  8th  order  terms  may  assume  a  greater  importance 
than  in  single  component  perovskite  ferroelectrics.  However,  the  8th  power 
equations  do  not  explicitly  address  the  disordered  nature. 

la  such  treatments,  the  composition  dependence  of  the  free  energy  co- 
el  ticients  is  not  implicit  in  the  model,  the  locus  of  the  morphotropic  phase 
boundary  therefore  being  an  input  to,  rather  than  output  from,  the  exercise. 
Neither  are  the  mesoscopic  features  of  solid  solutions  represented,  an  ex¬ 
ample  being  the  local  cation  displacements  in  PZT.  It  would  appear  that  the 
“local  symmetry”  ot  the  lead  cation  is  monoclinic  in  both  the  rhombohedral 
and  tetragonal  phases.  For  the  rhombohedral  phase,  it  has  been  proposed 
that  the  displacement  of  the  Pb  cations  from  their  cubic  position  comprises  a 
homogeneous  [111]  displacement  superimposed  with  minor,  random  (001) 
shii'is  1 13|.  la  ihe  tetragonal  phase,  the  homogeneous  [001]  displacements 


are  augmented  with  minor  (110)  (or  equivalent  (111))  random  shifts  [14]. 
In  both  cases,  the  random  shifts  average  to  zero  so  that  the  macroscopic 
symmetries  remain  rhombohedral  and  tetragonal  respectively.  The  random 
Pb  displacements  might  be  considered  as  either  static  shifts  or  dynamic 
fluctuations  about  the  homogeneous  displacement.  It  is  not  known  whether 
they  are  distributed  totally  randomly,  or  whether  their  directions  are  partially 
correlated  in  clusters  and  whether  their  magnitudes  are  equal  on  all  sites. 
Nevertheless,  it  seems  clear  that  the  influence  of  PbTiC>3  exists  well  into 
the  rhombohedral  phase  field  and  that  of  PbZrC>3  into  the  tetragonal  phase. 
Moreover,  in  general,  the  Ti  and  Zr  ions  do  not  have  equal  displacements  in 
either  phase  [13,  14]. 

The  minor  displacements  can  be  assumed  to  be  of  some  significance  in 
the  appearance  of  the  monoclinic  phase.  For  example,  it  might  be  postulated 
that  the  transition  from  the  tetragonal  to  monoclinic  phase  in  PZT  occurs 
due  to  the  condensation  of  the  random  (1 1 1)/{1 10)  displacements,  which 
combined  with  the  homogeneous  [001]  shifts  give  rise  to  a  macroscopic 
monoclinic  symmetry.  As  the  Pb  cation  shifts  are  often  cited  as  being  the 
main  component  of  dipole  moment  in  Pb  perovskite  ferroelectrics,  it  may  be 
possible  to  construct  thermodynamic  expressions  to  reflect  this  postulate,  in 
which  both  the  homogeneous  and  random  contributions  to  the  polarization 
are  represented,  or  in  which  the  symmetry  of  the  two  end  members  is  allowed 
to  influence  the  stable  state  well  inside  the  solid  solution. 

Investigations  of  the  composition  dependence  of  the  optical  phonon  fre¬ 
quencies  in  PbZri_xTix03  tend  to  support  the  concept  of  their  being  two 
contributions  to  the  polarization  in  the  solid  solution.  Infra-red  measure¬ 
ments  by  Sivasubramanian  [15],  on  the  tetragonal  side  of  the  phase  diagram 
0.475  <  x  <  1,  show  that  the  E(ITO)  mode  softens  with  increasing  Zr  con¬ 
centration.  In  PbTi03,  this  mode  corresponds  to  the  motion  of  the  Ti  ion 
against  the  oxygen  octahedra  and  has  been  identified  as  the  soft  mode  re¬ 
sponsible  for  the  ferroelectric  transition  [16].  Its  decrease  in  frequency  with 
increasing  Zr  concentration  suggests  that  the  transition  from  the  tetragonal  to 
the  rhombohedral  (or  monoclinic)  phase  is  also  driven  by  an  instability  in  this 
mode.  However,  Raman  measurements  [17]  on  a  similar  set  of  compositions 
at  low  temperatures  (4.2  to  300  K)  reveal  that  for  x  <  1  there  is  a  splitting 
of  the  three  lowest  E(TO)  modes,  not  observed  in  PbTi03.  The  frequency 
difference  due  to  the  splitting  increases  with  increasing  Zr  concentration 
and,  in  the  case  of  the  E(ITO)  mode,  both  branches  decrease  in  frequency 
with  increasing  Zr  concentration  towards  the  phase  boundary  region.  The 
experimental  data  appear  to  indicate  that  the  splitting  is  present  at  least  up  to 
room  temperature  and  that  both  branches  soften  with  increasing  Zr  content. 


but  there  is  no  Jala  above  room  temperature  to  show  whether  both  branches 
soften  through  Tc.  The  more  recent  Raman  study  by  Lima  [18]  at  7  K  con¬ 
centrates  less  on  the  lowest  frequency  mode  due  to  problems  of  resolving 
the  spectra,  but  shows  how  the  higher  modes  vary  in  frequency  as  the  mor- 
pholropic  region  is  traversed.  Although  a  full  experimental  survey  has  yet  to 
be  completed,  the  implication  is  that  the  splitting  of  the  low-lying  TO  modes 
may  contribute  to  the  polarization  in  a  less  than  straightforward  manner. 

Here  a  thermodynamic  model  is  presented  in  which  the  free  energy  of 
FZT  is  expressed  as  an  expansion  in  terms  of  two  polarization  contributions, 
representing  the  different  influences  of  the  two  end  members  of  the  solution 
series.  It  is  based  upon  a  linear  combination  of  the  free  energy  expansions 
of  the  end-members  of  the  solid  solution,  together  with  linear  and  quadratic 
coupling  terms.  The  total  polarization  is  a  function  of  the  two  end  member 
contributions  and  the  composition  variable. 

THEORY 

In  general,  the  free  energy,  AG,  of  an  infinite,  homogeneous,  single  domain 
lerroelcctric  with  a  cubic  prototype  structure,  as  in  the  perovskites,  is  given 
as  an  expansion  of  the  orthogonal  components  of  the  polarization  P  in  which 
only  the  coefficient  <*200  necessarily  temperature  dependent,  with  the  form 
^200  =  1*201)0'  fuk 


P)  =  «200(P?  +  P2  +  P3)  +  0(400 (Pj  +  Pj  +  Pb) 

+  «220(P?Pl  +  P^B  +  P3P')  +  «600(Pt  +  P2  +  P3) 
+  <*42o[P?(Pl  +  P§)  +  Pi(Pl  +  P?)  +  PUP?  +  P22)] 
+  £*22iP?P!p3 . -B,p,  -  E2P2  -  E3P3, 


where  L,s  the  applied  electric  field.  Expansions  terminated  at  the  6th  pow 

polunzaUo"  have  generally  been  found  sufficient  to  provide  an  accura 

tec„p,u,„  Ot  inosl  ferroelectric  perovskiie  systems.  However  “°Z  “ns 

fond  "  01  SOl,d  solul,ons  discussed  above,  a  solid  solution  enen 

function  is  proposed  of  the  form:  * 


AGss(x)  _  xGprtoWp.  p)  +  (1  -  x)GP2(£mnp,  q)  +  Gc(ynmp,  p,  q)  (2) 


where  x  is  the  mole  fraction  of  PbTi03.  The 
the  contribution  to  the  polarization  from  the 


vectors  p  and  q  represent 
two  end-members,  PbTi03 


and  PbZrC>3,  respectively.  The  forms  of  GPT(tfmnp,  p)  and  GPz(/3mnp,  q)  are 
identical  to  that  of  Eq.  1,  hence  the  symmetry  requirements  of  the  cubic 
prototype  are  maintained.  The  sets  of  coefficients,  amnp  and  /Jmnp,  are  those 
of  the  end  members  of  the  solid  solution,  in  which  only  To  is  a  function  of 
composition  and  is  common  to  both  sets,  in  the  form: 


To(x)  =  x(Tpx  —  Tpz)  +  Tpz  (3) 

where  Tpr  and  TPz  represent  the  Curie  temperatures  of  PbTi03  and  PbZr03 
respectively. 

The  coupling  energy,  for  terms  up  to  the  4th  power  of  polarization,  is  of 
the  form: 

Gc(Kmnp.  p.  q)  =  yiooCpiqi  +  p2q2  +  p3q3)  +  X20o(p?qi  +  pW2  +  P3q3) 

+  X220 (p? (q2  +  q3)  +  p2(q3  +  q>)  +  p3(q?  +  q2))  (4) 

Terms  in  p^qj,  where  i  ^  j,  are  disallowed  by  symmetry,  as  the  free  energy 
should  be  independent  of  the  sign  of  pjqj.  However,  terms  in  pjqi  are  nec¬ 
essary  to  differentiate  between  the  parallel  and  anti-parallel  states  of  pj  and 
qi.  Table  I  lists  the  non-degenerate  states  resulting  from  combination  of  the 
components  of  p  and  q,  using  the  notation  of  Vanderbilt  [12]  for  the  result¬ 
ing  phases.  Due  to  the  coupling  energy  Gc,  the  components  of  p  and  q  are 

TABLE  I  Classification  of  phases  for  non-degenerate 
combinations  of  pj  and  qj.  Only  components  constrained  to  be 
zero  are  shown  as  such.  If  p*  and  qt  are  both  shown  as  non-zero, 
either,  but  not  both,  may  assume  a  zero  value 
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noi  constrained  to  the  6lh  power  solutions  of  the  end-member  contributions. 
For  directions  in  which  the  component  of  the  total  polarization,  Pi,  is  not 
constrained  by  symmetry  to  be  zero,  either  of  the  contributions,  pj  or  qit  is 
allowed  to  be  zero  as  long  as  p;  +  qi  ^0.  It  should  be  noted  that  all  of  the 
phases  identified  by  Vanderbilt  [12]  for  the  conventional  12th  power  Devon¬ 
shire  model  are  possible  solutions  of  the  current  approach,  including  three 
non-degenerate  triclinic  phases. 

For  comparison  with  quantitative  available  data,  it  is  necessary  to  evaluate 
the  total  polarization  P,  which  is  a  function  of  p  and  q.  For  this  purpose,  it 
is  assumed  here  that  the  relationship  is  of  the  simple  form: 

P  =  kp(x)  p  +  k4(x)  q.  (5) 

Two  possible  assumptions  are  (i)  kp  =  kq  =  1  or  (ii)  kp  =  x  and  kq  =  1  —  x. 
The  first  assumption  implies  that  p  and  q  are  assigned  homogeneously  across 
all  unit  cells,  independent  of  B-site  occupation,  i.e.  the  dipole  moment  of  all 
unit  cells  is  equal.  However,  it  is  clear  from  the  structure  data  of  Corker  [13] 
and  Fruntti  [141  that  this  is  unlikely,  as  the  Zr  and  Ti  displacements  are  often 
different  within  a  single  composition.  Assumption  (ii)  implies  that  p  and  q 
are  distributed  according  to  the  B-site  occupancy  of  the  unit  cell.  This  would 
mean  that  in  the  example  addressed  below,  the  dipole  associated  with,  say, 
Zr  cells  in  the  tetragonal  phase,  would  be  zero.  This  is  also  not  supported 
by  the  structure  data  [13,  14].  Hence,  the  evaluation  of  kp  and  kq  must  at 
present  be  undertaken  by  comparison  with  experimental  data. 


CALCULATIONS  FOR  PZT 

The  model  presented  here  is  aimed  at  the  correct  prediction  of  the  mor 
pilot  topic  phase  boundary  and  the  appearance  of  perovskite  ferroelectrii 
states  with  non-equal  polarization  components  i.e.  monoclinic  and  triclini 
phases.  Hence,  the  orthorhombic  antiferroelectric  phase  in  PZT  is  ignored 
and  lor  the  purpose  of  these  calculations,  lead  zirconate  is  assumed  to  b 
a  .  homhohcdrul  ferroelectric.  Similarly,  the  octahedral  tilt  transition,  whicl 
uilleientiates  the  high  and  low  temperature  rhombohedral  phases  of  PZT  i 
not  addressed  here. 

Calculations  of  the  temperature  and  composition  dependence  of  p  and  < 
and  hence  (he  phase  boundaries  were  carried  out  using  the  above  model  fo 
^eto  applied  held.  Choosing  the  A  end-member  in  Eq.  2  to  be  PbTi03  thi 
sohmon  lor  x  =  1  should  yield,  for  T  <  TA,  p,  =  0.  p2  =  0,  p3  =£  0  q  =  0 

'  °  “  lC"a8H,lal  ferroelectric,  whereas  x  =  0  should  yield  solutions  for  i 


TABLE  II  Free  energy  coefficients  for  PbTi03  (amnp) 
and  hypothetically  ferroelectric  PbZrCh  (/Jmnp) 


mnp 

^mnp 

Annp 

200 

3.74  X  103  (T  -  To) 

2.82  x  105  (T  -  T0) 

400 

-7.9  x  107 

5.12  x  108 

220 

7.5  x  108 

-6.5  x  108 

600 

2.61  x  10s 

V5.93  x  10“ 

420 

6.3  x  10s 

v2  x  109 

222 

-3.66  x  109 

-9.5  x  109 

rhombohedral  ferroelectric,  i.e.  for  T  <  Tb,  p  =  0,  qi  =  q2  =  q3  ^  0. 
The  amnp  coefficients  for  PbTi03  used  for  these  calculations  are  those  given 
by  Haun  [4]  and  are  shown  in  Table  II.  The  /?mnp  coefficients  for  a  hypothet¬ 
ically  ferroelectric  PbZr03  have  been  estimated  by  extrapolating  the  Haun 
coefficients  for  compositions  with  x  >  0  back  to  x  =  0.  In  some  cases,  notably 
P200  and  ^400,  the  extrapolated  value  is  identical  to  that  of  the  original  Haun 
coefficient,  whereas  the  remaining  coefficients  tend  to  a  linear  dependence 
on  x  except  close  to  x  =  0  and  x  =  0.5. 

The  calculations  were  carried  out  by  numerical  minimization  of  AGss(x) 
with  respect  to  p  and  q  using  the  FindMinimum  function  of  Mathematica 
4.0.  Only  positive  values  of  pi  and  q*  were  considered,  hence  in  this  context 
the  F100  term  IS  redundant  and  for  simplicity  its  value  was  set  to  zero,  leaving 
only  y2oo  and  ^220  as  the  unknown  parameters. 

The  problems  of  identifying  the  absolute  free  energy  minimum  in  six 
dimensional  polarization  space  are  not  trivial,  hence  solutions  were  sought 
which  corresponded  to  the  most  likely  phases  in  Table  L  For  each  temperature 
and  PbTi03  molar  fraction,  the  local  minima  in  AGss(x)  were  systematically 
identified  for  pi  >  Oandq,  >  0,  with  the  constraints  shown  in  Table  I,  for  the 
cubic,  tetragonal,  orthorhombic,  rhombohedral,  monoclinic  Ma  and  mono¬ 
clinic  Me  phases.  The  energies  of  the  minima  were  then  compared  to  define 
the  stable  state. 

The  strain  matrix,  o',  is  calculated  according  to: 

cry  =  Qijki  PkPi  (6) 

For  reasons  of  symmetry,  there  are  only  three  independent  components  of  the 
electrostriction  coefficient  tensor  Q,  which  in  reduced  notation  are  Q 1 1 ,  Q 12 
and  Q44.  Their  values  as  a  function  of  x  were  derived  by  linear  interpolation 
from  those  of  the  end-members,  published  by  Haun  [4]. 


Calculation  of  the  angular  dependence  of  the  piezoelectric  coefficient, 
which  is  of  relevance  to  understanding  the  properties  of  polycrystalline  bod¬ 
ies,  has  been  carried  out  using  the  current  model.  The  electric  field,  E,  is 
applied  at  angle  6  to  the  [001]  axis  and  0  to  the  [100]  axis  in  the  (001)  plane. 
The  components  of  the  field  parallel  to  the  three  principle  axes  may  then  be 
calculated  and  included  in  the  free  energy  (Equation  2)  as  terms  inEj  pj  and 
Ej  qj.  The  values  of  p(E)  and  q(E)  under  applied  field  can  thus  be  found,  as 
may  the  strain  components,  cry  (Equation  6),  with  the  assumption  that  P(E) 
=  kp  p(E)  4~  kq  q(E).  Transformation  of  the  strain  tensor  to  calculate  the 
strain  parallel  to  the  applied  field,  o^,  is  accomplished  via: 

a0(p  =  a j  cos2  9  —  cos  0  sin  9  fa  cos  0  +  05  sin  0) 

+  sin2  0(<j|  cos2  04-0-4  cos  0  sin  0  4-  02  sin2  0),  (7) 

from  which  the  piezoelectric  coefficient  can  be  calculated  from  the  difference 
between  the  zero-field  and  field-induced  strains. 

RESULTS  AND  DISCUSSION 

The  stabilities  of  the  phases  were  determined  as  a  function  of  }^oo  and 
}  '220 •  In  the  examples  discussed,  only  stable  phases  corresponding  to  cubic, 
tetragonal,  rhombohedral  and  monoclinic  MA  were  found.  Phase  diagrams 
as  a  function  of  x  and  T  are  shown  for  various  combinations  of  y200  and  K220 
in  Figs.  2,  3,  and  4. 

For  y2(H)  —  >'220  5  0.  there  is  only  a  single  phase  boundary  at  To(x),  be¬ 
tween  cubic  and  monocliuic  MA,  for  all  PbTiC>3  concentrations  except  the 


Figure  2.  Phase  diagram  for  y2oo  =  /220  — >•  00. 


Figure  3.  Phase  diagram  for  y2 oo  =  Y220  =  10s. 

end-members.  For  non-zero  values,  ^200  and  K220  must  exceed  a  threshold 
value  to  avoid  discontinuities  as  x  — >  0  and  x  —>  1 . 

When  both  froo  and  y220  exceed  2  x  I08  the  appearance  of  the  monoclinic 
phase  is  completely  inhibited  see  (Fig.  2).  The  locus  of  the  morphotropic 
boundary  is  independent  of  y2oo  and  y220  and,  well  below  T0,  it  is  also 
almost  independent  of  temperature,  moving  to  slightly  higher  values  of  x  with 
decreasing  temperature,  and  extrapolating  to  x  =  0.45  at  0  K.  On  approaching 
To,  the  boundary  curves  sharply  towards  lower  values  of  x;  the  point  at  which 
the  morphotropic  phase  boundary  finally  coincides  with  the  cubic  phase 
boundary  is  at  x  =  0.  Hence  a  cubic-tetragonal  transition  occurs  on  cooling 


Figure  4.  Phase  diagram  for  y2 00  =  6  x  108,  y22o  =  108. 


for  ail  compositions  with  x  >  0,  however  for  compositions  up  to  x  =  0.35  the 
tetragonal  phase  is  only  stable  for  up  to  5  K  below  To,  before  transforming 
to  the  rhombohedral  phase.  The  paraelectric-ferroelectric  phase  transition  is 
1st  order  for  all  values  of  x. 

For  intermediate  values  of  the  coupling  coefficients,  the  MA  phase  ap¬ 
pears  in  the  centre  of  the  phase  diagram.  On  increasing  the  y  values  from 
zero,  there  is  an  almost  symmetrical  reduction  in  the  width  of  the  monoclinic 
phase  field,  yielding  to  the  rhombohedral  phase  at  low  x  and  the  tetragonal 
phase  at  high  x  values.  There  is  also  a  reduction  in  the  maximum  temperature 
of  nionoclinic  stability  which  occurs  close  to  x  =  0.5.  Thus  as  the  coupling 
energy  is  increased,  the  morphotropic  phase  boundary  between  the  rhom¬ 
bohedral  and  tetragonal  phases  is  revealed,  with  the  MA  phase  appearing 
as  a  wedge,  with  its  apex  close  to  the  morphotropic  phase  boundary.  The 
example  shown  in  Fig.  3  is  for  }'2oo  =  K220  =  108,  whilst  that  of  Fig.  4  is  for 
)'200  =  6  x  I08  and  >'220  =  108,  where  the  inequality  in  the  coupling  coeffi¬ 
cients  results  in  the  monoclinic  phase  favouring  the  “tetragonal”  side  of  the 
morphotropic  phase  boundary. 

In  the  absence  of  any  reliable  measurements  of  polarization  on  single 
crystals  of  PZT,  the  values  of  kp  and  k^  may  be  estimated  by  comparison  of 
the  values  of  pj  and  qj  calculated  here  with  the  values  of  P*  calculated  from  the 
Haun  model,  which  are  thought  to  be  an  accurate  representation  of  the  single 
crystal  polarization.  Figure  5  shows  the  ratio  of  P3  from  the  Haun  model  to 
PjT  q3  from  the  two-parameter  model  as  a  function  of  composition,  both 
calculated  at  300  K.  Only  points  in  the  rhombohedral  and  tetragonal  fields 


Figure  5.  Ratio  of  the  polarization  along  [001]  from  the  Haun  model  to  that  ealeu- 
lated  From  the  current  model. 
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Figure  6.  Calculated  piezoelectric  coefficient  as  a  function  of  the  applied  field  at 
250  K  in  units  of  pm  V~J  for  (a)  Pb(Zr0.6Tio.4)03  and  (b)  Pb(Zr04Tio.6)03;  the  [001] 
axis  is  vertical. 

were  considered  so  that  p3  and  q3  were  never  simultaneously  non-zero.  In  the 
rhombohedral  phase,  the  ratio  P3/q3(—  kq)  is  seen  to  be  constant  with  a  value 
of  0.67.  In  the  tetragonal  phase  the  ratio  P3/p3(=  kp)  falls  from  a  maximum 
of  1  at  x=  1  to  a  value  of  approximately  0.7  at  x  =  0.5.  The  trend  line  is 
quadratic  withkp  =  — 0.822x2  4-  1.8324x  and  coincides  with  the  value  0.67 
at  x  =  0.46,  which  might  be  interpreted  as  the  phase  boundary  composition. 

Examples  of  the  angular  dependence  of  the  piezoelectric  coefficient  for 
Pb(Zro.4Tio.6)03  and  Pb(2To.6Tio.4)03  are  shown  in  Fig.  6  and  are  similar 
to  those  calculated  by  Du  et  al.  [16]  from  the  Haun  model.  It  is  possible 
to  determine  which  of  the  two  parameters,  p  and  q,  contributes  most  to  the 
piezoelectric  coefficient  as  their  contributions  to  the  induced  strain  can  be 
examined  separately.  In  general  terms,  the  induced  strain  is  proportional  to 
P(E)2-P(0)2,  which,  assuming  that  P(E)  =  kp  p(E)  4-  kq  q(E),  expands  to: 

k^p(E)2  +  2kpkqp(E)q(E)  +  k2q(E)2  -  k^p(0)2  -  2kpkqP(0)q(0)  -  k2q(0)2. 

(8) 

However,  in  the  tetragonal  phase,  q(0)  =  0,  hence  there  are  only  three  terms: 
(i)  p(E)2  -  p(0)2,  (ii)  q(E)2,  and  (iii)  p(E)  q(E).  A  corresponding  argu¬ 
ment  applies  to  the  rhombohedral  phase,  with  p(0)  =  0.  Interestingly,  for  all 
the  compositions  and  temperatures  examined  numerically,  it  is  the  term  in 
p(E)  q(E),  which  dominates  the  value  of  the  piezoelectric  coefficient.  This 
implies  that  the  large  piezoelectric  coefficients  found  in  the  MPB  region 
are  due  to  the  field-ordering  of  the  random  components  of  the  polarization, 


“amplified"  by  the  homogeneous  component.  A  similar  argument  may  also 
applied  to  the  permittivity.  As  the  morphotropic  boundary  is  approached,  the 
random  components  of  polarization  (or  cation  displacement)  assume  greater 
importance,  culminating  in  the  well-reported  peak  in  permittivity  and  piezo¬ 
electric  coefficient  at  the  boundary.  According  to  the  model  proposed  here, 
the  phenomenon  is  therefore  essentially  due  to  the  local  random  nature  of 
the  PZT  solid  solution. 

CONCLUSIONS 

A  novel  thermodynamic  model  has  been  proposed  to  describe  the  properties 
ofPb(Zr,Ti)03  solid  solutions,  particularly  in  the  region  of  the  morphotropic 
phase  boundary.  Based  on  an  expression  which  sums  the  two  end  member 
contributions  to  the  free  energy  in  proportion  to  their  concentration,  the 
model  predicts  the  position  of  the  morphotropic  phase  boundary  at  x  =  0.45. 
This  should  be  compared  to  the  experimental  locus  of  x  =  0.48.  The  mor- 
pholropic  phase  boundary  is  independent  of  any  of  the  adjustable  parameters 
in  the  model  and  hence  the  good  agreement  with  experimental  data  appears 
to  confirm  the  suitability  of  the  approach  taken.  The  range  of  existence  of 
the  monoclinic  phase  is  controlled  by  the  coupling  energy  between  the  two 
end-member  contributions  to  the  polarization,  p  and  q.  At  zero-field,  p  and 
q  coexist  only  in  the  monoclinic  phase.  Although  the  model  allows  for  the 
appearance  of  the  three  possible  perovskite  monoclinic  phases,  only  the 
phase  was  predicted  by  the  model  for  PZT,  again  in  keeping  with  experiment. 
In  the  examples  here,  only  two  coefficients  were  used  to  describe  the  cou¬ 
pling  energy.  By  increasing  the  coupling  between  the  parallel  components  of 
p  and  q,  whilst  keeping  the  orthogonal  coupling  energy  at  moderate  levels, 
the  Ma  phase  is  restricted  to  the  PbTi03-rich  side  of  the  morphotropic  phase 
boundary,  as  suggested  by  experimental  data  [19]. 

Examination  of  the  contributions  to  the  piezoelectric  coefficient  from  the 
two  polarization  parameters  reveals  that  the  cross  term  between  the  two  pa¬ 
rameters  is  dominant,  implying  that  the  field  stabilization  of  the  random  fluc¬ 
tuations  is  the  origin  of  the  large  intrinsic  piezoelectric  coefficients  in  PZT. 
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Abstract  -  A  thermodynamic  model  k  proposed  to  account 
for  the  monoclmk  phase  observed  in  the  solid  solution 

system  PbCZr^TiJO,.  The  free  energy  comprises  three 
sals  of  tanas:  two  expansions  of  polarisation,  one 
npmeotint  oaeh  of  tho  end  members,  weighted  lineeriy 
with  respect  to  conkpasdtofi,  and  terms  rupee  tenting  the 
costing  between  the  two  contributions.  Employing 
previously  pubHshad  eosffMeott  for  PtTiOj  and  vahat 
extrapolated  from  Zr-ricb  PZT  compositions  for  PbZiO* 
calculations  predict  foe  composition  of  foe  morphotropic 
phase  boundary  between  foe  ihomboheftel  and  tetragonal 
phases  lies  at  x=0.45  and  is  virtually  independent  of 
temperature  and  foe  strength  of  coupling  between  the  two 
order  pnoatart,  The  reage  of  existence  of  foe 

monocHnic  phase  depends  upon  the  strength  of  the 
coupling.  but  in  general,  is  broadest  at  low  temperature 


INTRODUCTION 

Lead  zircoostotilazMe,  Pb(Zr^T\A  or  PZT  (11,  has 
been  widely  studied  due  to  foe  outstanding  piezoelectric 
properties  of  compositions  close  to  x  »  0.5.  The  peak  in 
the  piezoelectric  coefficient  for  these  compositions  has 
often  been  associated  with  foe  pretence,  at  x  *  04t»  of  a 
matphocropic  (toe  boundary  between  the  zirconium-rich 
rhombohadral  perovakite  and  the  titanium-rich  tetragonal 
perovskhe  phases.  Hsua  [2-6]  successfully  modelled  foe 
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Fig  1 .  Revised  phase  diagram  of  PZT,  by  combining 
those  of  Jaffe  [1]  and  Noheda  [S]. 


parameter  ts  expressed  in  tenm  of  a  power  series  of  foe 
pwdfrcabk  comperwntt  of  foa  polartzetioo,  A 
sufficiently  accurate  ft  to  experimental  data  has  usually 
been  achieved  with  expressions  expanded  up  to  the  6* 
powar  of  polarisation. 

The  recent  discovery  of  a  monoctmac  plume  at  foe 
morphotropic  phase  boundary  in  PZT  [*}  (Figure  1)  and 

foe  ©oofirtnatioo  that  ofoar  ayiMni  with  morphotropic 

phase  boundaries  that  exhibit  very  large  piezoelectric 
effects  [9]  she  exhibit  other  phases  close  to  foe  phase 
boundary  [10  A  If],  has  fed  to  renewed  interest  in  how 
best  to  represent  ferroelectric*  with  morphotropic  phase 
boundaries  in  foe  context  of  the  Landau-Devonshhe 
model.  The  monoclinic  phase  in  PZT  was  foe  first  known 

phase  to  a  ferroelectric  perovskite  in  which  foe  non-zero 
components  of  foe  polarization  with  respect  to  the  cubic 
axes  are  not  equal  under  conditions  of  zero  applied  field 
and  stress.  Such  phases  are  not  solutions  of  foe 
conventional  6*  power  Devonshire  model,  however 
Vandtsbttt  [12]  has  shown  foal  foe  inclusion  oft*  power 
terms  in  foe  perovskite  free  energy  expansion  can  lead  to 
foe  possibility  of  monocliwc  phases,  whilst  foe  inclusion 
of  IT*  power  term*  allows  triclinic  aohirinvvL  The 
suggestion  is  foot  due  to  foe  kophch  disorder  to  solid 
solutions  such  as  PZT,  the  I*  order  terms  may  assume  a 
greater  importance  than  in  single  component  perovskhe 
fcrroekctrica.  However,  the  t*  power  equations  do  not 

explicitly  address  the  dbonfcied  nature. 

hi  such  treatments,  foe  composition  dependence  of  foe 
free  energy  coefficients  is  not  implicit  in  foe  model,  foe 
locus  of  foe  morphotropic  phase  boundary  therefore  being 
an  input  to.  rather  than  output  from,  the  exercise.  Neither 
we  foe  mesoscopic  features  of  solid  solutions  represented, 
fat  PZT,  close  to  the  morphotropic  phase  boundary,  there 
may  expected  to  be  local  competition  between  tetragonal 
and  rbombohedral  polar  shifts,  producing  additional  strain 
contributions  to  foe  free  energy. 

Moreover,  foe  conventional  approach  does  not 
expSchly  take  account  of  foe  details  of  local  c*xm 
displacements  fin  this  system,  ft  would  appear  that  the 
local  symmetry  of  foe  lead  cation  la  tnowocHnic  in  both 
the  rhonfoohcdral  and  tetragonal  phases.  For  foe 
ihombohedral  phase,  it  his  been  proposed  that  the 
displacement  of  the  Pb  catkin  from  their  cubic  position 
comprises  a  homogeneous  [111]  duphoemcot 
superimposed  with  minor,  random  <00 1>  shifts  (13].  In 
foe  tetragonal  phase,  the  homogeneous  (001] 
displacements  arc  augmented  with  minor  <110>  (or 
equivalent  <11 1>)  random  shifts  [14].  Is  both  cases,  foe 
random  shifts  average  to  zero  so  that  foe  macroscopic 
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symmetries  remain  rhombohedral  and  tetragonal 
respectively.  The  random  Pb  displacements  might  be 
considered  as  either  static  shifts  or  dynamic  fluctuations 
about  the  homogeneous  displacement  It  is  not  known 
whether  they  are  distributed  totally  randomly,  or  whether 
their  directions  are  partially  correlated  in  clusters  and 
whether  their  magnitudes  are  equal  on  all  sites. 
Nevertheless,  it  seems  clear  that  die  influence  of  PbTiC>3 
exists  well  into  the  rhombohedral  phase  field  and  that  of 
PbZrC>3  into  die  tetragonal  phase.  Moreover,  in  general, 
the  Ti  and  Zr  ions  do  not  have  equal  displacements  in 
either  phase  [13  &  14]. 

The  minor  displacements  can  be  assumed  to  be  of 
some  significance  in  the  appearance  of  the  monoclinic 
phase.  For  example,  it  might  be  postulated  that  the 
transition  from  the  tetragonal  to  monoclinic  phase  in  PZT 
occurs  due  to  the  condensation  of  the  random 
<111>/<I10>  displacements,  which  combined  with  the 
homogeneous  [001]  shifts  give  rise  to  a  macroscopic 
monoclinic  symmetry.  As  die  Pb  cation  shifts  are  often 
cited  as  being  the  main  component  of  dipole  moment  in 
Pb  perovskite  ferroelectrics,  it  may  be  possible  to 
construct  thermodynamic  expressions  to  reflect  this 
postulate,  in  which  both  the  homogeneous  and  random 


contributions  to  the  polarization  are  represented,  or  in 
which  the  symmetry  of  the  two  end  members  is  allowed  to 
influence  the  stable  state  well  inside  the  solid  solution. 

The  example  explored  here  is  one  in  which  the  free 
energy  is  expressed  as  an  expansion  in  terms  of  two 
polarization  contributions,  representing  the  different 
influences  of  the  two  end  members  of  the  solution  series. 
It  is  based  upon  a  linear  combination  of  the  free  energy 
expansions  of  the  end-members  of  the  solid  solution, 
together  with  linear  and  quadratic  coupling  terms.  The 
total  polarization  is  a  function  of  the  two  end  member 
contributions  and  the  composition  variable. 

THEORY 

In  general,  the  free  energy,  AG,  of  an  infinite, 
homogeneous,  single  domain  ferroelectric  with  a  cubic 
prototype  structure,  as  in  the  perovskites,  is  given  as  an 
expansion  of  the  orthogonal  components  of  the 
polarization  P  in  which  only  the  coefficient  is 
necessarily  temperature  dependent,  with  the  form 
<*2oo  *  <*2oo(T-To): 


ACKa^  P)=  a*,  (P,2  +  P22  +  P,2)  +  a4O0  (P*  +  P24  +  Pi*)  +  am  (P,2  P22  +  P22  P32  +  P32  P,2)  + 

(Pi6  +  P26  +  P36)  +  a«o  Pi4  (P/  +  P,2)  +  P24  (Pi2  +  P.2)  +  P34  (Pi2  +  P22)]  + 

«B*P|  Pi*P»* .  -E,P,-E2P2-E3P3,  (1) 

where  E  is  the  applied  electric  field.  Expansions  encapsulate  the  features  of  solid  solutions  discussed 
terminated  at  the  6*  power  of  polarization  have  generally  above,  a  solid  solution  energy  function  is  proposed  of  the 
been  found  sufficient  to  provide  an  accurate  description  of  form: 
most  ferroelectric  perovskite  systems.  However,  to 


AG*(x)  -  x  GrKoc^p,  p>Kl-x)  GpzOaop,  q^efrn.*,  p,  q) 


where  x  is  the  mole  fraction  of  PbTiCH.  The  vectors  p  and 
q  represent  the  dissimilar  polarizing  tendencies  of  the  two 
end-members,  PbTiD,  and  PbZrCh,  respectively.  The 
forms  of  Gprfcxn^,  p)  and  GpzOmap,  q)  are  identical  to  that 
of  Equation  1,  hence  the  symmetry  requirements  of  the 
cubic  prototype  are  maintained.  The  sets  of  coefficients, 
ctnBp  and  pony,  are  those  of  the  end  members  of  the  solid 
solution,  in  which  only  T0  is  a  function  of  composition 


(2) 

and  is  common  to  both  sets,  in  die  form: 

To(x)  =  x  (TprpTpz)  +  Tp z  (3) 

where  Tpy  and  Tn  represent  die  Curie  temperatures  of 
PbTiOj  and  PbZrD,  respectively. 

The  coupling  energy,  for  terms  up  to  the  4*  power  of 
polarization,  is  of  the  form: 


Gc(r«q»P.<l)=  rioo(Piqi  +  p2q2  +  P3q3)  +  r2oo(pi2qi2  +  p22q22  +  P32q32)  + 

Y220  (Pi2  (qi2  +  qj2)  +  P22  (qs2  +  qi2)  +  Ps2  (qi2  +  qz2))  (4) 


Terms  in  pi%  where  i^j,  are  disallowed  by  symmetry, 
as  the  free  energy  should  be  independent  of  the  sign  of 
piqj.  However,  terms  in  p^  are  necessary  to  differentiate 
between  the  parallel  and  anti-parallel  states  of  pi  and  q,. 
Table  I  lists  the  non-degenerate  states  resulting  from 
combination  of  die  components  of  p  and  q,  using  the 
notation  of  Vanderbilt[12]  for  the  resulting  phases.  Due  to 
the  coupling  energy  Ge,  the  components  of  p  and  q  are 
not  constrained  to  the  6*  power  solutions  of  the  end- 


member  contributions.  For  directions  in  which  the 
component  of  the  total  polarization,  Pb  is  not  constrained 
by  symmetry  to  be  zero,  either  of  the  contributions,  p*  or 
qi,  is  allowed  to  be  zero  as  long  as  pj  +  qi  ^  0.  It  should  be 
noted  that  all  of  the  phases  identified  by  Vanderbilt  [12] 
for  the  conventional  12th  power  Devonshire  model  are 
possible  solutions  of  the  current  approach,  including  three 
non-degenerate  triclinic  phases. 
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TABLE  I.  Classification  of  phases  for  non-degenerate  combinations  of  pi  and  Only  components  constrained  to 
be  zero  are  shown  as  such.  If  Pi  and  qi  are  both  shown  as  non-zero,  either,  but  not  both,  may  assume  a  zero  value. 


Pi 

Pr 

P? 

Qi 

* 

03 

Phase 

0 

0 

0 

0 

0 

0 

Cubic  -  C 

0 

0 

P3 

0 

0 

03 

Tetragonal  -  T 

0 

P3 

P3 

0 

03 

03 

Orthorhombic  -  O 

P3 

Ps 

P3 

<b 

03 

03 

Rhombohedral  -  R 

Pi 

Pt 

P3 

oi 

qi 

03 

Monoclinic  -  MA 

Pi 

Pi 

-P3 

oi 

qi 

03 

Monoclinic  -  MB 

Pi 

P3 

-P3 

qi 

<b 

<b 

Triclinic  -  TV 

Pi 

“P3 

-P3 

qi 

03 

03 

Monoclinic  -  MA 

0 

Pi 

P3 

0 

02 

03 

Monoclinic  -  Me 

0 

P2 

-P3 

0 

02 

03 

Monoclinic  -  Me 

Pi 

P2 

P3 

qi 

02 

03 

Triclinic  -  Tr 

_ El... 

~rP?  . 

Jl 

-02 

—92 _ 

Triclinic  -  Tr 

For  comparison  with  quantitative  available  data,  h  is 
necessary  to  evaluate  the  total  polarization  P,  which  is  a 
function  of  p  and  q.  For  this  purpose,  it  is  assumed  here 
that  the  relationship  is  of  the  simple  form: 

P  *  kp(x)  p  +  kq(x)  q.  (5) 

The  simplest  assumptions  are  that  either  (i)  kp  * 
or  (ii)  kp  =  x  and  =*  1-x.  The  first  assumption  implies 
that  p  and  q  are  assigned  homogeneously  across  all  unit 
cells,  independent  of  B-site  occupation,  Le.  the  dipole 
moment  of  all  unit  cells  is  equal.  However,  it  is  clear  from 
the  structure  data  of  Corker  [13]  and  Frantti  [14]  that  this 
is  unlikely  as  the  Zr  and  Ti  displacements  are  often 
different  within  a  single  composition.  Assumption  (ii) 
implies  that  p  and  q  are  distributed  according  to  the  B-site 
occupancy  of  the  unit  cell.  This  would  mean  that  in  the 
example  addressed  below,  die  dipole  associated  with,  say, 
Zr  cells  in  the  tetragonal  phase,  would  be  zero.  This  is 
also  not  supported  by  the  structure  data  [13  &  14].  The 
evaluation  of  kp  and  kq  will  therefore  be  undertaken  by 
comparison  with  experimental  data. 

CALCULATIONS  FOR  PZT 

The  model  presented  here  is  aimed  at  the  correct 
prediction  of  the  morphotropic  phase  boundary  and  die 
appearance  of  perovskite  ferroelectric  states  with  non¬ 
equal  polarization  components  i.e.  monoclinic  and 
triclinic  phases.  Hence,  the  orthorhombic  anti  ferroelectric 
phase  in  PZT  is  ignored,  and  for  the  purpose  of  these 
calculations,  lead  zirconate  is  assumed  to  be  a 
ihombohedral  ferroelectric.  Similarly,  die  octahedral  tilt 
transition,  which  differentiates  die  high  and  low 
temperature  rhombobedrai  phases  of  PZT,  is  not 
addressed  here. 

Calculations  of  the  temperature  and  composition 
dependence  of  p  and  q  and  hence  the  phase  boundaries 
were  carried  out  using  die  above  model  for  zero  applied 
field.  Choosing  the  A  end-member  Equation  2  to  be 
PbTiOj,  the  solution  for  x  *  1  should  yield,  for  T  <  Ta, 
Pi  =  P2  *  ft*  p3  /  0,  q  ”  0,  Le.  a  tetragonal  ferroelectric. 


whereas  x  =  0  should  yield  solutions  for  a  rhombohedral 
ferroelectric,  i.e.  for  T  <  Tb,  p=<),  qi  -  q2  *  q3 1-  0.  The 
drop  coefficients  for  PbTiCh  used  for  these  calculations 
are  those  given  by  Haun  [4]  and  are  shown  in  Table  II. 
The  pm,,  coefficients  for  a  hypothetically  ferroelectric 
PbZri>3  have  been  estimated  by  extrapolating  the  Haun 
coefficients  for  compositions  with  x  >  0  back  to  x  =  0.  In 
some  cases,  notably  P200  and  P400,  the  extrapolated  value  is 
identical  to  that  of  the  original  Haun  coefficient,  whereas 
the  remaining  coefficients  tend  to  a  linear  dependence  on 
x  except  close  to  x  =  0  and  x  =  0.5. 

The  calculations  were  carried  out  by  numerical 
minimization  of  AG^x)  with  respect  to  p  and  q  using  the 
FindMinimum  function  of  Mathematica  4.0.  Only 
positive  values  of  pi  and  %  were  considered,  hence  in  this 
context  the  y100  term  is  redundant  and  for  simplicity  its 
value  was  set  to  zero,  leaving  only  Y200  and  7220  as  foe 
unknown  parameters. 

The  problems  of  identifying  foe  absolute  free  energy 
minimum  in  six  dimensional  polarization  space  are  not 
trivial,  hence  solutions  were  sought  which  corresponded 
to  foe  most  likely  phases  in  Table  II.  For  each  temperature 
and  PbTi(>3  molar  fraction,  the  local  minima  in  AG«(x) 
were  systematically  identified  for  pi  >  0  and  qi  >  0,  with 
foe  constraints  shown  in  Table  I,  for  foe  following  phases: 
cubic,  tetragonal,  orthorhombic,  rhombohedral, 
monoclinic  MA  and  monoclinic  Me.  The  energies  of  the 
minima  were  then  compared  to  define  the  stable  state. 


TABLE  II.  Free  energy  coefficients  for  PbTiOj  (a^) 
and  hypothetically  ferroelectric  PbZrCh  (p^) 


mnp 

Qtimp 

Poop 

200 

3.74  x  105  (T  -  To) 

2.82  x  10* (T -  T0) 

400 

-7.9  x  107 

5.12x10* 

220 

7.5  x  10* 

-6.5x10* 

600 

2.61  x  10* 

5.93  x  10* 

420 

63  x  10* 

2x10* 

222 

-3.66x10* 

-9.5  x  10’ 
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Lattice  parameters  and  permittivity  were  calculated  for 
the  case  of  PbZro.52Tio.4aO3.  The  spontaneous  strain 
matrix,  o,  is  calculated  according  to: 


<Jjj  “  Qijkl  PkP»  (6) 

For  reasons  of  symmetry,  there  are  only  three 
independent  components  of  fee  electrostriction  coefficient 
tensor  Q,  which  in  reduced  notation  are  Qu,  QI2  and  Q44. 
Their  values  for  PbZr0j2Tio4i03  were  derived  by  linear 
interpolation  from  those  of  fee  end-members,  published 
by  Haun  [4]:  Qu  *  0.0667,  QI2  *  >0.0197  and  Q44  = 
0.0567.  These  were  used  in  preference  to  values 
interpolated  from  die  full  set  of  coefficients  given  at 
intervals  of  Ax  -  0.1  in  fee  same  reference.  The  adopted 
procedure  assumes  that  the  electrostriction  coefficients  are 
“well-behaved”  on  passing  through  the  morphotropic 
phase  boundary  with  no  anomalous  values. 

The  spontaneous  strains  are  converted  to  the  lattice 
parameters  referenced  to  the  pseudo-cubic  axes,  again 
using  reduced  notation,  as  follows: 

ap*bp  =  ao(oi*f  1),  Cp  -  Sq  (o3  +  1),  (7) 

yp  *  n/2  -  tan  *04,  -  nJ2  +  tan* *05 

where  ap,  bp,and  Cp  are  fee  unit  cell  lengths,  yp  is  the  angle 
between  the  a-  and  b-axes,  and  pp  is  the  angle  between  fee 
c-axis  and  fee  x-y  plane,  ao  is  fee  equivalent  cubic  lattice 
parameter  extrapolated  for  the  full  temperature  range  from 
the  data  of  Nobeda  [8]  in  the  paraelectric  phase,  giving: 

ao  -  4.0584  +  2.9156  x  10"5  T  A.  (8) 


For  the  monoclinic  phase  (T  <  300  K),  fee 
transformation  to  the  monoclinic  lattice  parameters 
(denoted  by  subscript  m)  is  given  by: 


•/iFF” f*0)- 

(rX)  c-=c.‘ 

=  arccos|^cos(pp  )j 


a.  =2a 
b_  =  a.  tan] 


(9) 


The  small  signal  permittivities,  s33  and  en,  were 
calculated  from  the  increment  in  the  polarization,  Pit  on 
applying  a  field  E*  of  1  kV  m*1,  for  i  =  3  and  1 
respectively. 


RESULTS 

The  stabilities  of  fee  phases  were  determined  as  a 
function  of  Y200  and  7220-  In  fee  examples  discussed,  only 
stable  phases  corresponding  to  cubic,  tetragonal, 
rhombohedral  and  monoclinic  MA  were  found.  Phase 
diagrams  as  a  function  of  x  and  T  are  shown  for  various 
combinations  of  Y200  anti  Y220  in  Figs.  2, 3,  and  4. 


Fig.  2.  Phase  diagram  for  Y200  "  Y220  °°- 

For  Y200  =  7220  ^  0,  there  is  only  a  single  phase 
boundary  at  To(x),  between  cubic  and  monoclinic  Ma,  for 
all  PbTi03  concentrations  except  the  end-members.  For 
non-zero  values,  Y200  and  Y220  must  exceed  a  threshold 
value  to  avoid  discontinuities  as  x  — >  0  and  x  — ►  1. 

When  both  Y200  and  y^  exceed  2  x  10*  fee  appearance 
of  fee  monoclinic  phase  is  completely  inhibited  see  (Fig. 
2).  The  locus  of  the  morphotropic  boundary  is 
independent  of  7200  and  Y220  and*  well  below  To,  it  is  also 
almost  independent  of  temperature,  moving  to  slightly 
higher  values  of  x  with  decreasing  temperature,  and 
extrapolating  to  x  -  0.45  at  0  K.  On  approaching  To,  the 
boundary  curves  sharply  towards  lower  values  of  x;  the 
point  at  which  fee  morphotropic  phase  boundary  finally 
coincides  with  fee  cubic  phase  boundary  is  at  x  *  0.  Hence 
a  cubic-tetragonal  transition  occurs  on  cooling  for  all 
compositions  wife  x  >  0,  however  for  compositions  up  to 

x  ~  0.35  fee  tetragonal  phase  is  only  stable  for  up  to  5  K 
below  To,  before  transforming  to  the  rhombohedral  phase. 
The  paraelectric-ferroelectric  phase  transition  is  l*1  order 
for  all  values  of  x. 
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Fig.  4.  Phase  diagram  for  Y200  “  6  x  10s,  Y220  *  10*. 

For  intermediate  values  of  the  coupling  coefficients, 
the  Ma  phase  appears  in  the  centre  of  the  phase  diagram. 
On  increasing  the  y  values  from  zero,  there  is  an  almost 
symmetrical  reduction  in  the  width  of  the  monoclinic 
phase  field,  yielding  to  the  rhombohedral  phase  at  low  x 
and  die  tetragonal  phase  at  high  x  values.  There  is  also  a 
reduction  in  die  maximum  temperature  of  monoclinic 
stability  which  occurs  close  to  x  *  0.5.  Thus  as  the 
coupling  energy  is  increased,  the  morphotropic  phase 
boundary  between  the  rhombohedral  and  tetragonal 
phases  is  revealed,  with  the  MA  phase  appearing  as  a 
wedge,  with  its  apex  close  to  die  morphotropic  phase 
boundary.  The  example  shown  in  Figure  3  is  for  Y200  “  Y220 
*  10s,  whilst  that  of  Figure  4  is  for  y^oo  =6x10*  and 
Y220  -  10s,  where  the  inequality  in  the  coupling 
coefficients  results  in  the  monoclinic  phase  favouring  die 
‘‘tetragonal”  side  of  the  morphotropic  phase  boundary. 

The  lattice  parameters  were  calculated  for 
PbZro.32Tio.4sO3.  In  keeping  with  Figure  4  the  value  of  Y200 
was  fixed  at  6  x  10*  to  ensure  that  die  MA  phase  existed 
only  on  the  Ti-rich  side  of  the  MPB.  The  value  of  Y220  was 
adjusted  so  that  the  transition  temperature  for  the 
tetragonal  monoclinic  transition  was  300  K,  as  in  the  data 
of  Noheda  [8].  The  values  of  kp  and  k,  in  Equation  5  were 
than  adjusted  to  fit  the  magnitude  of  die  calculated  lattice 
parameters  to  die  measured  values  [8].  An  acceptable  fit 
was  obtained  for  kp  *  0.73,  k,  -  0.77.  A  comparison  of  the 
calculated  and  measured  values  is  shown  in  Figure  5. 

The  above  values  of  kp  and  k,  where  then  used  to 
calculate  the  relative  permittivities,  833  and  En  (Figure  6). 

DISCUSSION  AND  CONCLUSIONS 

A  novel  thermodynamic  model  has  been  proposed  to 
describe  the  properties  of  Pb(Zr,Ti)03  solid  solutions, 
particularly  in  the  region  of  the  morphotropic  phase 
boundary.  Based  on  an  expression  which  sums  the  two 
end  member  contributions  to  the  free  energy  in  proportion 
to  their  concentration,  the  model  predicts  the  position  of 
the  morphotropic  phase  boundary  at  x  *=  0.45.  This  should 


Fig.  5.  Comparison  of  calculated  and  measured  lattice 
parameters  for  PbZro.52Tia.aO3. 


be  compared  to  the  experimental  locus  of  x  *  0.48.  The 
morphotropic  phase  boundary  is  independent  of  any  of  the 
adjustable  parameters  in  the  model  and  hence  the  good 
agreement  with  experimental  data  appears  to  confirm  the 
suitability  of  die  approach  taken.  The  range  of  existence 
of  the  monoclinic  phase  is  controlled  by  the  coupling 
energy  between  the  two  end-member  contributions  to  the 
polarization,  p  and  q.  At  zero-field,  p  and  q  coexist  only 
in  the  monoclinic  phase.  Although  the  model  allows  for 
the  appearance  of  the  three  possible  perovskite  monoclinic 
phases,  only  the  MA  phase  was  predicted  by  die  model  for 
PZT,  again  in  keeping  with  experiment  In  the  examples 
here,  only  two  coefficients  were  used  to  describe  die 
coupling  energy.  By  increasing  the  coupling  between  die 
parallel  components  of  p  and  q,  whilst  keeping  the 
orthogonal  coupling  energy  at  moderate  levels,  the  MA 
phase  is  restricted  to  the  PbTKVrich  side  of  die 
morphotropic  phase  boundary,  as  suggested  by 
experimental  data  [15]. 

The  comparison  of  the  predicted  and  measured  lattice 
parameters  for  the  specific  case  of  PbZro.52Tio.aO3  is 
generally  good.  Parameters  within  the  model  have  been 
adjusted  to  achieve  a  tetragonal  monoclinic  transition 
temperature  of  300  K  (Y220)  and  to  scale  die  data  correctly 
(kp  and  k^).  Clearly,  there  is  some  discrepancy  between 
the  detailed  behaviour  below  300  K,  which  is  most 
apparent  in  the  c*  parameter. 
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Fig.  6.  Relative  permittivities,  e33  and  En,  for 
PbZr0  52Tio  4*03  calculated  as  function  of  temperature. 

The  calculated  permittivity  for  the  same  composition 
shows  a  small  anomaly  in  G33  at  the  tetragonal  monoclinic 
transition,  but  a  peak  of  over  10,000  is  apparent  far  8n. 
This  is  not  too  surprising  given  that  in  the  monoclinic 

phase,  die  P3  component  is  only  slightly  different  from  die 
saturated  component  in  die  tetragonal  form,  whereas  there 
is  for  greater  compliance  in  die  minor  components,  P]  and 

P2- 

A  feature  of  the  current  model  is  that  die  paraelectric- 
ferroelectric  transition  is  always  predicted  to  be  first  order 
and  the  transition  always  passes  through  the  tetragonal 
phase,  independent  of  composition.  When  To-T  is  small, 
the  a*oo  and  f3<oo  terms  dominate  Gn  and  G* z  respectively, 
hence  x  Gn<T-»To)  will  always  be  less  than 
(l-x)Gpz(T~>To),  even  for  small  x.  This  aspect  is  not 
reflected  in  die  experimental  data  where  die  morphotropic 
phase  boundary  exhibits  minimal  curvature  towards  lower 
x  with  increasing  temperature  and  the  paraelectric- 
ferroelectric  transition  has  only  been  confirmed  as  first 
order  for  compositions  in  which  x  is  substantially  greater 
than  0.5. 

A  second  difficulty  concerns  the  significance  of  p  and 
q,  Here  they  are  interpreted  as  the  contributions  to  die 
macroscopic  polarization,  P,  from  die  two  end  members 
of  the  solid  solution.  However,  as  outlined  above,  die 
mathematical  relationship  between  p,  q  and  P  is  not 


obvious.  A  simple  relationship  has  been  assumed  here 
(Equation  5),  which  on  comparison  of  the  calculated 
lattice  parameters  with  experimental  data,  shows  that 
around  x  *  0.5,  P  »  0.75  (p  +  q).  However,  it  must  be 
assumed  that  for  x  **  0,  P  *  q  and  for  x  -  1,  P  «  p.  Hence, 
the  form  of  k/x)  and  k^x)  is  unlikely  to  be  linear  in  x. 
Both  this  problem  and  that  of  the  first-order  nature  of  the 
paraelectric-ferroelectric  phase  transition  might  be 
alleviated  by  taking  an  approach  that  is  more  directly 
analogous  to  the  oh  initio  method  of  Bellaiche  [16],  in 
which  the  Hamiltonian  includes  a  term  for  the  average 
structure  and  a  summation  of  terms  representing  die  local 
environment 

The  above  difficulties  clearly  reflect  the .  rather 
empirical  nature  of  the  proposed  model.  However,  bearing 
in  mind  that  the  true  structure  of  PZT  does  not  yield  to  a 
conventional  treatment,  based  on  a  rigorous  group 
theoretical  approach,  the  alternative  model  perhaps 
provides  some  directions  towards  a  more  comprehensive 
model. 
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A  two  order  parameter  thermodynamic  model  is  proposed  which  accounts  for  a  number  of  features  of  the  solid  solution  system 
PbCZri-jtTixK^  not  predicted  by  the  conventional  6th  order  model.  The  free  energy  comprises  three  sets  of  terms:  two 
expansions  of  polarization,  representing  contributions  to  the  polarization  of  each  of  the  solid  solution  end-members,  weighted 
linearly  with  respect  to  composition,  plus  terms  representing  the  coupling  between  the  two  contributions.  Employing 
previously  published  coefficients  for  PbTi03  and  values  extrapolated  from  Zr-rich  PZT  compositions  for  PbZi03,  calculations 
predict  that  the  composition  of  the  morphotropic  phase  boundary  between  the  rhombohedral  and  tetragonal  phases  lies  at 
jc  =  0.45  and  is  independent  of  the  strength  of  coupling  between  the  two  order  parameters.  The  model  allows  for  the  existence 
of  a  monoclinic  phase,  the  composition  range  of  which  depends  upon  the  strength  of  the  coupling  terms.  A  satisfactory  fit  to 
the  lattice  parameters  of  Pb(Zro.52Tio.48)03  is  obtained  over  a  wide  range  of  temperature.  [DOT.  10.1143/JJAP.42.7418] 

KEYWORDS:  ferroelectricity,  PZT,  thermodynamic  theory 


1.  Introduction 

Lead  zirconate-titanate,!)  PtyZrj-xTyC^  or  PZT,  is  a 
solid  solution  between  the  two  compounds  lead  titanate 
(PbTiCy  and  lead  zirconate  (PbZrOa).  Both  are  of  the 
perovskite  structure  at  room  temperature,  with  PbTiC>3  being 
tetragonal  and  ferroelectric,  whilst  PbZr03  is  of  orthorhom¬ 
bic  symmetry  and  is  anti  ferroelectric.  The  compounds  form 
a  complete  perovskite  solid  solution  (Fig.  1),  with  no 
discernible  ordering  of  the  Ti  and  Zr  on  the  B-site.  The 
system  has  been  widely  studied  due  to  the  outstanding 
piezoelectric  properties  of  compositions  close  to  jc  =  0.5. 
The  peak  in  the  piezoelectric  coefficient  for  these  compo¬ 
sitions  has  often  been  associated  with  the  presence,  at 
x  =  0.48,  of  a  morphotropic  phase  boundary  (MPB)  between 
the  zirconium-rich  rhombohedral  perovskite  and  the  titani¬ 
um-rich  tetragonal  perovskite  phases. 

Haun^  successfully  modelled  the  PZT  phase  diagram 
using  the  thermodynamic  approach  of  the  Landau-Devon- 
shire  model^  in  which  the  free  energy  associated  with  the 
polar  order  parameter  is  expanded  in  terms  of  a  power  series 
of  the  pseudo-cubic  components  of  the  polarization.  Haun 
determined  the  coefficients  of  the  free-energy  expansion 


Fig.  1,  Revised  phase  diagram  of  PZT,  by  combining  those  of  Jaffe1*  and 
Noheda.21 


experimentally  for  each  of  1 1  discrete  compositions  across 
the  phase  diagram  from  the  temperature  dependence  of  the 
first  and  higher  order  dielectric  stiffnesses.  Determination  of 
the  electrostriction  coefficients  allowed  a  calculation  of  the 
intrinsic  piezoelectric  coefficients,  which  were  shown  to 
peak  at  the  morphotropic  phase  boundary  in  keeping  with 
experiment.  The  locus  of  the  morphotropic  phase  boundary 
was  an  input  to,  rather  than  output  of,  the  treatment,  which  is 
not  predictive  in  terms  of  phase  stability  as  function  of 
composition. 

The  Pb(Zni/3Nb2/3)03-PbTi03  system,  which  exhibits  a 
similar  morphotropic  phase  boundary  at  around 
O^PtyZn^N^yOs-O.lPbTiOs,  has  also  been  modelled 
with  some  success  by  assuming  a  linear  mixing  of  the  free 
energy  coefficients  of  the  end-members  according  to 
composition.8’^ 

Within  recent  years,  three  factors  have  prompted  a 
reassessment  of  the  mechanisms  of  piezoelectricity  in 
PZT;  these  are:  (i)  the  detection  of  a  monoclinic  phase  at 
the  morphotropic  phase  boundary; i0)  (ii)  experimental  and 
theoretical  evidence  for  local  structures  in  the  rhombohedral 
and  tetragonal  phases  in  which  the  polar  displacement  of  the 
lead  ion  comprises  a  homogeneous  component  correspond¬ 
ing  to  the  macroscopic  symmetry  and  local  random  “off- 
axis”  components  which  average  to  zero;1 1-1 3)  and  (iii)  the 
splitting  of  low-lying  transverse  optic  phonon  modes  as  a 
function  of  the  solid  solution  composition. 1415)  The  existing 
thermodynamic  theory  for  PZT  is  not  consistent  with  these 
observations. 

The  monoclinic  phase  in  PZT  is  one  of  the  few  phases  in 
ferroelectric  perovskites  in  which  the  non-zero  components 
of  the  polarization  with  respect  to  the  cubic  axes  are  not 
equal  under  conditions  of  zero  applied  field  and  stress.  The 
most  common  ferroelectric  phases  in  perovskites,  as  exem¬ 
plified  by  barium  titanate,  are  tetragonal,  orthorhombic  and 
rhombohedral,  with  the  polarization  vector  parallel  to  [001], 
[Oil]  and  [111]  respectively.  The  orthogonal  components  of 
polarization  for  these  phases  have  values  that  are  either  zero 
or  equal  to  each  other  (see  Table  I).  These  phases,  plus  the 
paraelectric  (cubic)  phase,  are  the  only  solutions  of  the 
conventional  Devonshire  model  expanded  up  to  the  6th 
power  of  polarization,  which  has  generally  proved  sufficient 
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Table  I.  Possible  symmetries  in  perovskite  fenoelectrics;  phases  above 
the  dotted  line  are  solutions  to  the  6th  power  Devonshire  theory. 


Reference  Crystal  Class 

Symmetry 

Polarization  state 

Polar 

Direction 

C 

Cubic 

m3m 

11 

^0 

11 

& 

H 

0 

paraelcctric 

T 

Tetragonal 

4mm 

F,=/>2=0,P3*0 

(001) 

O 

Orthorhombic 

mm2 

P2  =  0,  i>,2  =  Pi2  *  0 

(101) 

R 

Rhombohedral 

3m 

Pi2=P22  =  Pi2*0 

(111) 

ma 

Monoclinic 

m 

Pi2=p22*o,p}*o 

Pi2  >  Pt2 

Mb 

Monoclinic 

m 

Pi2  =  P22  *0,P}*0 
Pi2  <  Pi 2 

Me 

Monoclinic 

m 

Pi2*P22*0,P3=0 

Tr 

Triclinic 

1 

Pi1*  Pi2*  Pi2*  0 

for  a  satisfactory  fit  to  experimental  data.  In  the  monoclinic 
phase  of  PZT,  the  pseudo-cubic  components  of  polarization 
along  [100]  and  [010]  are  non-zero  and  equal,  but  are 
somewhat  smaller  than  that  along  [001].  Hence,  the  polar¬ 
ization  vector  lies  in  the  (110)  plane  at  an  arbitrary  angle  to 
the  [001]  direction.  The  transition  from  tetragonal  to 
monoclinic  is  of  2nd  order,  with  the  angle  of  the  polarization 
deviation  from  the  [001]  axis  increasing  continuously  from 
zero.  The  polarization  in  the  monoclinic  phase  might 
therefore  be  considered  as  the  homogeneous  superposition 
of  a  small  ihombohedral  component  upon  the  existing 
tetragonal  one. 

Although  monoclinic  phases  are  not  solutions  of  the  6th 
power  Devonshire  theory,  Vanderbilt16)  has  shown  that  the 
inclusion  of  8th  power  terms  in  the  perovskite  free  energy 
expansion  does  lead  to  the  possibility  of  monoclinic  phases. 
However,  the  form  of  the  8th  power  equations  does  not 
account  for  the  observed  splitting  of  the  soft-mode  nor  does 
it  explicitly  address  the  disordered  nature  of  PZT  as  does,  for 
example,  the  first-principles  model  of  Bellaiche.17* 

The  question  addressed  by  this  paper  is  as  follows:  is  it 
possible  to  construct  a  thermodynamic  theory  for  PZT,  (i)  in 
which  the  composition  dependence  of  the  free  energy  is 
implicit,  with  reference  only  to  the  end-members,  (ii)  which 
is  consistent  with  the  observation  of  soft  mode  splitting  and 
(ii i)  which  can  account  for  the  monoclinic  phase  as  a 
superposition  of  the  tetragonal  and  ihombohedral  distor¬ 
tions? 

The  first  principles  model  of  Grinberg13*  supports  the  view 
that  Pb  ion  shifts  are  the  main  contribution  to  the  macro¬ 
scopic  polarization  in  PZT.  The  local  deviations  of  the  Pb 
ion  displacements  from  the  macroscopic  polar  direction  are 
shown  to  encompass  a  range  in  magnitude  and  direction 
which  are  dependent  upon  the  local  Zr  and  Ti  configurations. 
However,  as  an  approximation,  the  spectrum  could  be 
reduced  to  two  types  of  perturbation,  depending  on  whether 
the  local  environment  is  Ti-  or  Zr-rich,  i.e.  predominantly 
PbTiC>3  or  PbZr03.  Hence  the  thermodynamic  model 
proposed  here  is  based  upon  two  polarization  order- 
parameters.  Each  parameter  is  related  to  the  polarization 
contribution  of  one  of  the  two  end-members  of  the  solution 
series.  The  model  comprises  a  summation  of  the  free  energy 
expansions  of  the  end-members  of  the  solid  solution, 
weighted  with  respect  to  composition,  together  with  linear 


and  quadratic  coupling  between  the  two  polarization  terms. 

It  is  shown  here  that  this  model  results  in  an  almost 
temperature  independent  rhombohedral-tetragonal  phase 
boundary,  close  to  the  composition  seen  experimentally.1* 
The  range  of  existence  of  the  neighbouring  monoclinic 
phase  is  dependent  upon  the  strength  of  coupling  between 
the  two  order  parameters. 

2.  Theory 

In  the  conventional  Landau-Devonshire  model,7*  the 
elastic  Gibbs  free  energy,  AG,  of  an  infinite,  homogeneous, 
single  domain  ferroelectric  with  a  cubic  prototype  structure, 
as  in  the  perovskites,  is  given  as  an  expansion  of  the 
orthogonal  components  of  the  polarization,  P,  in  which  only 
the  coefficient  of2oo  is  necessarily  temperature  dependent, 
with  the  form  <2200  —  otf20o(T  —  To): 

A Gict^  P)  =  a2oo(Pi2  +  Pi 2  4-  P32) 

+  <*400(P  l4  +  PlA  +  P34) 

+  &220(P l2 Pi2  +  P2P2  4“  P32Pl2) 

4-  a<50oCPi6  4-  P26  4-  P36) 

4-  ctmlPiW  4-  Pi2)  4-  P2V32  4-  Pi2) 
4-P34(Pi24-P22)] 

4-  dniP^P^P2 . 

—  E\P\  —  E2P2  “  P3P3,  (1) 

where  E  is  the  applied  electric  field.  To  encapsulate  the 
features  of  solid  solutions  discussed  above,  an  energy 
function  for  the  solution  between  end-members  A  and  B  is 
proposed,  of  the  form: 

AGss(x)  =  xGpX&mnpi  p)  4*  (1  *—  x)(jTQ^nmj}^  q) 

4 -Gc(Ynwp>P>q) 

where  x  is  the  mole  fraction  of  end-member  A.  The  forms  of 
GpSctmnp, p)  and  GB(Pmnp^q)  are  identical  to  that  of  eq.  (1), 
hence  the  symmetry  requirements  of  the  cubic  prototype  are 
maintained.  The  sets  of  coefficients,  and  are  those 
of  the  end  members  of  the  solid  solution,  in  which  only  To  is 
a  function  of  composition  and  is  common  to  both  sets,  in  the 
form: 

T0(x)  =  x(Ta  -  7b)  4*  Tb  (3) 

where  TA  and  Tb  represent  the  Curie  temperatures  of  the  two 
end-members. 

The  coupling  energy,  for  terms  up  to  the  4th  power  of 
polarization,  is  of  the  form: 

Gc(Ymnp>p,q)  =  Ymipiqi  4-  piqi  4-  piqi) 

4-  K20o(Pi2<7i2  4-  p2q2  4-  p2q2) 

4-  X22o(pi2(<?22  4-  q2)  4-  Pi(q2  4-  q2) 

+  P32(?i2  +  q-i))  (4) 

This  is  consistent  with  the  treatment  by  Salje18)  for 
ferroelastic  materials  with  two  order-parameters,  in  which 
bi-linear  and  bi-quadratic  coupling  terms  appear  to  be  the 
most  important.  Higher  order  coupling  terms  are  only 
excluded  here  for  a  reduction  in  complexity.  However, 
terms  in  piqj,  where  i  ^  j\  are  disallowed  by  symmetry,  as 
the  free  energy  should  be  independent  of  the  sign  of  piq^  but 
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Table  IL  Classification  of  phases  for  non-degenerate  combinations  of  p\ 
and  q\.  Only  components  constrained  to  be  zero  are  shown  as  such.  If  p, 
and  qi  are  both  shown  as  non-zero,  either,  but  not  both,  may  assume  a 
zero  value. 


Pi 

pi. 

Pi 

4\ 

4i 

4i 

Phase 

0 

0 

0 

0 

0 

0 

Cubic  -  C 

0 

0 

Pi 

0 

0 

4i 

Tetragonal  -  T 

0 

Pi 

Pi 

0 

4i 

4i 

Orthorhombic  -  O 

Pi 

Pi 

Pi 

4i 

4i 

4i 

Rhombohedral  -  R 

Pi 

Pi 

Pi 

qi 

4\ 

4i 

Monoclinic  -  Ma 

Pi 

Pi 

-Pi 

41 

4i 

4i 

Monoclinic  -  MB 

Pi 

Pi 

-Pi 

41 

4i 

4i 

Triclinic  -  Tr 

Pi 

- Pi 

-Pi 

4i 

4i 

4i 

Monoclinic  -  MA 

0 

Pi 

Pi 

0 

4i 

4i 

Monoclinic  -  Me 

0 

Pi 

- Pi 

0 

4i 

4i 

Monoclinic  -  Me 

Pi 

Pi 

Pi 

4i 

4i 

4i 

Triclinic  -  Tr 

Pi 

Pi 

-Pi 

4\ 

41 

4i 

Triclinic  -  Tr 

terms  in  piqt  are  necessary  to  differentiate  between  the 
parallel  and  anti-parallel  states  of  pt  and  qi. 

Table  II  lists  the  non-degenerate  states  resulting  from 
combination  of  the  components  of  p  and  q,  using  the 
notation  of  Vanderbilt16^  for  the  resulting  phases.  Due  to  the 
presence  of  the  coupling  energy  Gc,  the  components  of  p 
and  q  are  not  constrained  to  the  6th  power  solutions  of  the 
end-member  contributions.  For  directions  in  which  the 
component  of  the  total  polarization,  P,,  is  not  constrained 
by  symmetry  to  be  zero,  either  of  the  contributions,  pi  or  qx, 
is  allowed  to  be  zero  as  long  as  pi  +  q\  0.  It  should  be 
noted  that  all  of  the  phases  identified  by  Vanderbilt16)  for  the 
conventional  12th  power  Devonshire  model  are  possible 
solutions  of  the  current  approach,  including  three  non- 
degenerate  triclinic  phases. 

The  relationship  between  the  total  polarization,  P,  and  p 
and  q  is  not  explicit  in  the  model,  but  is  assumed  here  to  be 
of  the  simple  form: 

P  =  h(x)p  -f  k(x)q ,  (5) 

two  possible  examples  of  which  are  (i)  h{x)  =  k( x)  =  1  and 
(ii)  h(x)  =  x  and  k(x)  =  (1  —  x).  The  first  assumption  implies 
that  p  and  q  are  summed  homogeneously  in  all  unit  cells, 
independent  of  B-site  occupation,  i.e.  the  dipole  moment  of 
all  unit  cells  is  identical.  However,  it  is  clear  from  the 
structure  datalu2)  and  first  principles  calculations 13)  that  this 
is  unlikely  for  PZT.  Assumption  (ii)  implies  that  p  and  q  are 
distributed  according  to  the  B-site  occupancy  of  the  unit  cell. 
This  would  imply  that  the  dipole  associated  with,  say,  Zr- 
rich  regions  in  the  tetragonal  phase,  is  zero.  Again,  this  is  not 
supported  by  the  structure  refinements  and  calculations, 
therefore  the  approach  taken  here  is  to  attempt  to  discern  the 
values  of  h(x)  and  k(x)  by  comparison  with  previous  data  for 
the  system. 

3.  Calculations  for  PZT 

The  model  is  aimed  at  determining  the  states  of  PZT  close 
to  the  morphotropic  phase  boundary,  therefore  the  ortho¬ 
rhombic  antiferroelectric  phase  in  PZT  is  ignored,  and  for 
the  purpose  of  these  calculations,  lead  zirconate  is  assumed 
to  be  a  ihombohedral  ferroelectric.  Similarly,  the  octahedral 


Tabic  m.  Free  energy  coefficients  for  PbTiQj  and  hypothetically 

ferroelectric  PbZrOj 


mnp 

&mnp 

A n«p 

200 

3.74  x  105  (T  —  To) 

2.82  x  105  (T  -  To) 

400 

-7.9  x  107 

5.12  x  108 

220 

7.5  x  108 

-6.5  x  108 

600 

2.61  x  L08 

5.93  x  108 

420 

6.3  x  108 

2x  109 

222 

-3.66  x  109 

-9.5  x  109 

tilt  transition,  which  differentiates  the  high  and  low  temper¬ 
ature  rhombohedral  phases  of  PZT,  is  not  addressed. 

Calculations  of  the  temperature  and  composition  depend¬ 
ence  of  p  and  q  and  hence  the  phase  boundaries  were  carried 
out  using  the  above  model  for  zero  applied  field.  Choosing 
the  ‘A’  end-member  in  eq.  (2)  to  be  PbTi03,  the  solution  for 
x  =  1  should  yield,  for  T  <  rA,  Pi  =  P2  =  0,  p3  ^  0, 
q  =  0,  i.e.  a  tetragonal  ferroelectric,  whereas  x  =  0  should 
yield  solutions  for  a  rhombohedral  ferroelectric,  i.e.  for 
T  <  rB,  p  =  0,  q\  =  qi  —  q$  ^  0.  The  coefficients  for 
PbTi03  used  for  these  calculations  are  those  given  by  Haun4) 
and  are  shown  in  Table  HI.  The  Pmnp  coefficients  for  a 
hypothetically  ferroelectric  PbZr03  have  been  estimated  by 
extrapolating  the  Haun  coefficients  for  compositions  with 
x  >  0  back  to  x  =  0.  In  some  cases,  notably  fhm  and  ^400, 
the  extrapolated  value  is  identical  to  that  of  the  original 
Haun  coefficient,  whereas  the  remaining  coefficients  tend  to 
a  linear  dependence  on  x  except  close  to  x  =  0  and  x  =  0.5. 

The  calculations  were  carried  out  by  numerical  minimi¬ 
zation  of  AGssW  with  respect  to  p  and  q  using  the 
FindMinimum  function  of  Mathematica  4.0.  Only  positive 
values  of  pi  and  qi  were  considered,  hence  in  this  context  the 
yioo  term  is  redundant  and  for  simplicity  its  value  was  set  to 
zero,  leaving  only  y^oo  and  yno  as  the  unknown  coefficients. 

The  problems  of  identifying  the  absolute  free  energy 
minimum  in  six  dimensional  polarization  space  are  not 
trivial,  hence  solutions  were  sought  which  corresponded  to 
the  most  likely  phases  in  Table  II.  For  each  temperature  and 
PbTi03  molar  fraction,  the  local  minima  in  AGss(x)  were 
systematically  identified  for  pf  >  0  and  qi  >  0,  with  the 
constraints  shown  in  Table  I,  for  the  following  phases:  cubic, 
tetragonal,  orthorhombic,  rhombohedral,  monoclinic  MA 
and  monoclinic  Me .  The  energies  of  the  minima  were  then 
compared  to  define  the  stable  state. 

Lattice  parameters  were  calculated  for  the  case  of 
PbZro.52Tio.4sO3.  The  spontaneous  strain  matrix,  <r,  is 
calculated  according  to: 

<*ij  =  QijkiPkPi  (6) 

For  reasons  of  symmetry,  there  are  only  three  independent 
components  of  the  electrostriction  coefficient  tensor  g, 
which  in  reduced  notation  are  gn,  gn  and  £44-  Their  values 
for  PbZro.52Tio.48 03,  Qn  =  0.0667,  Q12  =  —0.0197  and 
Qm  ~  0.0567,  were  derived  by  linear  interpolation  from 
those  of  the  end-members,  published  by  Haun  4)  These  were 
used  in  preference  to  values  interpolated  from  the  full  set  of 
coefficients  in  the  same  reference.  The  adopted  procedure 
assumes  that  the  “true”  electrostriction  coefficients  are 
“well-behaved”  on  passing  through  the  morphotropic  phase 
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boundary  with  no  anomalous  values. 

The  spontaneous  strains  are  converted  to  the  lattice 
parameters  referenced  to  the  pseudo-cubic  axes,  again  using 
reduced  notation,  as  follows: 

dp -bp-  ao(cri  +  1),  cp  =  oo(03  +  *)»  (7) 

Yp  =  7r/2  —  tan“l  <r4,  fip  =  tt/2  -f  tan"1  ^5 

where  ap,  6P, and  cp  are  the  unit  cell  lengths,  yp  is  the  angle 
between  the  a -  and  6-axes,  and  fip  is  the  angle  between  the  0 
axis  and  the  a~b  plane,  oq  is  the  equivalent  cubic  lattice 
parameter  extrapolated  for  the  foil  temperature  range  from 
the  data  of  Noheda10)  in  the  paraelectric  phase,  giving: 

a0  =  4.0584  +  Tx  2.9156  x  10"5  A  (8) 

where  T  is  temperature. 

For  the  monoclinic  phase,  the  transformation  to  the 
monoclinic  lattice  parameters  (denoted  by  subscript  m)  is 
given  by: 

am  =  2aj,j  ^2(1  +  tan2(yp/2)), 

bm  =  am  tan(yp/2),  cm  =  cp  cos(y8p  -  n/ 2),  (9) 

fin  =  arccos|^cos(/3p)  j  cos(j/p/2)j 

4.  Results 

Phase  stabilities  were  determined  as  a  function  of  >^oo  and 
X220.  In  the  examples  discussed,  only  phases  corresponding 
to  cubic,  tetragonal,  rhombohedral  and  monoclinic  MA  were 
found.  Phase  diagrams  as  a  function  of  jc  and  T  are  shown  for 
various  combinations  of  y^oo  and  >*220  in  Figs.  2,  3,  and  4. 

For  Y2oo  =  ym  <  0,  there  is  only  a  single  phase  boundary 
at  7o(jt),  between  the  cubic  and  monoclinic  MA  phases,  for 
all  PbTiC>3  concentrations  except  the  end-members,  y^oo  and 
Yi 20  must  be  in  a  certain  range  to  avoid  discontinuities  as 
x  0  and  x-+  1. 

When  both  K200  and  yi 20  exceed  2  x  108,  the  appearance 
of  the  monoclinic  phase  is  completely  inhibited  (see  Fig.  2). 
The  locus  of  the  morphotropic  boundary  is  independent  of 
K200  and  Y2 20  and,  well  below  7o,  it  is  also  almost 
independent  of  temperature,  moving  to  slightly  higher 
values  of  x  with  decreasing  temperature,  and  extrapolating 


to  x  =  0.45  at  0  K.  On  approaching  7o,  the  boundary  curves 
sharply  towards  lower  values  of  x;  the  point  at  which  the 
morphotropic  phase  boundary  finally  coincides  with  the 
cubic  phase  boundary  is  at  x  =  0.  Hence  a  cubic-tetragonal 
transition  occurs  on  cooling  for  all  compositions  with  jc  >  0, 
however  for  compositions  up  to  jc  =  0.35  the  tetragonal 
phase  is  only  stable  for  up  to  5  K  below  To,  before 
transforming  to  the  rhombohedral  phase.  The  paraelectric- 
ferroelectric  phase  transition  is  1st  order  for  all  values  of  x. 

For  intermediate  values  of  the  coupling  coefficients,  the 
Ma  phase  appears  in  the  centre  of  the  phase  diagram.  On 
increasing  the  y  values  equally  from  zero,  there  is  an  almost 
symmetrical  reduction  in  the  width  of  the  monoclinic  phase 
field,  yielding  to  the  rhombohedral  phase  at  low  jc  and  to  the 
tetragonal  phase  at  high  x  values.  There  is  also  a  reduction  in 
the  maximum  temperature  of  monoclinic  stability,  which 
occurs  close  to  jc  =  0.5.  Thus  as  the  coupling  energy  is 
increased,  the  morphotropic  phase  boundary  between  the 
rhombohedral  and  tetragonal  phases  is  revealed,  with  the  MA 
phase  appearing  as  a  wedge,  with  its  apex  close  to  the 
morphotropic  phase  boundary.  The  example  shown  in  Fig.  3 
is  for  Y2qo  —  Y220  =  108,  whilst  that  of  Fig.  4  is  for  >*200  = 
6  x  108  and  K220  =  108,  where  the  inequality  in  the  coupling 


7422 


Jpo.  J.  Appl.  Phys.  Vol.  42  (2003)  Pt  1,  No.  12 


A.  J.  Bell  and  E.  Furman 


Fig.  5.  Calculated  values  of  py  and  q\  as  a  function  of  composition  at 
250K  for  noo  =  6  x  108,  yrm  =  10s. 


Temperature  /  K 

Fig.  6.  Comparison  of  calculated  and  measured  lattice  parameters  for 
PbZro.52Tio.48Q3. 


coefficients  results  in  the  monoclinic  phase  favouring  the 
“tetragonal”  side  of  the  morphotropic  phase  boundary.  The 
corresponding  values  of  p$  and  q\  as  a  function  of 
composition  at  250  K  are  shown  in  Fig.  5. 

The  lattice  parameters  were  calculated  for  PbZro.52- 
Tio.4s03.  In  keeping  with  Fig.  4  the  value  of  ^ 00  was  fixed 
at  6  x  108  to  ensure  that  the  Ma  phase  existed  only  on  the 
Ti-rich  side  of  the  MPB.  The  value  of  ym  was  adjusted  so 
that  the  tetragonal  to  monoclinic  transition  temperature  was 
300  K,  as  in  the  data  of  Noheda.I9)  The  values  of  h(x)  and 
k(x)  in  eq.  (5)  were  then  adjusted  to  fit  the  magnitude  of  the 
calculated  lattice  parameters  to  the  measured  values.  An 
acceptable  fit  was  obtained  for  =  0.73*  fc(jt)  =  0.77.  A 


comparison  of  the  calculated  and  measured  lattice  param¬ 
eters  is  shown  in  Fig.  6. 

5 .  Discussion 

The  model  presented  can  reproduce  some  aspects  of  the 
PZT  system  with  reasonable  accuracy.  It  predicts  the 
position  of  the  morphotropic  phase  boundary,  independent 
of  any  of  the  adjustable  parameters,  at  x  =  0.45,  whereas  the 
consensus  of  experimental  data  suggests  the  MPB  lies  at 
X  =  0.48. 

The  range  of  existence  of  the  monoclinic  phase  is 
controlled  by  the  coupling  energy  between  the  two  contri¬ 
butions  to  the  polarization,  p  and  q .  At  zero-field,  p  and  q 
coexist  only  in  the  monoclinic  phase.  Although  the  model 
allows  for  the  appearance  of  the  three  possible  perovskite 
monoclinic  phases,  only  the  MA  phase  was  predicted  by  the 
model  for  PZT.  Tn  the  examples  here,  only  two  coefficients 
were  used  to  describe  the  coupling  energy.  By  increasing  the 
positive  coupling  energy  between  the  parallel  components  of 
p  and  q  (i.e.  reducing  the  tolerance  of  the  system  to  large 
values  of  pi  and  q\),  whilst  keeping  the  orthogonal  coupling 
energy  at  moderate  levels,  the  MA  phase  is  restricted  to  the 
PbTiC>3-rich  side  of  the  morphotropic  phase  boundary.  This 
is  illustrated  by  the  example  of  Fig.  4  in  which  the  non-zero 
components  of  p  and  q  in  the  monoclinic  phase  were  ^3,  q\ 
and  q2  (=  q\ ). 

The  values  of  p  and  q  obtained  from  the  free  energy 
minimization  should  be  interpreted  as  macroscopic  average 
values.  In  the  tetragonal  phase,  whilst  p  is  non-zero,  q  is 
zero  and  vice-versa  in  the  rhombohedral  phase.  Only  in  the 
monoclinic  phase  are  they  simultaneously  non-zero.  It  may 
be  postulated  that  in  the  tetragonal  and  rhombohedral  phases 
the  non-zero  parameter  represents  the  polarization  due  to  the 
major,  homogeneous  cation  shifts  whilst  the  parameter  with 
zero  value  represents  the  average  of  the  minor,  random 
shifts.  Thus  the  transition  into  the  monoclinic  state  might  be 
regarded  as  an  order-disorder  transition,  in  which  either  the 
p  or  q  parameter  passes  from  the  disordered  to  the  ordered 
state. 

6.  Conclusions 

A  thermodynamic  model  has  been  proposed  to  describe 
the  properties  of  Pb(Zr,Ti)C>3  solid  solutions,  particularly  in 
the  region  of  the  morphotropic  phase  boundary.  The  model 
is  based  on  an  expression  in  terms  of  two  coupled  polar¬ 
ization  contributions,  p  and  qt  representing  the  preferred 
symmetries  of  compositions  towards  the  two  extremes  of  the 
solid  solution.  The  model  not  only  locates  the  composition 
of  the  morphotropic  phase  boundary  with  acceptable 
accuracy,  but  it  also  predicts  the  presence  of  a  monoclinic 
phase  around  the  morphotropic  phase  boundary  and,  with  a 
limited  number  of  adjustable  parameters,  agrees  well  with 
experiment.  The  comparison  of  the  predicted  and  measured 
lattice  parameters  for  the  specific  case  of  PbZro.52Tio.48O3  is 
generally  good.  However,  the  relationship  between  the 
macroscopic  polarization,  i>,  and  the  parameters  p  and  q  has 
only  been  determined  empirically  for  this  composition. 

Only  in  the  monoclinic  phase  are  both  p  and  q  non-zero; 
in  the  tetragonal  and  rhombohedral  phases  the  equilibrium 
value  of  one  of  the  parameters  is  always  zero.  It  has  been 
postulated  that  in  the  tetragonal  and  rhombohedral  phases 
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the  two  polarization  contributions  correspond  to  the  homo¬ 
genous  and  local  random  components  of  the  cation  displace¬ 
ments. 
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Strain-gradient- induced  polarization  or  flexoelectricity  was  investigated  in  unpoled  soft  lead 
zirconate  titanate  (PZT)  ceramic  where  the  texture  symmetry  °°"m  forbids  macropiezoelectricity. 
Even  under  high  strain  gradient  (Im"1)  the  induced  polarization  is  small  (1.6  fxCfm2)  at  20  °C. 
Higher  strain  gradients  induce  ferroelastic  poling  and  an  additional  extrinsic  contribution  to  the 
fiexoelectric  coefficient  jin  raising  the  value  from  0.5  to  2.0  /xC/m.  Cooling  through  the  Curie  point 
(Tc)  under  maximum  stress  (80  MPa)  where  the  peak  permittivity  (~~20  000)  could  raise  /x12  to  20 
/xC/m,  the  equivalent  electric  field  is  still  only  —  1  kV/m,  inadequate  to  achieve  significant 
ferroelectric  poling.  The  situation  may  be  different  in  thin  PZT  films  where  much  larger  strain 
gradients  can  occur.  ©  2003  American  Institute  of  Physics.  [DOI:  10.1063/1.1570517] 


The  well-known  elasto electric  coupling  effects  include 
piezoelectric  effect,  electrostrictive  effect,  and  Maxwell 
stress  effect.  Piezoelectric  ceramics1  and  more  recently 
single  crystals2  have  demonstrated  the  capability  for  broad 
applications  in  sensor,  actuator,  and  transducer  devices.  In 
soft  polymers,  it  is  shown  recently  that  Maxwell  stress  effect 
(attractive  forces  between  opposite  charges  on  the  elec¬ 
trodes)  can  generate  ultrahigh  strain  responses  and  exhibit 
great  potential  for  a  variety  of  electromechanical  device 
applications.3  All  the  earlier-mentioned  physical  effects, 
however,  generally  assume  the  situations  of  uniform  stress  or 
strain.  In  nature,  there  is  elastoelectric  coupling  caused  by 
inhomogeneous  deformation  where  stress  or  strain  gradient 
associated  polarization  effects  (fiexoelectric  effects)  need  to 
be  considered.  Overall,  mechanical  stress  or  strain  can  gen¬ 
erate  electric  polarization  in  a  deformable  dielectric  material 
through  the  following  two  mechanisms: 

dsjk 

P i~  dijk°‘jkJr  Pijkl  '  0) 

In  Eq.  (1)  and  thereafter  Einstein  summation  convention  is 
assumed  (ijykj-  1,2,3).  The  first  term  on  the  right-hand 
side  refers  to  the  well-known  direct  piezoelectric  effect, 
where  ajk  is  the  stress  uniformly  distributed  across  the 
sample  and  djjk  is  the  piezoelectric  coefficient,  a  third-rank 
polar  tensor.  The  second  term  on  the  right-hand  side  refers  to 
the  strain  gradient  ( dejk!dxi )  induced  polarization  and  fiijki 
is  the  fiexoelectric  coefficient,  a  fourth-rank  polar  tensor.  In 
centro symmetric  materials,  dijk=0,  so  the  piezoelectric  term 
in  Eq.  (1)  can  be  eliminated,  therefore 


pi  Hijki  dx » 


where  Pf  is  the  electric  polarization  induced  solely  by  strain 
gradient. 
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When  reviewing  the  history  of  fiexoelectric  investiga¬ 
tions,  it  is  noted  that  while  the  concept  was  originally  formed 
in  early  1960s,4,5  not  until  1981  was  such  effect  in  crystalline 
solids  given  the  name  “fiexoelectric.”  6  Based  upon  an  ionic 
model,7  Tagantsev  analyzed  the  fiexoelectric  effect  and  sug¬ 
gested  possible  larger  effects  in  ferroelectrics.  In  soft  poly¬ 
mers  Marvan  et  a/.8  observed  fiexoelectric  coefficients  of  the 
order  of  10' u- 10” 10  C/m.  Up  to  now  very  little  attention 
has  been  paid  to  test  the  magnitude  of  fiexoelectric  coeffi¬ 
cients  and  the  mechanism  remains  unclear.  Recently  we  mea¬ 
sured  the  fiexoelectric  coefficients  (/x12)  in  lead  magnesium 
niobate  (PMN)9,10  (a  well-known  relaxor  ferroelectric  mate¬ 
rial)  and  barium  strontium  titanate  (BST) 11,12  (a  normal 
ferroelectric  material).  Both  materials  were  tested  in  the 
phase  region  with  macroscopic  cubic  symmetry,  moreover 
the  samples  were  measured  in  the  form  of  a  cantilevered 
beam  so  that  any  remnant  piezoelectric  contributions  from 
the  top  and  bottom  halves  of  the  beam  would  cancel.  In  this 
letter,  we  investigate  the  fiexoelectric  effect  in  a  well-known 
lead  zirconate  titanate  (PZT)  piezoelectric  ceramic  in  the 
ferroelectric  phase  by  using  a  four-point  bend  fixture  to  gen¬ 
erate  a  uniform  strain  gradient 

In  this  work  several  interesting  questions  were  explored: 
(i)  using  a  four  point  bending  fixture  is  it  possible  to  induce 
and  measure  a  quasistatic  fiexoelectric  polarization  generated 
by  /x12;  (ii)  is  the  fiexoelectric  polarization  enhanced  or  in¬ 
hibited  by  the  onset  of  ferroelastic  domain  wall  motion 
which  will  be  evidenced  by  change  of  the  flexural  stiffness 
and  the  development  of  remnant  curvature  in  the  sample;  and 
(iii)  is  it  possible  in  soft  PZT  to  reach  levels  of  fiexoelectric 
induced  field  sufficient  to  pole  the  ceramic  into  a  piezoelec¬ 
tric  form. 

The  samples  used  were  unpoled  PZT~_5H  ceramics 
(doped  with  La  and  Sn)  fabricated  by  TRS  Ceramics  CoriP 
pany.  State  College,  Pennsylvania.  Dielectric  measurements 
performed  using  an  HP4284A  LCR  meter  show  a  weak  field 
permittivity  of  2200  at  20  °C  and  a  strong  but  rounded  di¬ 
electric  maximum  (£3=«20000)  without  little  dispersion 
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FIG.  2.  Schematic  illustration  of  strain  gradient  induced  quasi-static  polar¬ 
ization  measurement,  (a)  typical  four-point  bending  fixture;  (b)  stress  distri¬ 
bution;  (c)  schematic  of  strain  gradient  and  the  induced  polarization  along 
the  thickness  of  the  sample. 


FIG.  1.  (a)  Weak  field  dielectric  permittivity  (£—10  V/cm)  as  a  function  of 
frequency  and  temperature  in  the  soft  PZT-5H  sample  measured;  (b)  dielec¬ 
tric  hysteresis  (P  vs  E)  in  the  soft  PZT-5H  composition. 

over  the  frequency  range  of  1  kHz-1  MHz  [Fig.  1(a)],  sug¬ 
gesting  a  diffuse  phase  transition  without  strong  relaxor  char¬ 
acter.  Polarization  hysteresis  loops  measured  using  a  modi¬ 
fied  Sawyer-Tower  circuit  show  remnant  polarization  of 
35  fxd cm2  and  coercive  field  of  7  kV/cm  [Fig.  1(b)]  with  no 
discemable  bias.  Young’s  modulus  was  measured  to  be  70 
GPa  by  a  dynamic  resonance  method. 

A  uniform  strain  gradient  was  generated  using  a  four- 
point  bend  fixture  schematically  illustrated  in  Fig.  2(a). 
Samples  for  measurement  had  dimensions  60  mm  long,  7 
mm  wide,  and  3  mm  thick.  Surfaces  were  carefully  polished 
and  the  samples  were  annealed  at  700  °C  to  relieve  surface 
stresses.  Sputtered  gold  electrodes  with  dimensions  10 
X7mm2  were  applied  to  upper  and  lower  surfaces  in  the 
center  of  the  beam.  No  residual  piezoelectricity  could  be 
detected  by  Berlincourt  d 33  meter,  and  the  impedance  trajec¬ 
tory  from  1  kHz  to  1  MHz  was  free  from  evidence  of  piezo¬ 
electric  resonance.  The  bend  test  was  carried  out  according 
to  the  ASTM  C-l  161-94  by  using  an  Instron  machine  model 
4202  with  a  10  kN  load  cell.  Outer  and  inner  spans  of  the 
fixture  are  40  and  20  mm,  respectively.  The  generated  elec¬ 
tric  charge  was  detected  by  a  Keithley  65 17  electrometer  that 
can  resolve  10  fC  and  will  measure  up  to  2.1  /zC.  Before 
measurement  the  electrometer  was  carefully  calibrated  for 
voltage  burden  and  input  offset  current. 

The  stress  distribution  along  the  sample  length  direction 


(jj)  is  shown  in  Fig.  2(b).  Within  the  inner  span  the  stress 
crn(xi)  is  uniform  along  the  length,  while  along  the  thick¬ 
ness  direction  (*3)  the  stress  crn(;t3)  varies  and  there  is  a 
stress  or  strain  gradient  as  shown  by  Fig.  2(c),  which  can 
generate  electric  polarization  through  the  flexoelectric  effect. 
Because  the  bar  length  is  much  greater  than  the  bar  thick¬ 
ness,  we  omit  the  shear  stress  and  strain  and  only  consider 
the  principal  stress  and  strain.  Therefore,  for  simplification 
here  only  one  suffix  was  used  for  describing  the  stress  and 
strain  tensors.  The  absolute  value  of  surface  stress  was  cal¬ 
culated  using  the  following  equation: 

,  3FL 

(3) 

where  F  is  the  load,  w  is  the  width,  d  is  the  thickness,  and  L 
is  the  outer  span  of  the  bend  fixture. 

The  strain  gradient  in  the  thickness  direction  is  given  by 

<?Ei(x3)  list 


where  s  is  the  cross  head  speed  and  t  the  time  gone  by. 

Figure  3  presents  data  on  the  surface  stress  versus  strain 
curves  for  tests  carried  out  at  1  and  0.2  mm/min  crosshead 
speed.  As  shown  in  Fig.  3(a),  the  flexure  strength  is  mea¬ 
sured  to  be  85  MPa.  The  softening  of  the  samples  at  certain 
stress/speed  levels  (e.g.,  around  30  MPa  for  1  mm/min  cross¬ 
head  speed)  we  believe  corresponds  to  the  onset  of  ferroelas- 
tic  domain  motion.  As  expected  in  PZT  the  sample  becomes 
softer  at  lower  loading  rates  due  to  the  relaxational  nature  of 
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FIG.  3.  Surface  stress  as  a  function  of  strain  measured  at  two  different 
crosshead  speeds  on  PZT-5H  bar,  (a)  1  mm/min  crosshead  speed  (the  arrow 
indicates  the  level  of  stress  at  which  the  sample  broke);  (b)  0.2  mm/min 
crosshead  speed. 

ferroelastic  domain  switching.  That  ferroelastic  switching 
has  occurred  at  the  higher  loading  levels  is  evident  by  a  static 
remnant  curvature  in  the  sample  after  testing. 

Flexoelectric  polarization  versus  strain  gradient  is  dis¬ 
played  in  Fig.  4.  It  is  clear  that  the  behavior  is  not  linear, 
showing  a  low-gradient  and  high-gradient  linear  behavior  but 
of  different  slope.  The  changeover  occurs  at  the  gradient 
~  0.3  m” 1  corresponding  to  the  onset  of  ferroelastic  switch¬ 
ing.  At  very  small  strain  level,  the  sample  is  only  subjected 
to  elastic  deformation,  so  the  measured  electric  polarization 
response  represents  the  intrinsic  flexoelectric  effect.  The 
low-gradient  slope  corresponds  to  a  flexoelectric  coefficient 
jjlCI m,  and  the  high-gradient  slope  takes  fj,[2 
^2.0/ttC/m  showing  that  ferroelastic  domain  wall  motion 
aids  the  response. 

The  measured  flexoelectric  polarization  is  very  small 
compared  to  the  remnant  polarization,  only  1.6  /-iC/m2  at  a 
strain  gradient  of  1  m^1.  By  using  a  relative  permittivity 
value  of  er=2200  (at  1  kHz),  we  figure  out  that  a  strain 
gradient  of  1  m”1  is  equivalent  to  an  electric  field  of  100 
V/m,  which  is  obviously  too  small  compared  to  the  coercive 
field.  In  high-permittivity  ferroelectrics  we  may  expect  from 
elementary  theory  of  flexoelectricity7  that  be  propor¬ 
tional  to  dielectric  susceptibility  Xij  following  a  relation:11 

e 

(5) 

where  y  is  a  constant  of  value  close  to  unity.  For  the  PZT 
unpoled  ceramic  at  low  gradient  levels  the  normalized  flexo¬ 
electric  coefficients  ii\2ixn  at  20  °C  is  0.23  nC/m  giving 
y^O.57.  In  PMN  from  our  earlier  study,9,10  fi\2lx 22  ls  0.26 
nC/m  giving  yte0.65,  both  in  reasonable  accord  with  the 
elementary  theory.  For  BST  ceramic  at  25  °C,  however, 
/^,2/^22^3.72  nC/m  yielding  a  value  y=9.3u  much  higher 
than  those  in  the  lead-based  systems. 
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FIG.  4.  Polarization  vs  strain  gradient  for  an  unpoled  PZT-5H  sample  in¬ 
duced  during  a  four-point  bend  test  carried  out  at  a  crosshead  speed  of  0.2 
mm/min. 

It  may  be  noted  that  in  earlier  investigations  of  the  ther¬ 
mopolarization  effects  Strnkov  et  a/.13  found  exceedingly 
high  values  for  the  normalized  thermopolarization  coefficient 
(Z>°-)  in  triglycine  sulfate  which  they  attributed  to  the  order- 
disorder  nature  of  the  ferroelectric  phase  change  in  this  com¬ 
pound.  We  note  that  in  BaxSri_xTi03,13  there  is  strong  evi¬ 
dence  of  a  local  order-disorder  component  between 
polarization  vectors  in  the  unit  cell.  It  will  be  interesting  to 
measure  jil2  in  potassium  niobate  tantate,  where  again  order- 
disorder  has  been  identified,14  to  see  if  the  coefficients  are 
again  anomalously  large. 

In  conclusion,  the  PZT  ceramic  does  exhibit  modest 
flexoelectricity  and  ferroelastic  domain  wall  motion  en¬ 
hances  the  response,  but  it  is  not  possible  in  our  samples 
using  stress  levels  up  to  the  full  fracture  strength  to  induce 
ferroelectric  poling,  although  ferroelastic  poling  was  patently 
obvious. 

In  the  PZT  thin  films  epitaxially  grown  on  lattice- 
mismatched  substrates  the  strain  gradient  between  surfaces 
and  strain-relieving  dislocations  can  be  exceedingly  large 
and  the  flexoelectric  effects  could  be  of  major  importance  in 
these  systems. 
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Mechanical  strain  gradient  generated  electric  polarization  or  flexoelectric  effect  was  investigated  in 
unpoled  lead  zirconate  titanate  (PZT)  ceramics  in  the  ferroelectric  state  by  using  a  cantilevered 
beam  based  approach.  Flexoelectric  coefficient  /z12  at  room  temperature  was  measured  to  be 
1.4  fiC/m  in  the  PZT  ceramic  at  small  level  of  strain  gradient.  Temperature-dependent  experimental 
investigations  clearly  showed  that  high  dielectric  permittivity  in  the  ferroelec tries  enhanced 
flexoelectric  polarization:  essentially  a  linear  relation  was  found  to  exist  between  fixl  and  dielectric 
susceptibility  x  at  lower  permittivity  level  (2100-2800),  while  jjl12  versus  x  curve  started  to  deviate 
from  the  straight  line  at  the  ^—2800  and  nonlinear  enhancement  of  fin  with  x  was  observed,  with 
p,l2  value  reaching  9.5  at  ^ —  11  000.  The  nonlinearity  in  the  flexoelectric  effect  was  associated  with 
domain-related  processes.  It  is  suggested  that  flexoelectric  effect  can  have  a  significant  impact  on 
epitaxial  ferroelectric  thin  films  and  mesoscopic  structures.  ©  2005  American  Institute  of  Physics. 
[DOI:  10.1063/1.1868078] 


Flexoelectric  effect  is  the  coupling  between  mechanical 
strain  gradient  and  electric  polarization  and  can  be  described 
by 


Pi 


(i) 


where  Pt  is  the  flexoelectric  polarization,  the  flexoelec¬ 
tric  coefficient,  ejy-  the  elastic  strain,  and  xk  the  position  co¬ 
ordinate.  fii}ki  is  a  fourth-rank  polar  tensor  and  therefore  has 
nonzero  components  in  dielectric  solids  of  any  crystal  sym¬ 
metry.  Theoretical  estimations  predicted  that  flexoelectric  co¬ 
efficients  in  simple  dielectrics  are  generally  small  (of  the 
order  e! a  or  — 10“ 10  C/m,  where  e  is  the  electron  charge  and 
a  the  dimension  of  unit  cell).1”3  Earlier  experimental  work 
on  polymers4  supported  the  theoretical  predictions.  Recently 
we  developed  a  cantilevered  beam  based  approach5  to  per¬ 
form  reliable  measurements  of  flexoelectric  effect.  Experi¬ 
mental  investigations  using  the  approach  showed  that  flexo¬ 
electric  coefficients  in  the  ferroelectric  materials  are  many 
orders  of  magnitude  higher  (10-6  to  1CT4  C/m).6,7  The  can¬ 
tilevered  beam  based  approach  is  generally  used  for  dynamic 
flexoelectric  measurements  at  lower  level  of  strain  gradient. 
To  investigate  the  flexoelectric  effect  at  higher  strain  gradient 
we  developed  another  approach  based  on  four-point  bending 
and  used  the  approach  in  static  measurements  of  the  flexo¬ 
electric  coefficients  in  ferroelectric  materials.8 

Based  on  flexoelectric  effect,  ideas  for  developing  new 
types  of  piezoelectric  composites9  were  proposed,  where 
none  of  the  components  is  piezoelectric.  If  we  understood 
flexoelectric  effect  well  and  there  were  database  of  flexoelec¬ 
tric  coefficients  available,  a  range  of  properly  engineered 
flexoelectric  composite  structures  could  provide  completely 
new  piezoelectric  capability. 

Ferroelectric  thin  films  and  mesoscopic  structures  have 
exhibited  a  lot  of  potential  for  device  applications  such  as 
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FERAM  and  MEMS.10  Misfit  strain  certainly  exists  at  the 
interface  between  the  epitaxially  grown  thin  film  ferroelec¬ 
tric  and  the  substrate  or  electrode  material.  Nonuniform  re¬ 
laxation  of  the  misfit  strain  in  the  thin  films  can  lead  to  strain 
gradient,11  which  can  then  influence  the  dielectric  and  polar¬ 
ization  behaviors  of  the  ferroelectrics.  Experimentally,  some 
new  phenomena  recently  observed  in  the  thin  film  ferroelec¬ 
trics  were  suggested  to  be  due  to  the  consequences  of  flexo¬ 
electric  effect,  such  as  the  mechanical  stress  induced  imprint 
in  PZT  capacitor  structures12  and  pyroelectricity  in  highly 
stressed  quasiamorphous  BaTi03  films.13  Recently  a  phe¬ 
nomenological  model14  of  flexoelectricity  was  also  proposed, 
showing  that  flexoelectric  effect  could  play  an  important  role 
in  reducing  the  dielectric  maximum  in  ferroelectric  thin 
films.  Ferroelectric  PZT  based  ceramics  and  thin  films  are  a 
family  of  technologically  important  functional  materials.15  In 
this  letter,  we  report  experimental  investigations  of  the  flexo¬ 
electric  effect  in  ceramic  PZT. 

The  PZT  samples  are  the  same  as  those  used  before,8 
with  dimensions  60  mm  long,  7  mm  wide,  and  3  mm  thick. 
Flexoelectric  effect  was  investigated  using  the  cantilevered 
beam  based  approach  and  the  system  used  for  the  flexoelec¬ 
tric  measurement  is  identical  to  that  used  in  earlier  studies  of 
barium  strontium  titanate  (BST)  and  lead  magnesium  niobate 
(PMN).5*7  A  series  of  3-mm-diam  thin  sputtered  gold  elec¬ 
trodes  were  prepared  on  the  sample  surface  along  the  bar 
length.  The  ceramic  bar  sample  is  rigidly  clamped  at  one  end 
and  driven  into  transverse  vibration  at  1  Hz  by  a  small  mov¬ 
ing  coil  loudspeaker.  The  ac  mechanical  strain  as  a  function 
of  position  along  the  bar  was  measured  using  a  Micros¬ 
train™  DVRT  (differential  variable  reluctance  transducer) 
and  the  generated  current  was  measured  using  an  SR830 
DSP  Lock-in  amplifier.  The  measured  mode  shape  was  then 
used  to  calculate  the  strain  gradient  at  positions  of  the  elec¬ 
trodes  using  the  free  bar  model.5 

For  an  unpoled  PZT  ferroelectric  ceramic  at  the  morpho- 
tropic  phase  boundary,  although  individual  grains  may  have 
lower  symmetry  (tetragonal  or  rombohedral)  which  permits 
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Strain  gradient  (rn1) 

FIG.  1.  Room  temperature  flexoelectric  polarization  vs  strain  gradient  for 
unpoled  PZT  ceramic  measured  at  1  Hz  near  the  clamped  end  (xxlL 
=0.18)  of  the  cantilevered  beam. 


20  40  60  80  100  120  140  160  180  200 

Temperature  fC) 

FIG.  2.  Flexoelectric  polarization  (arbitrary  unit)  measured  as  a  function  of 
temperature  in  the  ferroelectric  state  for  the  unpoled  PZT  ceramic. 


piezoelectricity,  in  the  volume  of  the  ceramic  one  may  expect  _  £ 

a  macroscopic  symmetry  of  <» oom>  so  the  nonzero  compo-  ^ 12  ~~  ^ a 


(3) 


nents  for  the  flexoelectric  tensor  should  be  /*ni2,  1122, 

and  /Ai2i2>  or  in  matrix  notation  jilh  /x12*  and  /I44.  In  the 
unpoled  PZT  ceramic,  no  residual  piezoelectricity  could  be 
detected  by  Berlincourt  <i33  meter  and  no  evidence  of  piezo¬ 
electric  resonance  could  be  found  in  the  impedance  trajectory 
from  1  kHz  to  1  MHz.  By  using  the  cantilevered  beam  based 
approach,  any  remnant  piezoelectricity  from  the  top  and  bot¬ 
tom  halves  of  the  sample  bar  is  well  balanced  during  the 
flexoelectric  measurements.  That  some  unbalance  was  affect¬ 
ing  the  measured  value  was  however  ruled  out  by  the  simple 
experiment  of  inverting  the  sample  and  noting  that  the  mea¬ 
sured  signal  did  not  change  either  in  amplitude  or  phase. 
Clearly  in  ferroelectrics  the  free  surface  breaks  the  symmetry 
of  the  bulk  and  may  affect  the  polarization  behavior,  how¬ 
ever,  in  the  present  studies  the  effect  of  surface 
ferroelectricity16  is  unlikely  due  to  the  highly  conductive 
metal  electrodes  on  both  free  surfaces.  Thus  we  believe  that 
the  measured  electric  polarization  P3  is  solely  due  to  the 
strain  gradient  in  the  jc3  thickness  direction  and  can  be  writ¬ 
ten  as 


where  7  is  a  scaling  factor.  In  order  to  evaluate  the  impact  of 
dielectric  property  on  the  flexoelectric  effect,  dielectric  spec¬ 
tra  from  the  PZT  ceramic  were  measured  and  shown  in  Fig. 
3,  where  a  strong  but  rounded  dielectric  maximum  suggests  a 
diffuse  phase  transition  around  220  °C.  The  plot  of  /x12  ver¬ 
sus  x  is  shown  in  Fig.  4.  Comparison  of  experimental  data 
obtained  in  BST,  PZT,  and  PMN  reveals  that,  at  similar  level 
of  relative  dielectric  permittivity,  fil2  in  BST  is  roughly  one 
order  of  magnitude  higher  than  that  in  PMN  or  PZT.  It  re¬ 
mains  unclear  as  to  why  the  Pb-based  ferroelectrics  should 
have  lower  values  of  flexoelectric  coefficients. 

As  shown  in  Figs.  2  and  4,  the  initial  drop  of  iiX2  with  x 
is  unexpected  but  may  correspond  to  the  increase  of  loss 
tangent  (Fig.  3)  in  the  temperature  range  (24-50  °C).  A  lin¬ 
ear  relation  between  /z12  and  x  with  7=  1  is  seen  to  exist  in 
the  x  range  of  2100-2800  (or  temperature  range  of  40- 
74  °C).  The  fjbn  versus  x  curve  deviates  from  the  7=  1  line  at 
the  ^~2800  and  the  rise  of  fil2  with  x  becomes  essentially 
nonlinear,  with  iil2  reaching  9.5  at  ^  —  11  000. 


Nonlinearity  in  flexoelectric  coefficient  was  found  in 
<9e13  temperature-dependent  experimental  investigations  of  BST,7 

=  (2)  where  /il2  was  raised  to  100  /z C/m  when  approaching  the 

1  dielectric  peak  from  the  paraelectric  state.  Such  nonlinear 

Figure  1  presents  the  room  temperature  (24  °C)  flexo¬ 


electric  polarization  as  a  function  of  the  transverse  strain 
gradient  obtained  near  the  clamped  end  (xi/L=0.18)  of  the 
bar  by  using  the  cantilevered  beam  based  dynamic  approach. 
It  is  clear  that  the  generated  electric  polarization  is  linearly 
proportional  to  the  elastic  strain  gradient  and  the  slope  of  the 
line  gives  a  magnitude  of  /x12  of  1.4  pC/m,  close  to  the  early 
data  of  0.5  /xC/m8  obtained  by  static  measurements  using  a 
four-point  bend  fixture.  Figure  2  presents  the  measurements 
of  nn  in  the  PZT  ceramic  as  a  function  of  temperature  (24- 
180  °C).  When  heating  up  from  room  temperature  (24  °C), 
unexpectedly  the  /*12  initially  drops  before  it  becomes  stabi¬ 
lized  at  40  °C,  then  flexoelectric  polarization  basically  keeps 
flat  at  the  temperature  range  of  40-70  °C  but  starts  to  rise 
prominently  at  70  °C. 

Clearly  high  dielectric  permittivity  in  the  ferroelectric 
materials  can  enhance  flexoelectric  coefficients6,7  and  /z12 


Temperature  (°C) 

FIG.  3.  Dielectric  permittivity  and  loss  tangent  for  the  unpoled  PZT  ceramic 


can  be  related  to  dielectric  susceptibility  V  by  measured  as  function  of  temperature  at  a  frequency  of  1  kHz. 
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FIG.  4.  Flexoelectric  coefficient  /z12  vs  dielectric  susceptibility  x  f°r  the 
unpoled  PZT  ceramic.  The  dotted  line  shows  a  linear  relation  between  /tl2 
and  x  at  lower  level  of  dielectric  permittivity  with  a  y=l. 

enhancement  was  suggested  to  be  due  to  the  survival  of 
some  ferroelectric  domains  in  the  BST  ceramic  above  the 
Curie  transition  temperature  Tc.  In  the  previous  studies  of 
ceramic  PMN,6  nonlinear  enhancement  of  fxn  with  x  was 
also  observed  and  found  to  be  closely  associated  with  the 
pre-existing  polar  microdomains  in  this  prototypic  relaxor 
ferroelectric  material.17  Apart  from  the  nonlinear  behavior  of 
flexoelectric  coefficient  with  dielectric  permittivity,  previous 
static  investigations8  in  PZT  ceramic  revealed  nonlinear  re¬ 
lation  between  flexoelectric  coefficient  and  mechanical  strain 
gradient,  where  the  jump  of  /z12  from  0.5  /zC/m  at  lower 
strain  gradients  to  2  /xC/m  at  much  higher  strain  gradients 
was  shown  to  be  associated  with  the  onset  of  domain  wall 
motion  induced  by  the  very  high  level  of  inhomogeneous 
strain  achievable  with  the  four-point  bend  fixture.  As  dis¬ 
cussed  earlier,  nonlinearity  in  the  flexoelectric  effect  can  be 
tentatively  attributed  to  the  domain-related  processes  in  fer- 
roelectrics  and  the  scaling  factor  y  in  Eq.  (3)  depends  on  the 
materials  investigated,  the  level  of  dielectric  permittivity,  and 
the  magnitude  of  strain  gradient.  The  experimental  findings 
of  nonlinear  phenomena  are  in  good  agreement  with  the  the¬ 
oretical  studies  by  Catalan  et  a/.,14  where  it  was  shown  that 
flexoelectric  coefficient  is  only  a  linear  function  of  the  strain 
gradient  or  permittivity  when  the  induced  flexoelectric  polar¬ 
ization  is  small.  In  addition  to  the  extrinsic  influence  of 
domain-related  processes  discussed  earlier,  it  is  possible  that 
a  part  of  the  nonlinearity  in  the  flexoelectric  effect  may  also 
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be  of  intrinsic  origin  as  shown  in  the  theoretical  work  of 
Catalan  et  al 14 

In  epitaxial  ferroelectric  thin  films  and  mesoscopic  struc¬ 
tures,  Curie  phase  transition  temperature  Tc  or  dielectric 
peak  can  be  adjusted  to  around  room  temperature  where  ul- 
trahigh  dielectric  permittivity  becomes  available,  by  appro¬ 
priate  selection  of  lattice  misfit  and  film  thickness  as  shown 
in  recent  theoretical18  and  experimental19  studies.  Likewise, 
strain  gradient  in  these  ferroelectric  structures  can  also  be 
adjusted  by  tailoring  the  structure  dimensions  or  controlling 
the  relaxation  of  the  misfit  strain.  The  nonlinear  enhance¬ 
ment  of  flexoelectric  coefficients  with  the  dielectric  permit¬ 
tivity  and  strain  gradient  can  lead  to  a  significant  impact  of 
flexoelectric  effect  in  these  properly  engineered  thin  film 
ferroelectric  heterostructures. 

In  summary,  flexoelectric  effect  was  investigated  in  the 
ferroelectric  state  of  unpoled  lead  zirconate  titanate  (PZT) 
ceramics.  Temperature-dependent  flexoelectric  investigations 
showed  that  /z12  essentially  increases  with  relative  dielectric 
permittivity  and  the  nonlinear  phenomenon  was  found  at 
higher  level  of  dielectric  permittivity,  which  is  suggested  to 
be  associated  with  domain-related  processes. 
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Intermediate  states  in  Pb(B{/3B2/3)03-PbTi03 
ferroelectric  single  crystals 

A.  Amin  and  L.  E.  Cross 


Optical  and  dielectric  studies  of  oriented  and  poled 
single  crystal  compositions  on  both  sides  of  the  morpho- 
tropic  phase  boundary  in  (1  —xlPb/ZnujNbipJOT.'- 
xPbTiOT,  have  revealed  an  intermediate,  metastable, 
ferroelectric  orthorhombic  state  (F0)  in  a  narrow 
composition  range  between  the  ferroelectric  rhombo- 
hedral  (FR)  and  tetragonal  (FT)  phases,  which  has 
also  been  identified  from  pyroelectric  and  dielectric 
data.  Recently  published  structural  data  confirm  this 
orthorhombic  state  for  008  <  x  <  0-11.  Similarly,  dielec¬ 
tric  studies  and  electromechanical  measurements  on 
rhombohedral  Pb(Mglf3Nb2/3)03-PbTi03  single  crys¬ 
tals  have  revealed  an  intermediate,  metastable  F0  state 
at  high  temperature.  Structural  evidence  of  the  ortho¬ 
rhombic  phase  in  this  system  has  also  recently  been 
reported.  In  the  present  paper  it  is  shown  that  a 
metastable  orthorhombic  state  is  always  close  to  but 
unstable  with  respect  to  the  FR-FT  degenerate  states 
(morphotropy),  as  predicted  by  a  phenomenological 
model  for  A( B[ -XBX  )03  ferroelectrics.  Conditions  for 
stabilising  the  nearly  degenerate  FR,  Fr  and  F0  states 
in  the  Pb(BIll3BI2I/3)03~PbTi03  single  crystal  system 
are  discussed  in  terms  of  the  energy  function  tensor 
coefficients.  BCT/656 
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INTRODUCTION 

The  large  longitudinal  piezoelectric  coefficient  d33  ( ^  1 500 
pC  N_1)  and  electromechanical  coupling  factor  k33  (^0*92) 
of  <001)  oriented  and  poled  multidomain  single  crystals  of 
0-91Pb(Zn1/3Nb2/3)03-0-09PbTi03  were  first  reported  by 
Kuwata  et  alf2  who  examined  relaxor  ferroelectric  single 
crystals  by  a  standard  resonance-antiresonance  method 
under  a  low  drive  field  level,  significantly  below  that 
required  to  induce  saturation  in  electrically  induced  strain. 
In  the  relaxor-ferroelectric  system  (1  -.\')Pb(Bi/3B2/3)03- 
.\'PbTi03,  where  B1  =  Mg,  Zn  or  Sc  and  Bu  =  Nb,  lead  zinc 
niobate-lead  titanate  (PZN-PT)  and  lead  magnesium 
niobate-lead  titanate  (PMN-PT)  have  attracted  much 
attention  in  recent  years.3,4  Research  has  focused  on  com¬ 


positions  close  to  the  morphotropic  phase  boundary  (MPB) 
between  rhombohedral  and  tetragonal  ferroelectric  states. 
Unprecedented  strain  levels  in  excess  of  1*5%  have  been 
induced  in  rhombohedral  morphotropic  PZN-PT  92/8  and 
PMN-PT  70/30  single  crystals  with  an  electric  field 
£^3  MV  m”1  applied  along  the  [001]  axis.3,4  These  strain 
levels  are  about  an  order  of  magnitude  larger  than  in 
conventional  lead  zirconate  titanate  (PZT)  piezoceramics. 

Because  of  their  high  electromechanical  coupling,  single 
crystal  compositions  offer  the  potential  for  both  broadband 
and  high  acoustic  energy  density  when  compared  with 
standard  PZT8.  They  show  great  promise  in  medical  ultra¬ 
sound  imaging  systems,  underwater  sonar  projectors  and 
receiving  arrays.  There  is  the  potential  for  a  vast  improve¬ 
ment  in  the  axial  resolution  and  contrast  of  medical  ultra¬ 
sound  and  sonar  systems  as  a  result  of  the  broadband 
capabilities  of  single  crystals.  The  origin  of  this  enhanced 
electromechanical  behaviour  has  been  attributed  to  an 
electrically  induced  rhombohedral  ferroelectric  FR  (space 
group  R3m)  to  tetragonal  ferroelectric  FT  (space  group 
P4mm)  phase  transformation.3"7  The  FR  phase  has  a  much 
smaller  c/a  ratio  than  FT,  and  accordingly  under  field  the 
large  c/a  ratio  change  results  in  enhanced  electrically 
induced  strain. 

INTERMEDIATE  STATES 

Optical  studies  augmented  by  high  temperature  dielectric 
measurements  on  [001],  [110]  and  [111]  oriented  and 
poled  compositions,  on  both  sides  of  the  MPB  in 
(1  —  x)Pb(Zn1/3Nb2/3)03-xPbTi03  system,  have  revealed 
the  existence  of  an  intermediate,  metastable,  orthorhombic 
state  for  the  poled  multidomain  [110]  and  poled  mono- 
domain  [111]  orientations,  but  not  for  poled  multidomain 
[001]. 8  Synchrotron  X-ray  powder  diffraction  measure¬ 
ments  confirm  that  an  intermediate  ferroelectric  ortho¬ 
rhombic  state  F0  (space  group  Bmm2)  exists  for  a  narrow 
composition  range  (0-08  <x<0’ll)  with  near  vertical  phase 
boundaries  on  both  sides  of  the  MPB.9  For  x  =  0  08  and 
0*09  an  electrically  induced  orthorhombic  state  has  been 
reported.10"12 

Dielectric  and  optical  studies  of  [110]  oriented  and 
poled  ferroelectric  rhombohedral  (FR)  multidomain 
0'67Pb(Mg1/3Nb2/3)03-0*33PbTi03  crystals  revealed  a 
transformation  sequence  on  heating  to  an  intermediate 
ferroelectric  orthorhombic  state  (F0)  near  80°C,  followed 
by  F0-Ft  transition  at  100°C.  Transition  to  the  cubic  state 
occurred  at  145°C.13  On  the  other  hand  at  a  higher  poling 
field  (>5  kV  cm-1)  an  orthorhombic  monodomain  F0  state 
was  induced  at  room  temperature.  The  FQ-FT  transition  of 
the  monodomain  state  occurred  at  8(PC,  much  lower  than 
for  the  multidomain  state.  The  varied  temperature  range  of 
existence  of  the  intermediate  F0  state  suggested  metast¬ 
ability.13  No  evidence  of  an  intermediate  F0  state  was  found 
for  the  poled  multidomain  [100]  crystal. 

Further  evidence  of  an  intermediate  high  temperature 
F0  state  between  the  FR  and  FT  phases  has  recently  been 
reported  for  0-70Pb(Mg1/3Nb2/3)03-0-30PbTi03  single 
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1  Variation  of  piezoelectric  coefficient  d3 3  with  temperature 
for  poled  multidomain  [110]  oriented 
0*70Pb(Mg,/3Nb2/3)03-0*30PbTi03  single  crystal: 
arrows  indicate  ferroelectric  rhombohedral  FR- 
ferroelectric  orthorhombic  Fo  transition  (^85°C)  and 
ferroelectric  orthorhombic  F0-ferroelectric  tetragonal 
Ft  transition  («100°C);  Curie  temperature  of  this 
composition  is  « 125°C.  Inset  shows  variation  of  Young’s 
modulus  with  temperature  near  two  ferroelectric 
transitions  (data  from  Ref.  14) 


crystals.14  Figure  1  plots  longitudinal  piezoelectric  response 
versus  temperature  for  the  [110]  oriented  and  poled  multi- 
domain  crystal.  The  inset  illustrates  the  temperature  depen¬ 
dence  of  Young’s  modulus.  The  data  in  Fig.  1  are  taken 
from  Ref.  14,  excluding  the  responses  from  the  other  two 
orientations  for  clarity.  Note  the  excellent  piezoelectric 
properties  of  the  [110]  oriented  crystal.  Arrows  indicate 
the  Fr-F0  transition  near  85°C  and  the  F0-Fx  transition 
near  100°C,  with  the  Curie  transition  occurring  at  125°C14 
The  dielectric  data  are  qualitatively  similar  to  previously 
published  results.13  The  data  suggest  that  the  boundary 
between  the  FR  and  FM  phases  near  85°C  has  a  thin 
orthorhombic  sliver  extending  over  15°C.  This  supports  the 
recent  suggestion  of  mixed  states  in  this  system.15  No 
evidence  of  an  intermediate  F0  state  was  found  for  the 
poled  multidomain  [100]  crystal.  Simulations  from  first 
principles  indicate  that  the  transformation  under  electric 
field  between  FR  and  Fx  proceeds  by  rotation  of  the 
polarisation  between  <111)  and  <00 1>  via  an  intermediate 
monoclinic  or  orthorhombic  ferroelectric  state.7  These 
results  open  up  the  interesting  possibility  that  intermediate 
states  may  be  responsible  for  the  enhanced  electromechani¬ 
cal  performance. 

Recent  synchrotron  X-ray  powder  diffraction  experi¬ 
ments  have  revealed  a  new  phase  diagram  for  the 
(1  —  x)Pb(Mg1/3Nb2/3)03-xPbTi03  system.15  The  proposed 
phase  diagram  is  a  complex  landscape  of  phase  mixtures. 
Morphotropic  compositions  O31^x<0-33  have  a  mono¬ 
clinic  symmetry  Mc  (space  group  Pm)  with  the  unique 
axis  oriented  along  the  pseudocubic  [010]  direction.  For 
x=0*31  the  symmetry  is  Very  close  to  orthorhombic,  and  it 
is  possible  for  compositions  with  x  between  0-30  and  0-31 
to  be  entirely  orthorhombic.  In  this  case  the  transition 
between  rhombohedral  and  monoclinic  Mc  phases  will 
occur  via  the  orthorhombic  phase.  For  x<0-31  the  sym¬ 
metry  is  rhombohedral,  and  for  x^O-37  it  is  tetragonal. 

Powder  diffraction  patterns  for  the  relax  or-ferroelectric 
system  ( 1  —  x)Pb(BI1/3B,2I/3)03-xPbTi03  are  extremely  diffi¬ 
cult  to  analyse.915  Two  major  factors  contribute  to  this 
complexity,  and  to  the  difficulty  of  unequivocal  symmetry 
determination,  especially  for  compositions  close  to  the 
MPB.  The  first  is  inhomogeneity  due  to  the  concentration 


gradient,  which  is  typical  of  crystal  growth  from  the  melt 
of  incongruently  melting  compositions,  and  the  second  is 
the  broad  diffraction  maxima  characteristic  of  (disordered) 
relaxor-ferroelectric  systems.  This  complicates  peak  decon¬ 
volution  of  the  overlapping  or  partially  overlapping  reflec¬ 
tions  common  in  pseudosymmetric  structures,  for  example 
the  Fr,  Fx  and  F0  phases.  To  obtain  sharp  diffraction 
maxima,  an  electric  field  is  applied  to  the  crystal  to  induce 
an  ordered  state,  and  the  powder  pattern  is  then  recorded 
for  the  ordered  state.9  Conversely,  the  possibility  of  a  field 
induced  phase  transition  and  of  assigning  a  lower  symmetry 
for  mixed  phase  patterns  also  exists. 


ENERGY  FUNCTION 

To  account  for  the  recently  observed  monoclinic  phase  in 
PZT,  an  extension  of  the  Devonshire  theory  of  perovskite 
ferroelectrics  has  been  developed.16  In  this  model,  the 
ground  state  energy  is  represented  by  an  analytic  function 
defined  on  the  surface  of  a  sphere  and  truncated  at  the 
eighth  power  of  the  order  parameter.  The  order  parameter 
in  this  case  is  the  polarisation  vector  normalised  to  unity. 
The  parameter  space  is  thus  a  unit  sphere  spanned  by  the 
polar  and  azimuth  angles  a  and  /?  respectively,  where 
0  <  a  <  7T  and  0  <  ft  <  2n.  The  ferroelectric  phases  are 
obtained  by  numerically  searching  for  constrained  minima 
of  energy  as  a  function  of  orientation.  The  constraints  are 
taken  such  that  the  order  parameter  must  lie  along  a 
symmetry  axis  [001],  [111]  or  [Oil]  for  ferroelectric 
tetragonal,  rhombohedral  and  orthorhombic  phases 
respectively,  or  in  a  symmetry  plane  (mirror)  for  the  mono- 
clinic  ferroelectric  phases.  Three  monoclinic  solutions  are 
admissible.  These  are  denoted  by  MA  (Cm),  MB  (Bm)  and 
Mc  (Pm),  the  space  group  symbol  being  given  between 
parentheses.  The  polarisation  vector  is  constrained  in  a 
mirror  plane  along  [uuv]  with  u<v,  [uuv]  with  u>v,  or 
[Ouv]  for  each  of  the  monoclinic  species  respectively.  A  two 
dimensional  representation  of  the  ferroelectric  phase  dia¬ 
gram  in  terms  of  the  linearly  independent  parameters  a  and 
ft  is  obtained  by  mapping  the  unit  sphere  onto  a  plane.  All 
six  ferroelectric  phases  enumerated  above  are  accessible  by 
the  eighth  power  model.  This  model  can  only  describe 
transitions  among  the  ferroelectric  phases,  and  not  trans¬ 
itions  to  the  high  temperature,  higher  symmetry  phases.  It 
cannot  access  the  rhombohedral  (R3m)-rhombohedral 
(R3c)  transition,  associated  with  condensation  of  the  triply 
degenerate  Brillouin  zone  corner  R25  lattice  mode  observed 
in  the  PZT  system.  Besides,  temperature,  compositional, 
elastic,  and  electric  variables  are  not  explicit.  To  access  a 
triclinic  ferroelectric  phase,  a  twelfth  order  model  is 
required. 

The  observed  metastable  orthorhombic  state  in  the 
single  crystal  Pb(Zn1/3Nb2/3)03-PbTi03  system  and  the 
high  temperature  orthorhombic  state  in  single  crystal 
Pb(Mg1/3Nb2/3)03-PbTi03  can  be  readily  visualised,  and 
further  insights  can  be  gained  into  these  complex  systems 
within  the  limitations  of  a  sixth  order  theory.  The  energy 
landscape  in  terms  of  composition  and  temperature  has 
been  previously  developed  by  the  present  authors  and 
co workers 17-20  for  a  generalised  A(B,1„xB")03  MPB  phase 
diagram  between  ferroelectric-ferroelectric  and  ferro- 
electric-antiferroelectric  end  members.  For  the  PZT  system 
the  theory  predicts  the  position  of  the  MPB  as  a  function 
of  composition  and  temperature  and  the  rhombohedral 
(R3m)-rhombohedral  (R3c)  phase  boundary  associated 
with  the  triply  degenerate  Brillouin  zone  corner  R25  lattice 
mode.  Semiquantitatively,  the  single  crystal-single  domain 
dielectric,  piezoelectric  and  elastic  properties  are  accessible 
in  terms  of  temperature,  composition,  stress19  and  electric 
field.20  The  derived  single  crystal  polarisation  versus  temper¬ 
ature  results  for  the  end  member  PbTi03  are  for  example 
within  a  few  per  cent  of  the  published  experimental  results.17 
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Composition 


2  Generalised  A(B5-xB")03  phase  diagram  as  computed 
from  free  energy  minima  of  equation  (1),  showing 
equilibrium  states  and  near  vertical  MPB  between  FR 
and  Fx  states17 


It  must  again  be  emphasised  that  the  monoclinic  and 
triclinic  ferroelectric  solutions  are  outside  the  region  of 
validity  of  this  sixth  power  function. 

The  free  energy  function  Gt  for  a  simple  proper  ferroelec¬ 
tric  derived  from  a  prototypic  phase  with  symmetry  Pm3m, 
considering  only  Brillouin  zone  centre  lattice  modes,  takes 
the  following  form  for  electric  and  elastic  free  boundaries: 

=  a  t{P\  +  P\  +  P\)  +  an(Pi  +  Pi  +  Pt) 

+  aX2{P\P\  +  PlPl  +  P\P\)  +  a  lu(Pf  +  Pf  +  P63) 

+  a  iMPl  +  Pi)  +  Pi(Pl  +  P\)  +  Pt(Pl  +  Pi) 

+  a  l23PiP22Pi . (1) 

In  this  equation  the  function  is  truncated  at  the  sixth  order 
invariant  of  the  order  parameter,  i.e.  the  dielectric  polaris¬ 
ation  vector  P.  The  surviving  aijk  coefficients  are  determined 
by  the  generating  elements  of  the  point  group  m3m.  The 
coefficient  is  a  linearly  increasing  function  of  temperature 
passing  through  0  at  the  Curie-Weiss  temperature  6 .  The 
higher  order  aijk  coefficients  are  related  to  the  dielectric 
stiffness  and  higher  order  stiffness  coefficients  that  are 
required  to  account  for  dielectric  nonlinearity  of  the  ferro¬ 
electric  states. 

Careful  numerical  analysis  suggests  that  even  if  the  aijk 
coefficients  in  equation  (1)  are  given  a  linear  temperature 
dependence,  a  rhombohedral  (R3m)-tetragonal  (P4mm) 
degeneracy  (morphotropy) .  cannot  be  developed  from  a 
simple  Devonshire  function  (am  =  —  By  adding  a 

positive  aU2  term  to  destabilise  the  orthorhombic  phase 
and  a  negative  a123  to  restabilise  the  rhombohedral  phase,  a 
morphotropic  phase  boundary  (Fig.  2)  between  the 
rhombohedral  and  tetragonal  states  can  be  generated  over 
a  wide  temperature  range  without  any  interleaving  ortho¬ 
rhombic  region.  Numerical  values  of  the  energy  density 
function  Gx  are  illustrated  in  Fig.  3  for  PZT.  It  is  important 
to  reemphasise  that  the  orthorhombic  ferroelectric  phase 
(filled  circles  in  Fig.  3)  is  always  close  to  but  metastable 
with  respect  to  the  rhombohedral  and  tetragonal  phases. 


MOLE  FRACTION  PbTIO,  IN  PbZrOiPbTiO* 


3  Numerical  values  of  energy  density  function  (7,  versus 
composition  and  temperature  for  FR,  FT,  and  F0  states 
of  Pb(Zr,_xTix)03  (PZT)  system:  note  near  vertical 
MPB  (dashed  line)  near  Zr:Ti  ratio  of  unity  and  nearly 
degenerate  FR,  FT  and  F0  states  over  narrow  composition 
range  (jc%0*3)  at  around  300°C  (Ref.  17) 

In  summary,  recent  dielectric,  piezoelectric,  elastic,  optical 
and  structural  investigations  have  demonstrated  a  complex 
transformation  sequence  in  A(BI1_xBjcI)03‘  ferroelectric 
piezocrystals  that  is  dependent  upon  orientation,  domain 
states  and  electrical  history.  An  intermediate  Fc  state  was 
shown  to  exist  between  FR  and  FT,  the  temperature  range 
of  existence  of  which  varied,  suggesting  metastability. 
Furthermore,  the  results  indicate  that  the  margin  of  stability 
in  the  FR,  F0  and  FT  states  is  sufficiently  small  that 
changes  in  orientation,  electrical  history  or  mechanical 
constraints  may  affect  the  free  energy  balance.  The  data 
suggest  that  the  boundary  between  the  FR  and  FT  phases 
near  85°C  has  a  thin  orthorhombic  sliver  extending  over 
15°C.  This  result  supports  the  recent  suggestion  of  mixed 
states  in  this  system.15 

Phenomenological  models  that  explicitly  include  temper¬ 
ature,  compositional,  elastic  and  electrical  boundaries  are 
valuable  for  computing  phase  stability  and  physical  prop¬ 
erties  as  a  function  of  these  boundaries.  Deeper  insights 
can  be  gained  from  model  extension  to  account  for  the 
monoclinic  and  triclinic  solutions  where  needed,  and  to 
compute  the  single  crystal,  single  domain  macroscopic 
observables  and  the  influence  of  elastic  and  electrical 
boundaries. 
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The  dielectric  behavior  of  95.5%  Pb(Zn1/3Nb2/3)03-4.5%  PbTi03  single  crystals  oriented  along 
(001)  direction  with  and  without  dc  electric  field  has  been  studied  at  cryogenic  temperatures.  A 
pronounced  low-temperature  dielectric  relaxation  process  was  observed  below  200  K;  the  relaxation 
rate  follows  the  Arrhenius  law  (t0  =  ~-  1.0X  10“15  s  and  C/=0.24eV).  An  additional  dielectric 
anomaly  showed  up  around  250  K  at  10  kHz  under  a  dc  electric  field.  These  results  indicate  rather 
complicated  polarization  mechanisms  at  ciyogenic  temperatures  which  clearly  need  more  detailed 
study.  The  strain  levels  at  cryogenic  temperatures  suggest  that  this  material  is  very  promising  for 
space  applications,  in  which  the  performance  at  cryogenic  temperatures  is  critical.  ©  2003 
American  Institute  of  Physics.  [DOI:  10.1063/1.1541936] 


Since  the  discovery  of  the  electric-field-induced  ultra- 
high  strain  (>1%)  in  lead-based  single  crystals 
Pb(Zn1/3Nb2/3)03  and  Pb(Zn1/3Nb2/3)03-PbTi03(PZN- 
PT), 1,2  much  work  has  been  focused  on  the  growth  and  char¬ 
acterization  of  related  lead-based  single  crystals  as  well  as  on 
the  understanding  of  the  mechanisms  governing  their  ferro¬ 
electric  and  strain  behavior.  Moreover,  many  efforts  have 
been  devoted  to  establishing  complete  phase  diagrams  over  a 
wide  temperature  range  in  order  to  better  understand  the  re¬ 
lationships  between  the  phase  transformations  and  the  strain 
characteristics  under  electric  fields,3-6  in  which  x-ray  mea¬ 
surements  with  or  without  fields  were  performed.  It  is  well 
known  that  the  study  of  dielectric  spectra  is  a  powerful  tool 
in  determining  various  dielectric  peaks  as  well  as  dielectric 
relaxation  behavior,  which  is  critical  in  understanding  the 
dielectric  polarization  mechanisms  in  relation  to  piezoelec¬ 
tric  strains.  In  addition,  the  studies  of  dielectric  properties 
with  or  without  an  electric  field  will  also  provide  basic  ther¬ 
modynamic  parameters  for  phenomenological  calculations. 
Most  work  on  dielectric  properties  was  carried  out  above 
room  temperature  or  a  little  below.  Dielectric  spectra  at  even 
lower  temperatures,  especially  below  liquid  nitrogen  tem¬ 
perature,  are  rarely  studied,  in  particular  under  a  dc  electric 
field.7-10 

The  discovery  of  very  high  strain  in  PZN-PT,  single 
crystals  also  motivates  possible  application  of  the  material  in 
space  technology,  in  which  it  is  demanded  that  actuators 
have  a  reasonably  good  performance  at  cryogenic  tempera¬ 
tures,  for  example,  around  100  K  and  even  lower.  A  question 
is  raised,  “How  much  of  the  excellent  strain  at  room  tem¬ 
perature  remains  at  such  low  temperatures?” 

In  this  work,  we  study  the  dielectric  properties  from  300 
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K  down  to  12  K,  with  and  without  a  bias  electric  field.  The 
results  show  various  dielectric  anomalies  with  frequency  dis¬ 
persion  in  the  PZN-PT  system  with  or  without  dc  electric 
fields  below  room  temperature.  The  strain  levels  are  also 
measured  down  to  110  K;  preliminary  results  show  that  the 
PZN-PT  single  crystals  have  great  potential  for  applications 
in  space. 

Two  (OOl)-oriented  single-crystal  samples  with  the  com¬ 
position  of  95.5%  Pb(Znj/3Nb2/3)03-4.5%  PbTi03  (hereaf¬ 
ter  denoted  as  PZN-4.5  PT)  were  used  for  this  study.  The 
crystal  orientation  was  determined  by  the  Laue  back- 
reflection  technique  using  a  Northstar  real-time  orientation 
system.  The  dielectric  constant  and  loss  of  the  samples  were 
measured  using  an  HP4284A  LCR  meter  in  the  temperature 
range  12  to  300  K  with  a  cooling/heating  rate  of  1  K/min.  A 
dc  voltage  was  applied  to  the  sample  and  a  blocking  circuit 
was  adopted  to  isolate  the  high  dc  voltage  from  the  LCR 
meter. 

The  temperature  dependence  of  the  dielectric  constant 
(e)  and  loss  (tan  5)  as  a  function  of  frequency  of  PZN- 
4.5  PT  single  crystals  measured  in  a  cooling  cycle  is  shown 
in  Fig.  1.  From  300  to  12  K,  an  interesting  characteristic  is 
that  there  are  two  dielectric  relaxation  processes,  as  seen 
clearly  in  tan  8  One  is  a  set  of  broad  tan  S  peaks  (denoted  as 
peak  A)  around  — 100  K,  which  are  strongly  frequency  de¬ 
pendent.  The  other  is  around  250  K  (denoted  as  peak  B). 

For  peak  A,  the  relaxation  rate  for  the  polarization  de¬ 
rived  from  the  temperature  dependence  of  the  imaginary  part 
of  the  permittivity  for  the  sample  is  plotted  in  the  inset  of 
Fig.  1.  The  data  were  fitted  to  the  Arrhenius  relation 

r=  r0exp[  Ul  kBT]9  (1) 

where  ris  the  relaxation  time,  r0  is  the  pre-exponential  term, 
U  the  activation  energy  for  relaxation,  £B  the  Boltzmann’s 
constant,  and  T  the  temperature.  The  fit  parameters  are  U 
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FIG.  1.  Temperature  dependence  of  e  and  tan  S  of  PZN-4.5  PT  single  crys¬ 
tals  measured  in  a  cooling  cycle  at  10,  100,  and  1000  kHz  (e:  from  top  to 
bottom;  tan  &  from  bottom  to  top).  Inset  shows  the  relaxation  rate  {p)  as  a 
function  of  inverse  temperature  (circles:  experimental  data;  solid  line:  fitting 
to  the  Arrhenius  relation). 


FIG.  3.  Unipolar  strain  ( S )  vs  electric  field  ( E )  of  PZT-4.5  PT  single 
crystals  measured  at  0.2  Hz  at  cryogenic  temperatures:  296,  248,  229,  190, 
and  148  K,  labeled  as  1  to  5,  respectively.  Inset  shows  the  temperature 
dependence  of  the  strain  level  at  60  kV/cm. 


=0.24  eV,  and  r0  =  -  1.0X  10" 15  s.  The  fitting  curve  is 
shown  by  the  solid  line  in  the  inset  of  Fig.  1.  A  similar 
dielectric  relaxation  behavior  was  also  observed  in  other 
samples  with  the  activation  energy  of  0.22±0.2  eV. 

It  is  known  that  the  thermal/electric  field  history  has 
great  influence  on  the  dielectric  and  piezoelectric/ 
ferroelectric  behavior  of  PZT-4.5  PT  single  crystals  above 
room  temperature.  It  is  interesting  to  study  the  electric-field 
effect  on  the  dielectric  and  ferroelectric  properties  of  the  sys¬ 
tem  at  cryogenic  temperatures.  In  this  work,  we  measured 
the  temperature  dependence  of  the  dielectric  constant  and 
loss  with  and  without  dc  fields  (i.e.,  0-10  kV/cm).  The  tem¬ 
perature  dependence  of  e  and  tan  8  at  10  kHz  under  bias 
fields  of  0,  5,  and  10  kV/cm  measured  in  a  cooling  cycle  is 
shown  in  Fig.  2.  An  obvious  dielectric  peak,  peak  B,  in  tan  8 
shows  up  under  a  dc  electric  field  of  10  kV/cm  around  250 
K.  Although  no  obvious  dielectric  peak  in  e  is  present,  a 
change  in  the  slope  of  s  occurs  at  this  particular  temperature. 


In  fact,  even  without  dc  bias,  a  slim  trace  of  peak  B  can  be 
also  seen,  as  shown  in  Fig.  1.  It  should  be  pointed  out  that 
peak  B  is  clearly  seen  from  the  data  obtained  at  zero  electric 
field  after  application  of  a  dc  electric  field,  as  shown  in  Fig. 
2,  which  is  different  from  that  of  the  virgin  sample  in  Fig.  1. 
This  implies  the  effect  of  the  electrical  history  on  the  sample. 
In  addition,  peak  A  around  100  K  is  also  more  clearly  seen 
under  dc  bias. 

Both  peaks  A  and  B  are  also  observed  in  other  PZN-PT 
samples  with  different  PT  contents,  such  as  8%  PT  and  12% 
PT,  indicating  a  reproducible  behavior  for  the  PZN-PT  sys¬ 
tem. 

The  unipolar  strain  (S)  versus  electric  field  (£)  behavior 
for  PZT-4.5  PT  single  crystals  measured  at  cryogenic  tem¬ 
peratures  and  0.2  Hz  is  shown  in  Fig.  3.  The  inset  shows  the 
temperature  dependence  of  the  strain  level  at  60  kV/cm.  At 
room  temperature,  the  strain  level  of  ~  1  %  is  obtained  at  60 
kV/cm,  and  a  field-induced  phase  transition  is  observed.  All 
these  are  consistent  with  the  results  reported  in  the 
literature.2,3  With  decreasing  temperature,  the  strain  is  de¬ 
creased,  and  the  field-induced  phase  transition  gradually  nar¬ 
rowed,  disappearing  at  — 190  K.  At  a  temperature  as  low  as 
110  K,  the  strain  is  —0.45%  at  60  kV/cm,  which  is  still 
relatively  high,  indicating  excellent  potential  of  this  material 
for  practical  applications  at  cryogenic  temperatures. 

From  the  dielectric  spectra  obtained,  the  major  dielectric 
characteristics  for  the  PZN— 4.5  PT  samples  can  be  summa¬ 
rized  as  follows, 

(1)  Peak  A  around  100  K  with  obvious  frequency  dispersion 
is  observed.  The  relaxation  time  as  a  function  of  tem¬ 
perature  follows  the  Arrhenius  law  with  the  activation 
energy  t/=0.24eV  and  ro  =  ~1.0X  10" 15  s. 

(2)  There  is  a  peak  B  around  250  K,  which  can  be  seen  more 
clearly  under  dc  bias,  whose  peak  temperature  is  inde¬ 
pendent  of  dc  bias,  although  its  intensity  is  suppressed 
by  the  dc  bias. 


Peak  B  seems  to  be  a  more  “intrinsic”  behavior.  It  is 

FIG.  2.  Temperature  dependence  of  e  and  !S  reasonable  to  assume  it  exists  even  without  dc  bias.  It  shows 

tals  measured  at  10  kHz  under  dc  bias  field  0,  5,  and  10  kV/cm  in  a  cooling  •  . 

up  clearly  under  dc  bias,  because  the  dc  bias  greatly  sup- 
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presses  the  contribution  from  other  polarization  mechanisms 
in  the  sample.  In  addition,  it  is  also  found  that  the  tempera¬ 
ture  of  this  peak  is  independent  of  dc  bias  levels.  This  is  very 
similar  to  the  Bi-induced  dielectric  anomalies  in  SrTi03 :  the 
so  called  “dielectric  relaxor”  behavior,  in  which  the  e  maxi¬ 
mum  temperature  does  not  shift  with  dc  field,  but  the  ampli¬ 
tude  is  suppressed.11  All  these  characteristics  support  that 
this  anomaly  may  relate  to  an  “intrinsic”  behavior.  It  is  remi¬ 
niscent  of  the  results  reported  by  Kamzina  and  Krainik10  for 
pure  PZN  single  crystals,  that  there  is  a  dielectric  anomaly 
around  250  K  induced  by  application  of  a  dc  bias,  which  was 
attributed  to  a  real  ferroelectric -paraelectric  phase  transition. 
A  clear  understanding  of  this  anomaly  needs  further  study. 

For  peak  A,  in  fact,  a  low-temperature  dielectric  relax¬ 
ation  with  the  activation  energy  in  the  range  0.2-0.5  eV 
below  200  K  has  been  observed  in  many  solid  solutions, 
such  as  PZT,'2  Ba(Ti,Ce)03,13  and  (Sr,Ba)Nb206,14  etc. 
Usually,  it  is  attributed  to  domain-wall  motions/relaxation  at 
low  temperature. 

From  the  temperature  dependence  of  the  dielectric  con¬ 
stant  and  loss  at  very  low  temperatures  (Fig.  1),  it  can  be 
seen  that  around  100  K  and  even  below,  the  dielectric  con¬ 
stant  is  still  high  and  keeps  decreasing,  indicating  that  the 
polarization  mechanism  that  contributes  to  high  strain  is  not 
completely  frozen  out  as  yet.  This  is  different  from  the  pi¬ 
ezoelectric  copolymers,  in  which  the  strain  lost  its  major,  and 
the  material  only  retains  less  than  1%  of  the  value  at  room 
temperature.  Correspondingly,  the  dielectric  constant  and 
loss  of  the  copolymer  are  already  frozen  out  below  150  K; 
that  is,  the  temperature  dependence  of  s  and  tan  S  just  is 
flattened  below  that  temperature. 15 

A  study  of  the  dielectric  behavior  for 
95.5%  Pb(Zn|/3Nb2/3)03-4.5%  PbTi03  single  crystals  with 
and  without  a  dc  electric  field  shows  two  sets  of  dielectric 
anomalies  at  cryogenic  temperatures.  Peak  A  is  observed 
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around  100  K.  The  relaxation  time  follows  the  Arrhenius  law 
with  r0  =  ~1.0X10“15s  and  U=  0.24  eV,  as  observed  in 
many  solid  solution  systems.  It  may  be  attributed  to  domain- 
wall  motions.  Peak  B,  occurring  around  250  K  at  10  kHz 
under  a  dc  electric  field,  is  related  to  a  more  intrinsic  behav¬ 
ior.  The  strain  level  is  above  ~0.45%,  at  60  kV/cm,  down  to 
110  K,  indicating  this  material  is  very  promising  for  appli¬ 
cations  in  space. 

This  work  was  supported  by  a  grant  from  Office  of  Na¬ 
val  Research  under  Contract  No.  N000 14-98- 1-0530. 
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Polarization  relaxation  was  studied  in  Pb(Zn1/3Nb2/3)03-PbTi03  (PZN-PT)  single  crystals  that  show 
fatigue  anisotropy.  To  excite  prepoled  crystals,  a  modest  dc  voltage  (<l/2  of  the  coercive  field)  was 
applied  along  the  poling  direction.  Upon  removal  of  the  voltage,  the  polarization  decay  in  the  time 
domain  was  measured.  Experimental  data  were  modeled  with  a  stretched  exponential  function. 

Stretching  exponent  and  characteristic  time  (r^kt))  constants  for  polarization  relaxation 

were  determined  from  data  over  four  decades  in  the  time  domain  at  different  stages  of  bipolar 
cycling.  P(hkl)  values  after  101  cycles  were  0.146±0.002  and  0.247±0.0004  in  the  (001)  and  (111) 
orientations,  respectively.  The  0(\n)  constant  increased  up  to  0.453±0.104  after  105  cycles  in  (111) 
oriented  crystals  that  show  fatigue.  However,  much  less  change  is  observed  in  0(ooi)  as  a  function 
of  cycling  for  (001)  crystals.  Characteristic  time  constants  for  relaxation  (r^*/))  were  calculated  for 
(001)  and  (111)  orientations  as  0.401  ±0.048  s  and  57.46 ±0.10  s,  respectively.  These  results  suggest 
a  fester  polarization  relaxation  in  (001)  than  in  the  (111)  orientation  of  rhombohedral  PZN-PT 
ferroelectric  crystals.  ©  2004  American  Institute  of  Physics .  [DOI:  10.1063/1.1641183] 


I.  INTRODUCTION 

The  Pb(Zn1/3Nb2/3)03-PbTi03  (PZN-PT)  system  is  a 
complete  solid  solution,  based  on  the  perovskite  crystal 
structure,  of  relaxor  (PZN)  and  normal  (PT)  ferroelectric 
compounds.  Since  the  discovery  of  their  extraordinary  piezo¬ 
electric  properties  (</33~2500pC/N)  with  large,  electric- 
field-induced  strain  ( — 1,7%)  in  rhombohedral  crystals  ori¬ 
ented  along  [001]c  (where  the  subscript  denotes 
measurements  with  respect  to  the  cubic  axes),  numerous  in¬ 
vestigations  have  been  made  in  this  system.1"8  It  has  been 
widely  noticed  that  the  measured  properties  of  PZN-PT 
single  crystals  are  strongly  anisotropic.  Park  et  al  explored 
the  dependence  of  electric-field-induced  strains  on  crystallo¬ 
graphic  orientation  in  PZN-PT  crystals.  Large  piezoelectric 
strains  with  extremely  low  hysteresis  were  observed  in 
[001]c  poled  crystals,  compared  to  lower  strain  with  large 
hysteresis  in  the  [lll]c  orientation.  The  very  high  strains, 
coupled  with  the  low  hysteresis  of  [001]c  oriented  crystals, 
were  attributed  to  the  engineered  domain  configuration  in 
which  the  domain  walls  are  not  driven  by  applied  fields  in 
the  poled  state.3  Electric-field-induced  domain  configurations 
and  domain  switching  under  pulsed  electric  field  conditions 
have  been  investigated  by  Yu  et  a/.8  in  PZN-PT  crystals  ori¬ 
ented  at  various  crystallographic  directions.  Both  domain 
wall  mobility  and  activation  energies  were  reported  to  be 
higher  in  the  [001]c  direction  than  for  [  1 1  l]c .  The  observed 
anisotropy  in  the  switching  under  pulse  conditions  was  be¬ 
lieved  to  indicate  that  the  engineered  domain  state  along  the 
[001]c  direction  is  more  stable  than  that  in  the  [  1 1  l]c  direc¬ 
tion.  Consequently,  a  higher  activation  field  is  needed  in 
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[001]c  oriented  crystals.  However,  once  the  energy  barrier  is 
overcome  during  the  switching  process,  domain  walls  move 
faster  in  [001  ]c  oriented  crystals.6,8 

Many  ferroelectric  materials  show  a  progressive  loss  in 
switchable  polarization  on  cycling  through  the  hysteresis 
loop  repeatedly.  This  process  is  known  as  polarization  fa¬ 
tigue.  Fatigue  studies  in  PZN-PT  crystals  also  revealed 
anisotropy.9"11  It  was  found  that  [001]c  oriented  rhombohe¬ 
dral  PZN-PT  demonstrates  excellent  fatigue  resistance,  while 
[  1 1  l]c  oriented  rhombohedral  PZN-PT  fatigues  as  shown  in 
Fig.  1.  Fatigue  is  also  observed  in  [001]c  crystals  when  te¬ 
tragonal  material  is  generated  by  changes  in  temperature  or 
driving  field  strength.11  The  superior  fatigue  resistance  along 
the  [001]c  direction  is  consistent  with  the  observation  of 
higher  domain  wall  mobility  (at  least  in  the  rhombohedral 
phase  field).8  It  is  important  to  note  that  in  some  studies, 
fatigue  has  been  correlated  with  a  slowing  down  of  domain 
walls.12 

There  is  likely  to  be  a  complex  interaction  between  ex¬ 
isting  sources  of  depolarizing  fields  under  switching  condi¬ 
tions,  and  the  evolution  of  domains  and  polarization  after 
external  field  is  removed  (polarization  relaxation  from  an 
excited  state).  One  way  in  which  this  could  be  manifested  is 
in  time  variation  in  the  remanent  polarization  due  to  sponta¬ 
neous  polarization  reversal.  This  is  sometimes  referred  to  as 
retention  loss,  and  is  widely  studied  in  thin  films.13-16  Some 
of  the  recent  studies  on  polarization  relaxation  focused  on 
the  physical  origins  of  the  relaxation  kinetics  during  polar¬ 
ization  reversal.  Different  retention  studies  covered  different 
time  frames,  from  very  short  times  ( t<  1  s)  to  long  times 
(/>  1  s).  The  main  concern  in  practical  ferroelectric  memory 
applications  is  the  long-time  retention.  However,  the  short- 
time  retention  is  also  important  in  probing  the  polarization 
dynamics  that  govern  the  fast  decay  of  polarization  upon 
removal  of  the  electric  field.  Kang  et  al.  investigated  the 
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FIO.  1.  Fatigue  anisotropy  in  PZN-4.5PT  single  crystal  for  <001)  and  <111) 
directions,  comparing  remanent  polarization  behavior  for  Nth  cycle  versus 
10  switched  case. 

short-time  retention  loss  properties  and  reported  that  fatigued 
capacitors  showed  a  significant  loss  in  retained  polarization 
in  the  short-time  regime.  This  was  interpreted  in  terms  of 
increasing  depolarization  fields  with  fatigue.17  Thus,  it  ap¬ 
pears  that  collective  information  on  the  polarization  dynam¬ 
ics  both  in  the  presence  and  absence  of  external  fields  may 
be  helpful  to  characterize  and  understand  the  polarization 
switching  and  fatigue  anisotropy. 

The  purpose  of  this  study  is  to  gain  insight  into  the  an¬ 
isotropy  in  PZN-PT  and  related  materials  by  probing  the  re¬ 
laxation  of  dc  field-excited  domains  upon  removal  of  the 
electric  field.  We  utilize  the  classical  time  decay  approach  to 
understand  the  effect  of  local  perturbations  on  polarization  as 
a  function  of  orientation  and  cycling  history.18-22  Here,  the 
focus  is  the  relaxation  kinetics  of  polarization  in  prepoled 
crystals  after  exposure  to  a  subcoercive  field  level  dc  field 
perturbation  applied  along  the  poling  direction.  The  polariza¬ 
tion  decay  in  the  time  domain  was  measured  for  two  differ¬ 
ent  crystallographic  orientations  in  PZN-4.5PT  single  crys¬ 
tals  that  demonstrate  different  fatigue  behavior,  and  in 
crystals  that  have  experienced  different  switching  histories. 
The  aim  of  the  work  was  to  determine  if  there  are  differences 
in  polarization  relaxation  in  [001] c  and  [lll]c,  as  well  as 
to  determine  how  the  relaxations  change  as  a  function  of 
bipolar  cycling. 

The  dynamic  processes  occurring  in  complex  systems 
such  as  the  PZN-4.5%PT  relaxor  ciystals  encompass  differ¬ 
ent  length  and  time  scales  for  the  polarization.  Both  fast  and 
slow  rearrangements  take  place  within  the  microscopic,  me¬ 
soscopic,  and  macroscopic  organization  of  the  systems.  A 
common  aspect  of  relaxor  or  glassy  ferroelectric  systems  is  a 
certain  amount  of  order  at  an  intermediate,  so-called  mesos¬ 
copic,  scale,  while  they  are  disordered  at  the  microscopic 
scale  and  are  homogeneous  at  the  macroscopic  scale.23  A 
simple  exponential  relaxation  law  cannot  describe  the  relax¬ 
ation  phenomena  and  kinetics  in  such  materials.  There  are 
two  mechanisms  driving  nonexponential  relaxation  in 
quenched,  disordered  media  with  stretched  exponential  re¬ 
laxation.  One  is  related  to  random  pinning  sites.24  Secondly, 


Ozgul  et  al. 

in  disordered  systems  such  as  spin  glasses,  complex  relax¬ 
ation  is  caused  by  the  existence  of  frustrated 
antiferromagnetic/ferromagnetic-type  interactions,25-27 
which  are  a  direct  consequence  of  the  quenched  disorder.28  A 
polarization  analogy  to  spin  glasses  has  been  made  in  dipole 
glasses.  In  such  materials,  the  nature  of  the  polarization  dy¬ 
namics  is  influenced  by  a  variety  of  perturbation  sources 
(e.g.,  micro-  and  mesoscopic  chemical  heterogeneities,  ran¬ 
dom  site  occupancy,  and  dipolar  interactions).  During  relax¬ 
ation  from  higher  energy  levels,  domains  will  relax  to  lower 
energy  configurations  by  multiple  pathways  influenced  by 
perturbation  sources.  This  study  is  intended  to  determine 
how  polarization  relaxation  correlates  with  the  orientation- 
dependent  fatigue  behavior  in  PZN-PT  crystals. 

II.  EXPERIMENTAL  PROCEDURE 

Single  crystals  of  0.955Pb(Zn1/3Nb2/3)03-0.045PbTi03 
(PZN-4.5PT)  solid  solutions  were  grown  by  a  high- 
temperature  flux  technique.29  The  crystals  of  this  composi¬ 
tion  have  a  perovskite  structure  and  are  rhombohedral 
(pseudocubic)  at  room  temperature.  The  crystals  were  ori¬ 
ented  either  along  the  £11 1]^  or  along  the  [001]c  axes  to 
within  ±2°  by  using  a  Laue  Camera  (Multiwire  Laboratories 
Ltd,,  real-time  Laue  machine). 

For  electrical  characterization,  plate-shaped  samples 
were  cut  from  the  oriented  crystals  and  prepared  by  polishing 
with  silicon  carbide  and  alumina  polishing  powders  (down  to 
—1  /zm)  to  achieve  flat  and  parallel  surfaces  onto  which 
silver  paste  electrodes  were  applied.  Silver  paste  electrodes 
were  preferred  because  they  can  be  removed  easily  by  wash¬ 
ing  with  acetone  without  changing  the  nature  of  the  crystal/ 
electrode  interface  after  the  electrode  removal.  The  thickness 
of  samples  used  in  this  study  ranged  from  450  to  600  /mi. 
Coercive  fields  (measured  at  20  kV/cm  and  0.1  Hz)  for  the 
crystals  are  about  3  and  5  kV/cm  in  (001)  and  (111)  direc¬ 
tions,  respectively. 

High-field  measurements  utilized  a  modified  Sawyer- 
Tower  circuit.  A  high-voltage  amplifier  (Trek  Model  609C-6) 
was  used  in  both  poling  and  polarization  fatigue  measure¬ 
ments.  The  experimental  procedure  for  the  relaxation  studies 
is  shown  in  Fig.  2.  Crystals  were  exposed  to  electric  fields 
(E)  as  high  as  5  Ec  for  poling  to  achieve  saturated  polariza¬ 
tion  prior  to  time  domain  experiments.  For  fatigue  measure¬ 
ments,  electric  fields  were  applied  with  a  triangular  bipolar 
wave  form.  The  magnitude  and  the  frequency  of  the  applied 
field  were  20  kV/cm  and  10  Hz,  respectively,  for  [001]c 
crystals.  The  applied  field  was  higher  (40  kV/cm)  for  [  11  lie- 
samples  to  ensure  full  switching  due  to  the  increasing  coer¬ 
cive  field  with  fatigue.  During  measurements,  the  samples 
were  submerged  in  Fluorinert,™  an  insulating  liquid,  to  pre¬ 
vent  arcing. 

To  investigate  the  influence  of  continuous  bipolar  cy¬ 
cling  on  domain  dynamics,  the  relaxation  in  the  time  domain 
was  probed  for  both  orientations  at  different  stages  in  the 
cycling  history:  1,  10,  102,  103,  104,  and  105  cycles.  That  is, 
fatigue  experiments  were  interrupted  at  these  numbers  of 
cycles  to  track  the  influence  of  cycling-induced  changes  on 
domain  dynamics.  To  maintain  as  constant  as  possible  a  do- 
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FIG.  2.  Experimental  procedure  for  polarization  relaxation  measurements  in  prepoled  PZN-4.5PT  single  crystals. 


main  state  for  comparison,  ciystals  were  always  repoled  fol¬ 
lowing  bipolar  cycling.  Crystals  exposed  to  a  few  number  of 
electrical  cycles  (i.e.,  up  to  10  cycles)  are  called  “virgin.”  In 
the  case  of  no  electric  field  cycling,  unipolar  electric-field- 
poled  crystals  are  called  “uncycled.” 


It  should  be  mentioned  here  that  the  main  concern  in  this 
study  is  to  probe  domain  relaxation  due  to  perturbation  of  the 
domain  structure  with  a  comparatively  weak  excitation  field. 
Thus,  it  was  important  to  separate  this  type  of  response  from 
the  primary  relaxation  that  occurs  in  crystals  following  the 


excited  domain  state  by  100  V  DC 


polarization  relaxation  to  lower  energy  state 
through  intermediate  pathways 

- > 

(b) 


FIG.  3.  A  schematic  of  polarization  states  during  (a)  high-field  poling  and  partial  depoling  and  (b)  low-field  excitation  and  relaxation  in  ferroelectric  crystals. 


Downloaded  02  Dec  2005  to  146.186.113.243.  Redistribution  subject  to  AIP  license  or  copyright,  see  http://jap.aip.org/jap/copyrlght.jsp 


2634  J.  Appl.  Phys.,  Vol.  95,  No.  5,  1  March  2004 


Ozgul  et  at. 


high-field  poling  process  schematically  depicted  in  Fig.  3(a). 
When  a  field  much  higher  than  a  coercive  field  is  applied  to 
the  crystal,  domains  are  driven  to  a  high-eneigy  level  (2) 
from  the  lower  energy  level  (1).  Following  high-field  poling, 
crystals  will  experience  a  primary  relaxation  (partial  depol- 
ing)  back  to  a  lower  energy  state.  In  order  to  minimize  the 
complication  associated  with  these  two  different  processes, 
prepoled  crystals  were  given  at  least  12  h  to  stabilize  after 
the  primary  relaxation  process.  After  reaching  a  more  stable 
polarization  configuration,  a  weak,  constant  (—1.6  kV/cm 
depending  on  the  sample  thicknesses)  dc  field  was  applied 
for  about  2  min  in  the  poling  direction  to  excite  the  polariza¬ 
tion  to  a  higher  energy  state,  as  shown  in  Fig.  3(b).  After  2 
min,  the  voltage  source  was  turned  off  and  the  current  decay 
due  to  the  domain  relaxation  was  measured  by  a  HP4140B 
pA-meter/dc-voItage  source.  This  setup  was  able  to  measure 
the  current  after  about  0.5  s.  In  this  study,  current  was  moni¬ 
tored  over  four  decades  in  the  time  domain.  In  this  case,  the 
polarization  is  excited  from  a  lower  energy  state  (a)  to  a 
higher  energy  level  (b)  by  a  low  dc  field.  Upon  removal  of 
the  field,  domains  relax  back  to  lower  the  crystal’s  total  free 
energy  through  intermediate  pathways  via  domain  wall  mo¬ 
tion  as  well  as  mesoscopic  and  atomic  scale  rearrangements. 
These  experiments  were  performed  in  both  orientations  to 
relate  the  fatigue  anisotropy  to  domain  dynamics.  The  mea¬ 
sured  currents  in  the  time  domain  during  relaxation  reflect 
the  nature  of  the  domains  in  two  distinct  orientations. 

The  measured  current  versus  time  (I-t)  data  were  con¬ 
verted  to  polarization  versus  time  (P-t).  The  time  depen¬ 
dence  of  the  polarization  P(t)  was  fitted  to  a  stretched  ex¬ 
ponential  function  to  extract  the  constants  which  give 
information  about  the  relative  polarization  relaxation  in  both 
directions  as  a  function  of  the  cycling  history  (fatigue). 

Polarization  relaxation  experiments  are  very  dependent 
on  the  measurement  conditions  and  techniques  utilized  to 
probe  them.  Time-dependent  processes  can  be  influenced  by 
temperature,  the  total  time  over  which  the  properties  probed, 
and  the  magnitude  of  the  perturbation  used  to  make  the  mea¬ 
surements,  Consequently,  the  results  reported  here  corre¬ 
spond  to  the  specific  measurement  conditions  described. 

111.  RESULTS  AND  DISCUSSION 

A  typical  time  dependence  for  the  measured  current 
(I-t),  along  with  the  converted  polarization  data  P(t ),  mea¬ 
sured  at  room  temperature,  are  presented  in  Fig.  4.  In  this 
case,  a  1.6  kV/cm  (<  1/2  of  Ec)  constant  dc  field  acts  as  the 
driving  force  for  the  polarization  excitation.  The  polarization 
relaxation  (from  an  excited  state)  occurs  when  the  external 
field  is  removed,  as  depicted  in  Fig.  3(b).  The  measured  po¬ 
larization  decay  in  the  time  domain  (P-t)  data  was  fitted  to 
a  well-known  phenomenological  relaxation  function;  namely, 
the  Kohlrausch-Williams-Watts3031  (KWW)  function 

P(t)-Pi-P2^(t/r)fi],  (1) 

where  r  is  a  characteristic  relaxation  time,  and  (3  is  the 
stretched-exponent  coefficient  ranging  from  0  to  1.  Recently, 
this  function  has  been  used  to  analyze  different  types  of  re¬ 
laxation  data  in  the  time  domain.14,15,32  It  has  also  been  sug- 


(b)  Time  (s') 

FIG.  4.  Polarization  relaxation  probed  by  (a)  dc  measurements  in  the  time 
domain  and  (b)  calculated  decrease  in  polarization,  (time  scale  is  linear). 


gested  that  there  is  a  direct  connection  between  the  stretched 
exponential  behavior  and  the  fractal  nature  of  the  disordered 
material.33  In  recent  years,  dynamical  heterogeneity  was  es¬ 
tablished  as  a  possible  source  of  stretched  exponential 
relaxation.21  The  components  of  the  relaxation  response  in 
most  materials  are  related  to  the  size  distribution  of  relaxing 
entities  (domains).21  A  droplet  model  was  proposed  to  de¬ 
scribe  the  existence  of  polarization  irregularities  (i.e.,  polar 
droplets)  within  the  poled  state.34-36  Priya  et  al  have  studied 
the  mechanical  aging  behavior  in  relaxor-PT  crystals  close  to 
the  morphotropic  phase  boundary,  and  concluded  that  poled 
single  crystals  may  not  have  simple  domain  states,  but  rather 
may  have  many  irregular  small  polar  regions  within 
domains.37  These  domain  structures  are  common  in  relaxor 
materials  and  solid  solutions,  as  previously  reported.38,39 

A  typical  KWW  fit  of  experimental  P(t)  data  for  (001) 
and  (111)  orientation  is  shown  (the  solid  line  is  the  fitting, 
the  symbols  are  the  experimental  data)  in  Fig.  5.  For  both 
orientations,  the  magnitude  of  the  polarization  relaxation  de¬ 
tected  is  modest,  at  most  a  few  tenths  of  a  jnC/cm2  when  the 
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FIG.  5.  Probed  current  decay  in  time  (a)  and  calculated  polarization  relax¬ 
ation  (b).  (time  scale  is  logarithmic). 


fit  is  extrapolated  back  to  t~ Os.  Thus,  the  polarization 
changes  that  are  being  probed  in  this  measurement  do,  in 
fact,  correspond  to  small  perturbations,  rather  than  to  signifi¬ 
cant  rewriting  of  the  domain  state.  Furthermore,  a  relatively 
small  difference  in  the  magnitude  of  the  polarization  relax¬ 
ation  current  was  detected  in  (111)  oriented  crystals  before 
and  after  the  fatigue  process  (the  amount  of  relaxing  polar¬ 
ization  decreases  by  about  a  factor  of  two  after  105  cycles). 
The  fitting  parameters  ($  and  r  were  determined  from  the 
KWW  fits  as  a  function  of  the  number  of  bipolar  switching 
cycles  and  plotted  in  Fig.  6  for  both  orientations.  As  shown 
in  Fig.  6(a),  the  stretching  exponent  {$)  values  are  close  for 
both  orientations  in  virgin  crystals:  0.146±0.002  for  (001) 
and  0.247 ±0.0004  for  (111)  after  10  cycles.  However,  a  re¬ 
markable  increase  was  recorded  in  P(m)  to  0,453±0.104 
after  105  switching  cycles.  An  average,  low  /?<ooi)  value 
(0.1 545  ±0.0085)  was  maintained  throughout  the  cycling  ex¬ 
periments.  Greatly  different  characteristic  time  constants  (r) 
were  determined  for  the  (001)  and  (111)  orientations  as  0.401 
±0.048  and  57.46±0.10  s,  respectively.  There  is  not  much 
change  in  the  relaxation  time  constant  values  during  cycling, 
regardless  of  fatigue.  This  may  be  because  only  mobile  do¬ 
mains  contribute  to  the  measured  relaxation  data  in  these 
conditions.  These  results  are  consistent  with  both  the  fatigue 
anisotropy  observation  shown  in  Fig.  1  and  the  literature 


FIG.  6.  Polarization  relaxation  data  fit  to  the  stretched  exponential  function 
(KWW)  for  (a)  (001);  (b)  (111)  orientations,  (linear  time  scale). 


relating  fatigue  to  the  frozen  islands  of  nonswitching 
regions.12,40,41  As  fatigue  evolves  in  the  (111)  oriented  crys¬ 
tal,  it  seems  that  the  slow-moving  domain  walls  are  gradu¬ 
ally  frozen  out  and  that  their  contribution  to  the  total  switch¬ 
ing  processes  diminishes,  reducing  the  total  effective  volume 
of  switchable  ferroelectric  crystal.  (Ill)  oriented  crystals 
may  be  more  susceptible  to  fatigue  because  of  the  higher  r 
values.  An  increase  in  corresponds  to  a  less  stretched 
relaxation,  and  this  is  observed  in  the  fatiguing  orientation. 
On  the  other  hand,  the  low,  stable  /?<ooi)  values  for  (001) 
crystals  are  a  reflection  of  a  very  broad  time  distribution  for 
switching.  These  analyses  are  summarized  schematically  in 
Fig.  7,  where  g(r)  is  the  distribution  function  for  time  con¬ 
stants  r. 

Comparable  fitting  constants  were  reported  in  other 
ferroelectric  systems  for  long-time  relaxation 
experiments.15,17  Gruverman  et  al.15  reported  a  stretching  ex¬ 
ponent  in  tetragonal  Pb(Zr02Tio8)03  (PZT)  films  of  0.24. 
Kang  etal}1  investigated  short-time  retention  loss  in  PZT 
films  before  and  after  fatigue.  The  decay  of  retained  polar¬ 
ization  A P(t)  was  analyzed  in  terms  of  polarization  relax¬ 
ation  with  a  distribution  of  relaxation  time,  and  was  found  to 
follow  a  power  law  APCO-APoo+APofl+^/to)]-”, 
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FIG.  7.  Polarization  relaxation  data  fit  to  the  stretched  exponential  function 
(KWW)  for  (a)  <001);  (b)  (111)  orientations,  (logarithmic  time  scale). 


where  A and  (AP^+APq)  are  the  asymptotic  values  of 
A P  at  /=<»  and  0,  respectively,  and  t0  is  a  characteristic 
relaxation  time.  The  power  law  exponent  n  values  were  0.07 
and  0.23 ±0.02  before  and  after  fatigue,  respectively.  The 
relaxation  parameters  obtained  in  this  study  seem  to  be  in  a 
reasonable  range,  and  an  increase  with  fatigue  in  stretching 
exponent  is  consistent  with  the  study  by  Kang  et  al 

It  is  difficult  to  make  a  sensible  comparison  of  r  values 
with  data  from  the  literature,  since  relaxation  processes  in 
ferroelectrics  take  place  over  such  wide  time  scales  (r  values 
ranging  over  more  than  10  orders  of  magnitude  in  time  have 
been  reported)  and  the  field  excitation  is  applied  under  dif¬ 
ferent  conditions.  Any  set  of  measurements  samples  only  a 
portion  of  the  distribution.  Furthermore,  it  is  not  clear  that 
the  relaxation  times  associated  with  spontaneous  polarization 
reversal  (which  are  presumably  governed  by  net  internal 
fields)  will  be  identical  to  those  associated  with  polarization 
decay  after  weak  excitation,  or  those  associated  with  charac¬ 
teristic  time  constants  for  polarization  reversal  under  applied 
fields.  Consequently,  comparisons  should  only  be  made  for  r 
values  of  samples  exposed  to  the  same  type  of  excitation. 
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(a)  Number  of  cycles 


FIG.  8.  Stretched  exponential  constants  with  the  number  of  cycles  for  (001) 
and  (111)  orientations:  (a)  stretching  exponent  &  and  (b)  characteristic  time 
for  relaxation  r. 


Although  it  is  difficult  to  make  any  comparison  of  the  r 
values  measured  here  to  those  measured  elsewhere,  an  inter¬ 
nal  comparison  of  and  r<1H)  measured  under  the  same 
conditions  is  revealing  (see  Fig.  8),  It  is  clear  that  is 

<111> 


X  <001 >  X  <11 1>  T 


FIG.  9.  A  schematic  depiction  of  time  constant  distribution  function  for 
relaxation  in  (001)  and  (111)  orientations  based  on  fitting  data. 
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appreciably  smaller  than  which  is  consistent  with 

prior  reports  of  higher  domain  wall  mobilities  in  PZN-PT 
crystals.6,8 

The  major  difference  in  relaxation  response  of  crystals 
with  different  orientations  could  be  related  to  the  different 
domain  configurations  adopted.  In  poled  rhombohedral  crys¬ 
tals,  the  spontaneous  polarization  is  along  the  {1 1 1)  direction. 
Thus,  at  least  one  of  the  polarization  vectors  is  parallel  to  the 
poling  direction  in  (111)  oriented  crystals.  This  state  is 
sampled  during  the  (initial)  switching  and  poling  experi¬ 
ments,  as  evidenced  by  the  large  polarizations  observed  for 
(111)  crystals  (see  Fig.  9).  The  instability  of  this  domain  state 
is,  however,  widely  discussed,  and  at  lower  fields,  the  sample 
drops  back  to  a  multidomain  configuration.  Thus,  it  is  likely 
that  in  the  polarization  relaxation  experiments,  excitation  of 
the  polarization  is  from  the  three  canted  polarization  direc¬ 
tions  with  one  component  parallel  to  [111].  Note  that  the 
magnitudes  of  the  polarization  currents  in  this  case  are  small 
(comparable  to  those  of  the  (001)  crystals),  suggesting  that 
the  polarization  relaxation  measurement  is  not  probing  a 
switching/depoling  process.  It  is  also  interesting  that  there  is 
now  experimental  evidence  that  (111)  crystals  switch 
through  the  inclined  axes,  rather  than  by  180°  switching.42 
This  would  suggest  that  as  the  fatigue  process  proceeds,  the 
sample  has  been  cycled  through  states  in  which  domains 
continually  appear  and  disappear,  presumably  by  nucleation 
and  growth  processes. 

Pinning  of  domain  walls  in  this  case  would  be  consistent 
with  the  observed  fatigue  and  the  increase  in  since 

the  exponent  p  is  characteristic  of  the  domain  size  and  shape 
distribution  in  the  material.17 

In  contrast,  all  of  the  polarization  vectors  are  inclined  to 
the  poling  direction  in  (001)  crystals.  For  this  orientation, 
large  concentrations  of  domain  walls  (many  of  them 
charged)  are  built  in  during  the  poling  process,  and  the  re¬ 
sulting  domain  states  are  known  to  be  stable.  The  domain 
structures  in  (001)  and  (111)  orientations  are  known  to  differ. 
The  mechanics  of  the  polarization  reversal  process  are  not 
known  for  (001)  crystals,  but  in  any  event,  the  large  existing 
concentrations  of  domain  walls  (as  well  as  the  presence  of  a 
mechanism  for  absorbing  and  desorbing  charges  at  the 
charged  domain  walls)  may  be  related  to  the  increase  in  fa¬ 
tigue  resistance. 

IV.  CONCLUSIONS 

Single-crystal  PZN-PT  exhibits  a  stretched  exponential- 
type  relaxation  in  the  time  domain  at  room  temperature.  Re¬ 
laxation  data  fit  to  a  stretched  exponential  function  (KWW) 
with  small  ft  values  (<0.25),  implying  a  relatively  broad 
distribution  of  relaxation  times.  The  crystallographic  orienta¬ 
tion  dependence  of  the  slow  relaxation  process  has  been 
studied  as  a  function  of  bipolar  switching  cycles  (fatigue). 
Stretched  exponential  analysis  of  the  experimental  results  for 
(001)  and  (111)  orientations  at  different  levels  of  cycling 
demonstrated  an  orientation  dependence  of  the  stretching  ex¬ 
ponent  p  and  characteristic  relaxation  time  r.  Small  values 
for  P  and  r  were  determined  for  the  (001)  orientation,  and 
these  remained  relatively  stable  throughout  the  cycling. 


However,  P  increased  with  fatigue  for  the  (111)  orientation. 
The  characteristic  time  constant  r  is  much  higher  for  (111) 
crystals.  Both  r<Ui)  and  j>  are  independent  of  fatigue. 
These  observations  would  be  consistent  with  a  picture  in 
which  the  slow  domains  freeze  out  with  fatigue,  so  that  their 
contribution  to  the  relaxation  diminishes,  for  (111)  oriented 
crystals.  The  polarization  relaxation  is  much  faster  in  (001) 
oriented  PZN-PT.  Constant  and  relatively  small  £<ooi>  an(* 
r<ooi)  values  imply  a  broader  distribution  of  relaxation  pro¬ 
cesses  in  (001)  crystals,  peaked  at  a  shorter  time  than  for  the 
(111)  orientation.  Additional  work  to  clarify  the  relation  be¬ 
tween  the  polarization  relaxation  kinetics  and  differences  in 
the  domain  structures  is  needed  to  clarify  the  impact  of  ki¬ 
netics  on  the  fatigue  process. 
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Direct  evidence  of  ferroelastic  participation  in  180°  polarization  switching 
and  fatigue  for  111  oriented  rhombohedral  ferroelectric  0.955 
Pb(Zn1/3Nb2/3)03: 0.045  PbTi03  single  crystals 
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Direct  elastic  deformation  measurements  taken  during  polarization  reversal  in  111  oriented 
0.955  Pb(Znj/3Nb2/3)03-  0.045  PbTi03  single  crystals  confirm  the  observation  by  W.  Cao 
[Ferroelectrics  290,  107  (2003)]  that  a  ferroelastic  reorientation  of  the  domain  polarization  vectors 
is  involved,  and  the  observed  strain  changes  are  consistent  with  the  model  he  proposed.  First  cycle 
nonrecoverable  fatigue  occurring  in  our  crystals  and  the  associated  changes  of  shape  in  polarization 
and  strain  hysteresis  are  suggested  to  be  due  to  microcracking  associated  with  the  large  very  rapid 
strain  changes  associate  with  the  ferroelastic  switching  mode.  ©  2004  American  Institute  of 
Physics.  [DOI:  10.1063/1.1690869] 


Single  crystals  in  the  ferroelectric  rhombohedral  struc¬ 
ture  phase  of  the  lead  zinc  niobate:  lead  titanate  (PZN:PT) 
solid  solution  family  show  outstanding  but  highly  anisotropic 
piezoelectric  properties.  Samples  poled  along  the  prototypic 
001  direction  can  have  piezo-coupling  k33  up  to  95%  and 
d33~2500  pC/N. 1,2  Surprisingly  for  carefully  poled  samples 
at  composition  not  too  close  to  the  rhombohedralrtetragonal 
phase  boundary,  the  behavior  is  fully  intrinsic  and  associated 
with  easy  field  induced  tilting  of  the  macroscopic  Ps  vector 
in  the  poled  domain  away  from  the  equilibrium  (111)  orien¬ 
tation.  A  model  now  confirmed  by  phenomenological  calcu¬ 
lation  from  prototypic  Q{j 3  optical  observation4  and  first 
principles  ah  initio  calculation.5  As  would  be  expected  on 
this  model,  111  field  poled  crystals  at  the  same  composition 
have  much  more  modest  piezo  properties. 

Ferroelectric  switching  for  both  001  and  111  oriented 
fields  shows  excellent  “square”  dielectric  hysteresis  with 
similar  low  coercive  fields.  For  the  001  oriented  fields  both 
single  crystals6  and  epitaxial  oriented  thin  films7  show  re¬ 
markable  fatigue  free  switching  which  is  unexpected  as  the 
remanent  state  is  necessarily  polydomain.  Further,  for  111 
oriented  fields  where  the  crystal  can  be  driven  to  a  single 
macrodomain  state  there  is  evidence  of  strong  fatigue.  For 
this  111  oriented  switching,  however,  recently  Cao  et  al.s 
have  shown  from  acoustic  measurements  that  switching  is 
not  by  a  simple  180°  domain  wall  process,  but  most  likely 
involves  a  Ferroelastic  wall  motion. 

In  this  study  carried  out  on  (0.955)  Pb(Zni/3Nb2/3)03 
(0.045)  PbTi03  crystals  we  have  measured  both  longitudinal 
and  transverse  strain  and  dielectric  hysteresis  for  1 80°  polar¬ 
ization  switching  in  the  111  oriented  crystal.  The  results  are 
consistent  with  the  model  suggested  by  Cao  and  indicate  that 
reorientation  from  111  to  FFF  is  via  a  ferroelastic  mecha¬ 
nism.  Fatigue  studies  on  repeated  switching  also  suggest  that 
in  our  crystals  there  is  a  rapid  initial  damage  which  is  not 
recoverable  by  annealing  above  the  Curie  temperature  fol- 
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lowed  by  a  slower  more  substantial  fatigue  which  is  recov¬ 
erable  by  annealing.  Both  elastic  and  dielectric  response  are 
consistent  with  the  probability  that  the  unrecoverable  fatigue 
involves  a  microcracking  under  the  stress  concentrations  as¬ 
sociated  with  the  ferroelectric  component  of  switching.  The 
precise  mechanism  for  the  preferred  ferroelectric  trajectory  is 
not  clear;  however  the  massive  transverse  permittivity  and 
the  adaptive  state  models  of  Jin  et  al.9  underscore  the  pro¬ 
pensity  for  rotational  instability  of  the  domain  polarization 
vectors. 

Crystal  of  PZN:PT  at  the  0.955  Pb(Zn1/3Nb2/3)03  0.045 
PbTi03  composition  were  supplied  by  TRS  State  College. 
The  ferroelectric  Curie  temperature  of  162  °C  agreed  well 
with  that  expected  for  the  4.5%  PT  composition.  Samples 
with  111,  OH  and  112  orthogonal  faces  had  dimensions  in 
the  range  5X5X1  mm  (larger  111  face)  to  2X2X2  mm. 
Gold  electrodes  were  applied  to  the  111  faces  by  sputter 
coating.  Dielectric  hysteresis  was  measured  on  a  modified 
Sawyer-Tower  circuit  and  associated  strain  measured  using 
a  linear  variable  displacement  transducer  (LVDT)  driven  by 
a  lock-in  amplifier  (Stanford  Research  system  model  SR 
830).  Cyclic  triangular  bipolar  waves  at  0.2  Hz  with  field 
levels  up  to  10  kV/cm  were  supplied  from  a  Trek  609  C-6 
high  voltage  amplifier.  For  longitudinal  (33)  strain  measure¬ 
ment  the  LVDT’s  insulating  probes  contacted  directly  onto 
the  111  face  electrodes.  In  transverse  (11,22)  strain  measure¬ 
ments  the  sample  was  now  mounted  with  1  F0  or  112  faces 
contacted  by  LVDT  probes  and  electric  field  supplied  by  fine 
silver  wires  attached  by  air  drying  silver  paint  to  the  elec- 
troded  faces.  All  measurements  were  carried  through  in  an 
air  conditioned  room  at  20  °C. 

Initial  hysteresis  on  first  cycling  of  any  of  the  111  ori¬ 
ented  crystal  tested  showed  excellent  near  square  hysteresis 
as  in  Fig.  1(a)  with  an  abrupt  transition  into  the  poled  con¬ 
dition  for  both  positive  and  negative  switches.  Associated 
longitudinal  ferroelastic  hysteresis  is  evidenced  in  Fig.  1(b) 
carried  out  at  the  same  time  and  on  the  same  sample.  It  is 
interesting  to  note  that  the  peaks  in  the  longitudinal  strain 
occur  on  the  steeply  falling  part  of  the  polarization  curve. 
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FIG.  1 .  (Color)  (a)  First  cycle  of  hysteresis  in  polarization  vs  field  for  a  111 
oriented  0.955  Pb(Zn1/3Nb2/3)03 :  0.045  PbTi03  single  crystal  subjected  to  a 
triangular  bipolar  wave  of  10  kV/cm  at  0.2  Hz.  (b)  Longitudinal  strain  S33 
during  first  cycle  switchover  from  111  oriented  Ps  to  TIT  oriented  Ps . 

longitudinal  strain  being  all  negative.  Unfortunately,  with  our 
simple  setup  it  is  not  possible  to  measure  both  longitudinal 
and  transverse  strain  on  the  same  sample  at  the  same  time. 
For  the  crystal  for  which  data  are  presented  in  Fig.  2,  trans¬ 
verse  then  longitudinal  strain  were  measured  in  sequence 
[Fig.  2(b)]  with  the  polarization  loop  [Fig.  2(a)]  correspond¬ 
ing  to  the  second  cycle.  Again  the  strain  changes  are  nicely 
abrupt  but  there  is  a  minor  rounding  associated  with  the  ini¬ 
tial  fatigue  mechanism.  As  in  all  111  samples  measured  the 
transverse  strain  occurring  during  switching  is  all  positive 
but  the  maximum  transverse  strain  is  less  than  half  the  mag¬ 
nitude  of  the  maximum  longitudinal  change. 

In  a  number  of  crystals  tested  the  initial  shape  change 
both  in  dielectric  and  elastic  hysteresis  was  more  severe  as  in 
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FIG.  3.  (Color)  (a)  Early  cycle  of  hysteresis  in  111  oriented  0.955 
Pb(Zn1/3Nb2/3)03 : 0.045  PbTi03  single  crystal  under  10  kV/cm  triangular 
bipolar  wave  at  0.2  Hz.  (b)  Early  cycle  transverse  and  longitudinal  strain 
under  0.2  Hz  cyclic  triangular  wave  field  L  is  longitudinal  strain  Tx  Oil 
transverse  strain  T2  112  transverse  strain. 

Fig.  3  where  the  field  required  for  complete  switchover  is 
almost  doubled  and  dominantly  the  transverse  dilatation  per¬ 
sists  to  these  higher  field  levels.  As  has  been  reported  earlier7 
repeated  cycling  leads  to  more  severe  fatigue  and  loss  of 
remanent  polarization.  For  our  crystals  only  200  cycles  at  0.2 
Hz  with  peak  triangular  bipolar  wave  of  10  kV/cm  lowers 
remanent  polarization  from  42.2  to  —30  ficlcm2  and  raises 
coercivity  from  3.95  to  5.5  kV/cm.  In  agreement  with  the 
earlier  studies,  the  more  severe  fatigue  is  recoverable  by  an¬ 
nealing  at  200  °C  which  is  above  Tc ,  but  the  polarization 
does  not  “square  up”  and  the  curvature  of  both  longitudinal 
and  transverse  strain  persists.  It  has  been  suggested  by  Cao8 


E(kWcm) 


FIG.  2.  (Color)  (a)  Second  cycle  hysteresis  in  polarization  vs  field  for  a  Ill 
oriented  0.955  Pb(Zn,/3Nb2/3)03 : 0.045  PbTiOj  single  crystal  subjected  to  a 
triangular  bipolar  wave  of  10  kV/cm  amplitude  at  0.2  Hz.  (b)  Second  cycle 
longitudinal  strain  (blue)  and  first  cycle  transverse  strain  (red)  in  a  111 
oriented  0.955  PWZn^Nb^Oa  : 0.045  PbTi03  crystal  subjected  to  cycling 
under  a  triangular  bipolar  wave  of  10  kV/cm  of  amplitude  at  0.2  Hz. 


FIG.  4.  (Color)  (a)  Relationship  between  rhombohedral  and  prototypic  cubic 
coordinates  and  labeling  of  the  eight  possible  polar  directions  in  rhombohe¬ 
dral  domain  state  (taken  from  Cao  in  Ref.  18)  Proposed  switching  path  in 
111  oriented  0.955  Pb(Zn1/3Nb2/3)03 : 0.045  PbTi03 .  State  1  domains  may 
switch  to  any  one  of  states  2,  3,  4  before  switching  over  to  state  8  to 
complete  the  reversal  process,  (b)  Orientation  of  the  spontaneous  strain  el¬ 
lipsoids  during  the  switching  process  indicating  the  origin  of  Aj  A3  which 
contribute  the  large  contraction  along  111  during  switching,  (c)  Plan  view 
normal  to  111  indicating  the  origin  of  the  expansions  A4  A5  which  contrib¬ 
ute  to  the  transverse  elongation.  Note:  Eccentricities  have  been  grossly  ex¬ 
aggerated  for  illustrative  purposes. 


W.  Zhu  and  L.  E.  Cross 


••*390  Appl.  Phys.  U  Vol.  84,  No.  13,  29  March  2004 

that  for  PZN:4.5  PT  180°  switching  for  111  oriented  crystals 
is  a^two  stage  process,  as  indicated  in  Fig.  4  taken  from  his 
article.  The  suggestion  is  that  the  polarization  vector  in 
single  domain  (1)  in  the  111  poled  state  “folds  out”  into 
domains  with  orientations  (2),  (3),  (4)  before  folding  over 
into  the  reverse  state7  with  III  orientation.  X-ray  measure¬ 
ments  of  the  ferroelectric  spontaneous  sheer  in  the  4.5  PT 
composites  give  a  value  at  the  sheer  angle  a~  89.89° 10  cor¬ 
responding  to  an  equivalent  elongation  along  1 1 1  in  the  Ref. 
8  direction 

AL 

—  =  A,  =0.0019. 

Since  under  pure  sheer  deformation  volume  is  conserved,  the 
transverse  strain 

AT 

A2=—  is  given  by 

1 

Ai  +  2A2=0  A2=-~A1  =  -0.00095. 

An  estimate  of  the  maximum  possible  longitudinal  strain 
at  switchover  may  be  taken  as  the  loss  of  the  state  1  elonga¬ 
tion  Ax ,  together  with  the  contraction  from  the  zero  strain 
state  associated  with  the  domain  deformation  ellipsoids  for  2, 
3,  and  4  states  [Fig.  4(a)]  now  being  tilted  19.47°  above  the 
plane  perpendicular  to  1 1 1  through  the  center  of  these  strain 
ellipsoids  [Fig.  4(b)].  This  yields  an  additional  contraction  of 

A3= -0.0008 

so  that  the  total  contraction  as  evidenced  in  Fig.  4(b)  is 

-  A  j  +  A3  =  “  0.0019-  0.0008  =  -  0.0027. 

For  the  transverse  strain,  the  tilted  strain  ellipsoid  for  the 
single  domain  gives  a  maximum  dilation  [Fig.  4(c)] 

A4  =  0.0016. 

However,  there  are  now  three  equivalent  (2),  (3),  (4)  do¬ 
mains  making  angles  of  120°  with  each  other.  The  smallest 
dilation  will  take  place  at  the  intersection  of  the  three  ellip¬ 
soids.  Taking  note  of  the  fact  that  the  circular  section  in  Fig. 
4(c)  is  not  that  in  Fig.  4(b),  but  corresponds  to  the  smaller 
diameter  circular  section  of  the  full  spontaneously  deformed 
single  domain  yield 

A5=0.000  15. 

Since  A4  and  A5  are  not  collinear  we  make  a  crude  approxi¬ 
mation  to  the  transverse  strain 

A4  +  A5\  (0.0016+0.000 15) 

1  =0.000  87. 


For  the  virgin  crystals  measured,  the  average  maximum  lon¬ 
gitudinal  strain  was 

A  L 

-jr-  =  —  0.0020 ±  0.000 1 

and  the  average  maximum  transverse  strain  was 
AT 

—  =  +  0.0008  ±.0001 

T 

suggesting  that  switchover  is  close  to  a  two  stage  process, 
and  lending  strong  support  for  the  switching  model  of  Cao. 

Perhaps  the  most  interesting  feature  of  the  fatigue  behav¬ 
ior  in  our  crystals  is  the  changeover  in  the  character  of  the 
loops  after  just  the  first  cycle.  The  most  severe  we  have  seen 
is  the  change  exhibited  in  Fig.  3  where  the  normally  very 
steeply  changing  completion  of  switchover  is  now  strongly 
delayed  and  particularly  the  transverse  strain  persists  to 
much  higher  field  levels  Fig.  3(b).  Since  the  longitudinal 
contraction  appears  to  complete  while  the  transverse  expan¬ 
sion  persists,  there  must  be  an  increase  of  volume  and  we 
suspect  microcracking.  If  some  of  these  cracks  make  an 
angle  with  the  E  field  direction  which  appears  probable,  the 
field  seen  by  the  crystal  in  that  region  will  be  reduced  and 
switchover  of  polarization  to  the  final  III  orientation  de¬ 
layed. 

It  is  not  unlikely  that  the  propensity  for  cracking  may 
depend  in  a  complex  manner  on  finishing,  electroding, 
mounting  of  the  crystals  and  on  internal  impedance  and  other 
details  of  the  electric  field  supply.  Differences  in  the  severity 
of  the  behavior  have  been  observed  as  indicated,  but  so  far 
we  have  not  observed  any  sample  in  our  experimental  setup 
which  does  not  show  some  initial  degradation  of  the  form 
described. 
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Abstract.  The  hydrostatic  piezoelectric  coefficient  dh  of  Pb(ZrxTii_x)03  ceramics  (PZT),  and  of 
Pb(Zni/3Nb2/3)C>3-PbTi03  and  Pb(Yhi/2Nbi/2)03-PbTi03  (PZN-PT,  PYN-PT,  respectively)  single  crystals  with 
compositions  near  to  the  morphotropic  phase  boundary  (MPB)  have  been  measured  using  a  dynamic  hydrostatic 
method.  The  effects  of  DC  electric  field  and  static  component  of  hydrostatic  stress  on  dh  of  PZT  ceramics,  PZN-PT 
and  PYN-PT  single  crystals  were  studied.  Changes  of  the  piezoelectric  hydrostatic  coefficients  dh  caused  by  an 
electric  field  (DC  bias)  were  observed  along  with  pressure  and  temperature  dependencies.  The  measurement  of 
the  hydrostatic  piezoelectric  coefficient  dh  seems  to  be  promising  for  investigation  of  intrinsic  (single  domain)  and 
extrinsic  (domain-walls)  contributions  to  piezoelectric  behavior  of  single  crystals  and  ceramic  materials. 
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1.  Introduction 

It  is  well  known  that  perovskite  solid  solutions  have 
been  intensively  investigated  because  of  their  excel¬ 
lent  piezoelectric  and  dielectric  properties.  They  con¬ 
sist  of  simple  cubic  structure  at  high  temperatures  and 
variety  of  phases  with  polar  states  at  lower  temper¬ 
atures.  Among  these  materials  are  the  ferroelectric 
Pb(ZrxTIi  -x)03  (PZT)  solid  solutions  that  present  such 
properties  in  the  morphotropic  phase  boundary  (MPB) 
where  rhombohedral  and  tetragonal  phases  coexist 
practically  independently  of  temperature.  Excellent 
piezoelectric  properties  were  found  for  relaxor-based 
single  crystals  such  as  Pb(Zni/3Nb2/3)03-PbTi03  and 
Pb(Ybi/2Nb1/2)03-Pbri03  (PZN-PT,  andPYN-PT,  re¬ 
spectively)  near  MPB  [1-4].  This  is  the  basis  for 
practically  all  transducers  and  other  piezoelectric 
devices. 
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Recently,  considerable  research  has  been  focused  on 
the  piezoelectric  and  parameters  and  dielectric 
properties  of  PZT  bulk  ceramics  and  thin-films.  The 
temperature  dependence  of  piezoelectric  coefficients 
of  PZT  ceramics  using  resonance  method  was  reported 
in  [5].  The  effect  of  hydrostatic  stress  on  piezoelec¬ 
tric  coefficients  of  PZT  was  published  in  [9].  Besides 
J33  and  t/31,  no  much  attention  was  given  to  study  the 
hydrostatic  piezoelectric  coefficient  dh  variation  under 
different  influences  (as  DC  electric  field,  temperature, 
pressure).  Thus,  in  our  previous  papers  [6-8],  the  meth¬ 
ods  for  determining  PZT  material  parameters  as  piezo¬ 
electric  coefficients  ^33,  d^ ,  and  hydrostatic  piezoelec¬ 
tric  coefficients  dh,  gh  were  described.  In  this  article 
our  work  deals  with  dh  measurements  by  a  dynamic 
hydrostatic  method  combining  a  static  and  variable  su¬ 
perimposed  pressure.  The  effect  of  a  DC  electric  field 
and  temperatures  on  the  hydrostatic  piezoelectric  co¬ 
efficients  dh  of  PZT  ceramics  as  well  as  relaxor-based 
single  crystals  as  PZN-PT  and  PYN-PT  were  investi¬ 
gated  and  described. 
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2.  Hydrostatic  Piezoelectric  Coefficients 

The  determination  of  the  dynamic  hydrostatic  coeffi¬ 
cients  utilizes  a  poled  sample  along  z-axes  under  hy¬ 
drostatic  pressure  pTu  =  —p,  all  other  components  of 
stress  tensor  equal  to  zero).  Simultaneously,  the  applied 
electric  field  £3  =  0,  i.e.  a  short-circuit  condition. 

The  polarisation  P3  on  the  sample  induced  by  hy¬ 
drostatic  pressure  p  [7],  is 

i>3  =  (d3 1  +  dyi  +  ^33)  •  (— p) 

=  (d33+2d3iM“£)  =  4  •(-/>),  (1) 

where  dik  are  piezoelectric  coefficients,  d^i  —  d32l  and 
dh  is  hydrostatic  coefficient. 

Measurement  of  the  induced  charge  Q  on  the  area 
A  of  the  electrodes  and  pressure  p  permits  calculation 
of  the  hydrostatic  piezoelectric  coefficient  <4. 

The  dynamic  method  was  described  in  [7]  by 
equation 

dO  1  dp 

—  T  =  dh  *  — ,  where  Q  -  go  sin  cot.  (2) 
dt  A  dt 

For  the  combined  static  and  dynamic  (with  the  si¬ 
nusoidal  shape)  pressure  excitation 

p  =  p$  sin  cot  +  p statio  (3) 

we  can  express  the  dynamic  hydrostatic  piezoelectric 
coefficient  dh  as: 


where  70  is  the  amplitude  of  the  current,  po  is  the  am¬ 
plitude  of  the  dynamic  pressure  excitation,  and  co  is  the 
angular  frequency. 

The  hydrostatic  piezoelectric  coefficient  dh  can  be 
also  calculated  from  the  Eq.  (1), 

dh  =  d33  4-  2d3i.  (5) 

Equation  (5)  was  used  for  calculation  of  dh  values 
in  our  previous  paper  [6]. 

3.  Experiment 

The  dynamic  method  of  dh  measurement  is  based  on 
dynamic  oscillations  of  pressure  around  an  appointed 


static  pressure  mastered  by  an  auxiliary  hydraulic  pis¬ 
ton  at  a  frequency  of  about  0.5-1  Hz.  The  temperature 
control  is  realized  by  the  PID  temperature  controller 
with  a  resistive  heater  and  by  the  compressor  cooler. 
The  advantage  of  the  described  method  is  the  elim¬ 
ination  of  pyroelectric  effect  and  higher  accuracy  of 
the  measurement.  The  possibility  of  an  applied  electric 
field  influence  on  the  device  under  test  was  realized. 
The  block  diagram  of  the  experimental  set  up  is  shown 
in  Fig.  1. 

We  studied  PZT  ceramics  and  PZN-PT,  PYN-PT 
crystals.  The  PZT  samples  EPC856  and  APC856  soft 
PZT  ceramics  with  composition  of  Pb(Zr58T42)03 
were  provided  by  American  Piezoceramics  Interna¬ 
tional,  Mackeyville,  PA,  USA.  These  samples  were 
machined  in  disk-shaped  (1  mm  thick  and  diameter  of 
10  mm)  and  square  plates  4x4x2  mm3. 

The  Pb(Zni/3Nb2/3)03-PbTi03  and  Pb(Ybi/2- 
Nbi/2)03-PbTi03  (PZN-PT,  PYN-PT,  respectively) 
single  crystals  were  grown  at  the  Pennsylvania  State 
University,  Material  Research  Institute,  University 
Park,  PA,  USA.  The  used  PZN-4.5%  PT  (001)  and 
(111)  cuts,  and  PYN-40%  PT  (1 1 1)  cut  were  machined 
in  rectangular  plate-shaped  samples  with  side  length 
from  1.3  mm  to  6.8  mm  and  thickness  from  0.4  mm  up 
to  0.6  mm.  The  measurements  of  the  dynamic  hydro¬ 
static  piezoelectric  coefficient  dh  were  performed  in  the 
temperature  range  from  283  K  to  433  K  at  a  static  pres¬ 
sure  of  10  MPa  and  a  dynamical  pressure  of  100  kPa. 

All  samples  were  electroded  on  both  major  surfaces 
and  poled  by  the  manufacturers. 


4.  Results  and  Discussion 

The  dependencies  of  the  hydrostatic  coefficient  dh 
of  the  soft  PZT  ceramics  (EPC856)  on  DC  electric 
field  measured  at  various  temperatures  are  displayed 
in  Fig.  2.  The  typical  relative  decrease  of  dh ,  caused 
by  DC  electric  field  change  from  zero  to  700  kVm"*1 
was  smaller  than  20%  at  temperatures  between  283  K 
and  423  K.  The  dh(E)  increases  to  vicinity  of  Curie 
temperature  (Tc  =  423  K). 

Figure  3  shows  the  similar  dependencies  of  dh  vs. 
DC  electric  field  and  temperature  for  the  soft  PZT  ce¬ 
ramics  (APC  856).  The  comparison  of  the  dh  values  ob¬ 
tained  from  resonance  method  and  dynamic  hydrostatic 
method  shows,  that  the  second  method  gives  slightly 
smaller  values  by  up  to  10  %.  The  details  are  shown  in 
the  Table  1. 
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i.  A  measuring  set  of  high  pressure  dynamic  method  with  the  possibility  of  application  of  electric  field  onto  sample. 


Fig.  2.  Hydrostatic  coefficient  dh  of  PZT  ceramics  (EPC  856)  vs. 
DC  electric  field  and  various  temperatures.  Disk-shaped,  10  mm  in 
diameter,  1  mm  in  thickness. 

The  temperature  dependencies  of  the  dh  coefficients 
at  different  DC  electric  field  values  applied  to  APC856 
ceramics  are  shown  in  Fig.  4.  A  noticeable  influence  of 
the  electric  field  on  Curie  temperature  Tq  position  was 
exhibited. 

Apparently,  the  phase  transition  is  shifted  to  higher 
temperature  by  increasing  the  electric  field.  This  is  due 
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Fig.  3.  Hydrostatic  coefficient  dh  of  PZT  ceramics  (APC  856)  vs.  DC 
electric  field  and  various  temperatures.  Plate-shaped,  4x4x2  mm3. 

to  the  upkeep  of  polar  state  by  the  applied  electric 
field  which  delays  the  phase  transition  ferroelectric- 
paraelectric  to  high  temperature. 

The  samples  of  PZN-PT  single  crystals  exhibit  small 
variation  of  dh  with  hydrostatic  pressure  in  the  pressure 
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Table  1.  Hydrostatic  piezoelectric  coefficients  dh  of  PZT  ceramics 
at  room  temperature  (APC  International,  Mackeyville,  PA,  USA.) 

PZT  PZT 

(APC856)  (EPC856) 


Resonance  method 

[10“12  C*N-1]  558 

-^31 

[10— 12  C-N”1]  236 

dh 

[10-12C-N-11  86 

Dynamic  hydrostatic  dh 
method  [10“12  C-N-1]  79 
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Fig.  4.  Temperature  dependencies  of  hydrostatic  coefficient  dh  of 
PZT  ceramics  (APC  856)  at  various  DC  electric  field.  Plate-shaped, 
4x4x2  mm3. 


591 

237 

117 

111 


range  of  zero  to  70  MPa.  Decrease  of  dh  value  does  not 
exceed  5%  of  initial  values  as  shown  in  Fig.  5. 

In  Fig.  6  is  plotted  the  hydrostatic  coefficient  dh  of 
PZN-4.5%PT  (001)  cut  vs.  DC  electric  field  at  various 
temperatures.  We  can  observe  a  change  of  the  shapes 
of  curves  at  temperatures  above  373  K  and  low  DC 
electric  fields.  We  suppose  that  it  is  connected  with  the 
phase  transition  from  ferroelectric  phase  to  paraelectric 
phase.  In  vicinity  of  Curie  temperature  Tc  =  423  K  the 
value  of  dh  is  very  low  at  zero  DC  bias.  For  low  electric 
field  values  we  can  observe  a  sharp  peak  of  the  coef¬ 
ficient  dh  and  its  consequent  exponential  decay  with 
increasing  bias.  The  peak  position  shifts  toward  higher 
applied  electric  field  values  with  increasing  tempera- 


p [10* Pa] 


Fig.  5.  Hydrostatic  coefficient  dh  of  PZN~4.5%PT  (001)  cut  and 
PZN-4.5%PT  (111)  cut  vs.  hydrostatic  pressure.  Plate-shaped, 
6.8  x  6.8  x  0.48  mm3,  3.2  x  3.2  x  0.6  mm3,  respectively. 
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Fig.  6.  Hydrostatic  coefficient  dj,  of  PZN-4.5%PT  (001)  cut  vs. 
DC  electric  field  at  various  temperatures.  Plate-shaped,  2.0  x  4.2  x 
0.46  mm3. 


tures.  Furthermore,  the  peak  magnitude  undergoes  a 
decrease  with  increasing  electric  field.  The  piezoelec¬ 
tric  hydrostatic  coefficient  peak  could  be  due  to  the 
domain  subdivision  which  takes  place  when  the  tem¬ 
peratures  approach  Tc*  Thus,  an  increase  of  temper¬ 
ature  engages  an  increase  in  the  domain  density.  The 
influence  of  electric  field  on  domains  (domain  reorien¬ 
tation)  will  depend  on  their  density,  and  consequently, 
an  increase  in  domain  density  will  need  a  high  elec¬ 
tric  field  for  domain  alignement.  This  could  explain 
the  peak  shift  to  high  electric  fields  when  the  tempera¬ 
ture  increases  [10, 11].  The  exponential  decay  of  right 


Hydrostatic  Piezoelectric  Coefficient  <4  of  PZT  Ceramics  447 


500 

- 

PZN:4v5%  PT  <111> 

-A-233K 

— V—333  K  ' 

▼  400 

_ 

— 0 —  2t3K 

403K  _ 

z 

-0—313K 

— *1 —  413  K 

d 

A  — I — 3S3K 

-0-423K  ’ 

~  300 

o 

- 

‘ 

■cf  200 

■j 

■ 

1UU 

4-  t—h  1  -4- 
%&&&&<> 

0 

T 

-i  j  i _ 

Ji - 1 - = 

0  500  1000  1500 


E  [  10*  V.m'1  ] 

Fig.  7.  Hydrostatic  coefficient  dt,  of  PZN-4.5%PT  (111)  cut  vs. 
DC  electric  field  at  various  temperatures.  Plate-shaped,  3.2  x  3.2  x 
0.59  mm3. 
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Fig.  8.  Hydrostatic  coefficient  df,  of  PYN-PT  (111)  cut  vs.  IX  elec¬ 
tric  field  at  various  temperatures.  Plate-shaped,  2.2  x  1.3  x  0.5  mm3. 


5.  Conclusion 


side  of  the  peak  could  be  associated  to  a  relaxation 
phenomenon  stimulated  by  charge  carrier  screening 
in  the  material  bulk  due  to  the  increase  conductivity 
with  temperatures.  It  was  shown,  in  our  experimen¬ 
tal  measurement,  the  value  of  dh  equaled  to  47*10~12 
CN"1  for  T  =s  298  K  and  £  =  0  Vn r1.  This  value 
is  smaller  than  dh  »  250* lO-"12  CN“l,  which  was  cal¬ 
culated  from  Eq.  (5)  using  measured  d33  and  d3 1  [1]. 
To  put  in  evidence  this  difference  further  investiga¬ 
tion  of  the  effect  of  a  DC  electric  field  on  domain 
engineering  states  and  domain  wall  configuration  is 
required. 

The  hydrostatic  coefficient  dh  of  PZN-4.5%PT 
(111)  cut  vs.  DC  electric  field  at  various  temperatures 
is  plotted  on  Fig.  7.  In  this  case  we  can  observe  sig¬ 
nificantly  different  behavior  for  temperatures  in  vicin¬ 
ity  of  Curie  temperature  Tc  =  423  K.  The  coefficient 
dh  reaches  two  maxima  for  low  bias  fields.  Possibly 
it  could  be  caused  by  presence  of  intermediate  or¬ 
thorhombic  phase.  The  values  of  dh  at  low  tempera¬ 
tures  up  to  353  K  show  slight  change.  Experimental 
value  of  dh  =  68  •  10“12  CN~l  for  T  =  298  K  and 
E  —  0  Vm_i  is  in  good  agreement  with  [1]. 

The  hydrostatic  coefficient  dh  of  P  YN-40%PT  (111) 
cut  vs.  DC  electric  field  at  various  temperatures  is 
plotted  in  Fig.  8.  The  relative  change  of  caused 
by  DC  electric  field  of  1600  kVm"1  was  smaller 
than  6%  at  temperatures  between  303  K  and  393  K. 
The  dh  increases  with  increasing  temperature  up  to 
Tc  =  523  K. 


The  non-linear  effects  of  DC  bias  electric  field  on 
the  hydrostatic  piezoelectric  coefficient  dh  for  variable 
temperature  in  soft  PZT  ceramics,  and  PZN-PT  and 
PYN-PT  relaxor  based  single  crystals  were  studied. 
The  hydrostatic  coefficient  dh  was  determined  by  the 
dynamic  hydrostatic  method.  We  believe  that  the  dis¬ 
crepancy  between  the  values  of  the  d*,  determined  by 
Eq.  (5),  and  those  measured  directly  (using  dynamic 
hydrostatic  method)  is  mainly  caused  by  extrinsic  con¬ 
tributions  of  domain  walls.  This  hypothesis  supports 
the  significant  contribution  of  domain  walls  motion  ob¬ 
served  and  described  in  the  paper  [1]  on  domain  engi¬ 
neered  samples  of  PZN-4.5%PT  single  crystals,  (001) 
cut. 
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For  (001)c  oriented  BaTi03  thin  films,  it  has  been  found  that  epitaxial  constraint  can  result  in  a 
dramatic  increase  in  the  temperature  of  a  tetragonal  (T)  structural  phase  transition.  For  2000-A-thick 
films  grown  directly  on  SrTi03  substrates,  a  7— ►  cubic  (C)  phase  transition  was  found  on  heating  at 
>950  K,  where  the  lattice  constant  changed  smoothly  with  temperature.  It  was  also  found  for  films 
of  the  same  thickness  that  the  7— >  C  transition  is  nearly  restored  to  that  of  bulk  crystals  by  the  use 
of  a  buffer  layer,  which  relaxes  the  epitaxial  constraint.  These  results  provide  evidence  of  an 
epitaxially  induced  high  temperature  structural  phase  transition  in  BaTi03  thin  films,  where  the 
ferroelectric  (internal)  and  structural  (external)  aspects  of  the  phase  transition  are  decoupled.  © 
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The  classic  first  order  ferroelectric  phase  transformation 
is  the  cubic  (C)  to  tetragonal  (T)  one  in  BaTi03  bulk 
crystals,1,2  which  occurs  near  a  Curie  temperature  of  Tc 
=  405  K.  Subsequent  ferroelectric  transitions  occur  to  lower 
symmetry  orthorhombic  and  rhombohedral  phases  with  de¬ 
creasing  temperature.  Thin  films  of  BaTi03  have  also  been 
intensively  investigated  for  applications  in  ferroelectric 
memories,  integrated  capacitors,  and  electro-optic 

^  q  I  I 

modulators.  Both  epitaxial  and  polycrystalline  films 
have  been  studied. 

Investigations  of  epitaxial  BaTi03  films  from  on  SrTi038 
and  MgO9”11  have  reported  the  lack  of  structural  phase  trans¬ 
formations.  This  was  evidenced  by  the  lack  of  anomalies  in 
the  temperature  dependence  of  the  lattice  parameters.  In  par¬ 
ticular,  studies  by  Yoneda  et  al 9  have  found  the  lack  of  a 
7-*C  transition  on  heating  for  films  of  100,  1000,  and 
4000  A  thickness  grown  on  MgO — the  tetragonal  structure 
remained  present  at  all  temperatures  measured  upto  850  K. 
However,  corresponding  electrical  property  measurements 
revealed  a  maximum  in  the  temperature  dependent  dielectric 
constant  ( K )  at  7^410  K,  in  near  agreement  with  that  of 
bulk  BaTi03  single  crystals  and  ceramics.  These  results  have 
recently  been  interpreted  as  the  effect  of  epitaxial  constraint 
that  prevents  external  dimension  changes  of  the  film  due  to 
clamping,  while  permitting  a  ferroelectric  phase  transition 
with  internal  displacive  rearrangements.15 


“^Electronic  mail:  dviehlan@vt.edu 


Recently,  a  Landau-Ginzburg-Devonshire  (LGD)  ther¬ 
modynamic  theory  has  been  developed  by  Pertsev  et  al}6  to 
describe  the  effect  of  epitiaxial  constraint  on  the  structural 
phase  transitions  of  SrTi03  thin  films  grown  with  differing 
misfits  on  various  substrates.  SrTi03  bulk  crystals  undergo  a 
C  — ►  T  transition  on  cooling  at  7  =  105  K  which  is  ferroelastic 
in  nature,  but  not  ferroelectric.  Rather,  it  is  an  incipient  ferro¬ 
electric.  Using  the  LGD  theory,16  the  epitaxial  phase  diagram 
of  (001)c  oriented  SrTi03  as  calculated  for  films  grown  on 
various  cubic  substrates.  It  was  shown  that  decoupled  ferro¬ 
electric  and  ferroelastic  instabilities  could  be  driven  by  epi¬ 
taxial  constraint,  where  various  ferroelectric  phase  transi¬ 
tions  could  occur  as  internal  atomic  rearrangements  within 
ferroelastically  distorted  T  phases. 

The  thermodynamic  investigations  of  constrained  thin 
layers  of  SrTi0316  suggest  that  epitaxy  might  result  in  decou¬ 
pled  ferroelectric  and  ferroelastic  phase  transitions  from  pro- 
totypic  cubic  (m3m)  to  tetragonal  symmetries  in  many  per- 
ovskite  films.  However,  it  has  yet  to  be  shown  whether  there 
is  an  independent  ferroelastic  7— ►  C  structural  phase  transi¬ 
tion  in  epitaxially  constrained  ferroelectric  films  at  higher 
temperatures,  or  not.  In  this  letter,  we  will  report  a  high- 
temperature  T~ >C  phase  transition  (>900  K)  in  epitaxially 
constrained  (001  )c  oriented  BaTi03  films  grown  directly  on 
SrTi03  substrates.  It  will  also  be  shown  for  films  of  the  same 
thickness  that  the  7-+  C  phase  transition  temperature  can  be 
restored  to  that  of  bulk  single  crystals  by  the  use  of  a  buffer 
layer,  which  relaxes  the  epitaxial  constraint. 
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FIG.  1.  Reciprocal  space  scans  for  (001)c  oriented  BaTi03  thin  layers  taken 
along  the  (HOL)  zone  for  (a)  a  2000-A-thick  BaTi03  films  grown  directly 
on  (001)c  SrTi03  substrate;  and  (b)  a  2000-A-thick  BaTi03  films  grown  on 
a  600- A -thick  SrRu03  electrode  (or  buffer  layer),  which  was  grown  on  top 
of  a  (001)c  SrTi03  substrate.  The  values  of  ( HKL )  are  normalized  to  those 
of  BaTi03  single  crystals,  i.e.,  (H.K^L)^ tal  =  (l  ,1,1).  Intensity  is  given  on 
a  log  scale. 

We  have  grown  phase-pure  BaTi03  (BTO)  thin  films  of 
2000  A  thickness  by  pulsed  laser  deposition  onto  (001)c 
single  crystal  substrates  of  SrTi03.  Some  of  the  films  were 
grbwn  directly  on  SrTi03:  deposited  at  670  °C  in  an  oxygen 
ambient  of  20  mTorr  at  a  growth  rate  of  0.7  A/s.  Others 
films  were  grown  on  top  of  a  conducting  perovskite  oxide 
electrode.  SrRu03  (SRO)17  was  chosen  as  the  bottom  elec¬ 
trode  due  to  the  closest  lattice  mismatch  with  the  BTO  struc¬ 
ture:  films  of  SRO  of  500  A  were  deposited  at  600  °C  in  an 
oxygen  ambient  of  100  mTorr,  and  followed  by  the  BTO 
film  deposition  as  described  earlier.  X-ray  diffraction  studies 
were  performed  using  a  Phillips  MPD  system.  Careful  inves¬ 
tigations  were  performed  using  a  high  resolution  x-ray  dif¬ 
fractometer  equipped  with  a  two  bounce  hybrid  monochro¬ 
mator,  an  open  three-circle  Eulerian  cradle,  and  a  doomed 
hot  stage.  A  Ge(220)  cut  crystal  was  used  as  an  analyzer, 
which  had  a  theta-resolution  of  ~0.0068°  (or  0.43  arc  sec). 
The  x-ray  wavelength  was  that  of  Cu^A  =  1 .5406  A)  and 
the  x-ray  generator  was  operated  at  45  kV  and  40  mA.  The 
penetration  depth  in  BaTi03  at  this  x-ray  wavelength  is 
much  greater  than  that  of  the  film  thickness.  The  lattice  con¬ 
stants  of  (001)c  oriented  BaTi03  films  measured  as  a  func¬ 
tion  of  temperature  on  heating  over  the  range  of 
300-1000  K.  The  a -lattice  parameter  of  BaTi03  bulk  single 
crystals  at  7  =  300  K  is  a- 3.99  A,  and  thus  we  use  one  re¬ 
ciprocal  lattice  unit  as  corresponding  to  a*(=b*)=27r/a 
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FIG.  2.  Longitudinal  (200)  profile  for  a  2000-A-thick  BaTi03  films  grown 
directly  on  (001)c  SrTi03  substrate;  and  a  2000-A-thick  BaTi03  films  grown 
on  a  600-A-thick  SrRu03  electrode  (or  buffer  layer),  which  was  grown  on 
top  of  a  (001)c  SrTi03  substrate. 

=  1.547  A-1.  The  mesh  scans  presented  in  this  letter  are  plot¬ 
ted  in  this  reciprocal  unit. 

Figure  1  shows  reciprocal  lattice  scans  taken  along  the 
(HOL)  zone  for  (001)c  BaTi03  films  (a)  grown  directly  on 
SrTi03;  and  (b)  grown  on  a  SrRu03  buffer  layer.  The  mesh 
scans  reveal  two  peaks  for  the  film  grown  directly  on 
SrTi03,  whereas  three  peaks  can  be  seen  for  that  grown  on 
the  SrRu03  buffer  layer.  The  phases  corresponding  to  the 
peaks  are  identified  in  the  figures.  These  results  confirm  that 
the  SrRu03  buffer  layer  has  a  lattice  parameter  close  to  its 
bulk  value  (a=3.98  A),  and  that  the  BaTi03  layer  is  much 
better  lattice  matched  to  it  than  to  the  SrTi03  substrate.  Ac¬ 
cordingly,  for  a  given  thickness,  films  grown  on  a  SrRu03 
buffer  layer  will  be  significantly  more  elastically  relaxed.  It 
is  also  relevant  to  note  that  both  films  had  a  relatively  small 
mosaic  spread  of  —0.2°,  revealing  satisfactory  film  quality. 
Figure  2  shows  longitudinal  (200)  profiles  of  (001)c  BaTi03 
films  grown  directly  on  SrTi03,  and  on  a  SrRu03  buffer 
layer.  The  data  reveal  the  coexistence  of  c  and  a  domains  for 
films  grown  on  the  buffer  layer,  whereas  those  grown  di¬ 
rectly  on  SrTi03  consisted  mostly  of  only  a  domains. 

A  complete  structural  analysis  was  then  done  using  lon¬ 
gitudinal  (200)  and  (202)  profiles  for  both  types  of  films. 
Figure  3  shows  the  lattice  parameters  as  a  function  of  tem¬ 
perature  for  films  grown  (a)  directly  on  the  substrate;  and  (b) 
on  the  buffer  layer.  Comparisons  of  data  in  the  figure  will 
reveal  a  dramatic  difference  in  the  7— >C  phase  transition 
temperature  on  heating  between  these  two  types  of  films.  For 
2000-A-thick  films  grown  on  SrTi03  substrates,  the  T-*C 
phase  transition  occurred  on  heating  near  950  K  (maybe 
slightly  deviated  due  to  fitness  error).  However,  for  an  iden¬ 
tical  2000-A-thick  film  grown  on  SrRu03  buffer  layers,  the 
T— ►  C  phase  transition  occurred  near  470  K — which  is  close 
to  rc=405  K  of  bulk  BaTi03  single  crystals.1,2  These  results 
provide  evidence  of  an  epitaxially  induced  high  temperature 
structural  phase  transition  in  BaTi03  thin  films,  which  is 
relaxed  by  growth  of  an  identical  film  on  a  SrRu03  buffer 
layer. 

It  is  also  important  to  note  several  other  features  in  Fig. 
3.  First,  the  a-lattice  parameter  for  both  films  exhibited  no 
indications  of  a  phase  transition.  Rather,  the  lattice  param¬ 
eters  changed  gradually  with  temperature,  similar  to  a  simple 
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FIG.  3.  Temperature  dependent  lattice  parameters  for  (a)  a  2000-A-thick 
BaTi03  films  grown  directly  on  (001  )c  SrTi03  substrate;  and  (b)  a 
2000-A-thick  BaTi03  films  grown  on  a  600-A-thick  SrRu03  electrode  (or 
buffer  layer),  which  was  grown  on  top  of  a  (001)c  SrTi03  substrate.  Solid 
circles  and  squares  represent  increasing  temperature  and  open  circles  and 
squares  represent  decreasing  temperature. 


thermal  expansion.  This  is  because  any  changes  in  the 
a-lattice  parameter  with  temperature  must  be  invariant  with 
respect  to  corresponding  ones  in  the  substrate  or  buffer  layer. 
Second,  the  cl  a  ratio  is  nearly  equal  for  both  films  at  300  K: 
in  fact,  for  both  types  of  films,  the  relative  changes  in  c!a 
with  temperature  below  470  K  were  gradual  and  of  similar 
magnitude.  These  results  suggest  that  the  ferroelectric  phase 
transformation  temperature  is  nearly  the  same  for  both  films, 
as  Tc  is  generally  excepted  to  be  proportional  to  da}  This  is 
also  supported  by  previous  dielectric  and  structural  studies  of 
(001)c  BaTi03  thin  films  grown  on  Pt/MgO  that  demon¬ 
strated  the  lack  of  a  T— >C  phase  transition,  inspite  of  a 
maximum  in  the  dielectric  constant  at  470  K.9  Our  structural 


data  for  the  films  grown  directly  on  the  SrTi03  substrate  are 
qualitatively  similar — except  we  clearly  show  that  a  7— 
phase  transition  occurs  near  950  K.  In  our  case,  for  films 
grown  directly  on  the  substrate,  the  ferroelectric  (internal) 
and  structural  (external)  aspects  of  the  phase  transition  must 
clearly  be  decoupled. 

In  summary,  2000-A-thick  BaTi03  films  grown  directly 
on  (001)c  SrTi03  substrates  have  a  tetragonal  structure  until 
a  temperature  greater  than  900  K.  However,  an  identical  film 
of  the  same  thickness  grown  on  a  500-A-thick  SrRu03  buffer 
layer  on  top  of  the  substrate  undergoes  a  T-+C  transition  at 
470  K,  in  close  proximity  to  Tc=405  K  of  bulk  BaTi03 
crystals.  These  results  provide  evidence  of  an  epitaxially  in¬ 
duced  high  temperature  structural  phase  transition  in  BaTi03 
thin  films,  where  the  ferroelectric  (internal)  and  structural 
(external)  aspects  of  the  phase  transition  are  decoupled. 
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A  phenomenological  thermodynamic  potential  was  constructed  based  on  the  properties  of  bulk 
BaTi03  single  crystals.  An  eighth-order  polynomial  of  Landau-Devonshire  expansion  was 
employed.  It  reproduces  bulk  properties  including  the  three  possible  ferroelectric  transition 
temperatures  and  their  dependence  on  electric  fields,  as  well  as  the  dielectric  and  piezoelectric 
constants.  Different  from  the  existing  thermodynamic  potential,  it  is  applicable  to  predicting  the 
ferroelectric  phase  transitions  and  properties  of  BaTi03  thin  films  under  large  compressive  biaxial 
strains.  ©  2005  American  Institute  of  Physics.  [DOI:  10.1063/1.2042528] 


I.  INTRODUCTION 

Ferroelectric  BaTi03  ceramics  and  single  crystals  have 
been  extensively  studied  for  many  decades,  and  their  bulk 
ferroelectric  phase  transitions,  crystal  structures,  and  ferro¬ 
electric  properties  are  well  established.  Based  on  the  struc¬ 
tures  and  properties  of  bulk  ceramics  and  single  crystals,  a 
thermodynamic  description  of  the  Landau-Devonshire  theory 
was  developed  by  Bell  and  Cross.1  It  reproduces  quite  well 
most  of  the  single-crystal  structural  and  ferroelectric  proper¬ 
ties  including  the  various  ferroelectric  phase-transition  tem¬ 
peratures.  It  has  been  used  to  predict  the  effect  of  strain  on 
the  phase  transitions  and  ferroelectric  properties  of  BaTi03 
thin  films.2”"4  However,  it  is  only  applicable  to  BaTi03  films 
under  relatively  small  compressive  strains  (^0.4%). 2-4  For 
a  certain  temperature  range,  there  is  no  solution  for  equilib¬ 
rium  polarization  under  larger  compressive  strains,  i.e.,  no 
ferroelectric  phase  transition.  On  the  other  hand,  it  has  been 
shown  that  BaTi03  thin  films  may  be  subjected  to  much 
higher  compressive  strains,  >1%,  due  to  substrate  con¬ 
straints.  For  example,  it  is  recently  discovered  that  BaTi03 
films  can  be  compressively  strained  as  much  as  1.6%,  and 
the  cubic  to  tetragonal  ferroelectric  transition  temperature 
of  BaTi03  thin  films  can  be  increased  to  over  600  °C,5’6  a 
huge  shift  compared  with  bulk  BaTi03  single  crystals 
(-125  °C).7'10 

In  order  to  predict  the  phase  transitions,  domain  struc¬ 
tures,  and  ferroelectric  properties  of  BaTi03  thin  films  under 
such  large  compressive  constraints,  a  thermodynamic  poten¬ 
tial  is  required.  The  main  objective  of  the  paper  is  to  con¬ 
struct  such  a  thermodynamic  potential. 
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II.  PHENOMENOLOGICAL  THERMODYNAMIC 
POTENTIAL 

In  a  phenomenological  description  of  the  ferroelectric 
phase  transitions  in  BaTi03,  the  spontaneous  polarization  P 
=(PUP2,P3)  is  chosen  as  the  order  parameter.  The  Landau- 
Devonshire  potential  is  expanded  as  a  polynomial  of  the  po¬ 
larization  components  Pi(i=  1 ,2, 3). 11,12  In  this  work,  we  em¬ 
ployed  an  eighth-order  polynomial  for  the  Landau- 
Devonshire  potential.  Under  mechanical  stress-free  boundary 
condition,  it  is 

fo  =  «1  {P\  +  P\  +  Pi)  +  an(Pt+P*2+ P$) 

+  an{P\P22  +  P\P\  +  P]P2J  +  «m(P‘  +  P62  +  P\) 

+  *\n\.P]{P\  +  Pi)  +  P\(P\  +  P$)  +  P\iP\  +  ^)1 
+  a\nP\P\P\  +  “1111(^1  +  ^*2  +  ^3) 

+  «n.2[^22 + Pi) + p\(p\ + Pi) + p\tp\ + ^2)] 

+  «1122  (P\P*2  +  P*2P\  +  P\P\) 

+  «1  m{P\P\Pi  +  P\p\p\  +  PiPiPl) .  (1) 

where  all  of  the  coefficients  are  assumed  to  be  temperature 
independent  except  ax.  ax  is  linearly  dependent  on  tempera-, 
ture  and  obeys  the  Curie- Weiss  law.  The  coefficients  ob¬ 
tained  in  this  work  are  listed  in  Table  I.  axx,  aul,  and  ann 
were  fitted  to  the  ferroelectric  transition  temperature  from 
paraelectric  cubic  phase  to  ferroelectric  tetragonal  phase  and 
to  the  spontaneous  polarization  and  dielectric  constant  of  the 
tetragonal  phase.  a12,  au 2,  <*11121  and  axnl  were  fitted  to  the 
properties  of  ferroelectric  orthorhombic  phase  while  the  re¬ 
maining  two  coefficients  a123  and  aU23  were  determined 
from  the  ferroelectric  rhombohedral  phase.  It  should  be 
noted  that  the  existing  thermodynamic  potential1  employed  a 
sixth-order  polynomial.  In  order  to  fit  the  three  ferroelectric 
phase  transitions  as  well  as  the  ferroelectric  properties  as  a 
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TABLE  I.  Coefficients  of  Landau-Devonshire  potential  in  Eq.  (1)  where  T  is  temperature  in  °C. 


Coefficients 

This  work 

Bell  and  Cross8 

Units 

*i 

4.124  X  105(7’-115) 

3.34  XlO^r- 108) 

C-2m2N 

*n 

-2.097  X108 

4.69  X  106(7’- 120) -2.02  X  108 

C"4m6N 

*12 

7.974  X108 

3.23  X  108 

C-4m6N 

*111 

1.294  X109 

-5.52X  107(7"-120)+2.76x  109 

C~6m10N 

*112 

-1.950X  109 

4.47  X109 

C"6  m10  N 

*123 

-2.500  X109 

4.91  X109 

C~6ml0N 

*im 

'  3.863  X  10'° 

0.0 

C’8  ml4N 

*1112 

2.529  X  1010 

0.0 

C'8  m14  N 

*1122 

1.637  X10'° 

0.0 

C“8ml4N 

*1123 

1.367  X  1010 

0.0 

C~8ml4N 

“Reference  1. 


function  of  temperature,  the  coefficients  ah  au,  and  am 
were  assumed  to  be  dependent  on  temperature.  It  is  the  de¬ 
pendence  of  au  and  aul  on  temperature  that  limits  the  ap¬ 
plication  of  the  existing  potential  to  small  compressive 
strains.  For  comparison,  the  coefficients  from  Bell  and 
Cross1  are  listed  in  the  table  as  well. 

With  the  new  coefficients  listed  in  Table  I,  the  potential 
in  Eq.  (1)  yields  the  transition  temperatures: 

^c(Cubic«-+ Tetragonal)  “  1^5  C,  ^c(TetragonalwOrthorhombic)  “  ^ 

<md  = —7 1  C  for  BaTi03  single 

crystals  under  the  stress-free  condition.  The  ferroelectric 
properties  calculated  at  room  temperature  are  compared  in 
Table  II  with  available  experimentally  measured 
values10,13-20  and  those  from  the  coefficients  of  Bell  and 
Cross.1,21  The  polarizations  and  dielectric  constants  as  a 
function  of  temperature  obtained  from  the  potential  are 
shown  in  Figs.  1  and  2,  respectively. 

Since  only  the  second-order  coefficient  is  dependent  on 
temperature  linearly,  the  entropy  change  for  a  phase  transi¬ 
tion  at  the  transition  temperature  Tc  can  be  related  to  polar- 
ization  as  A5=[/^roduct(7’c)-^arent(rc)]/2e0C,11  where  e0  is 
the  permittivity  of  vacuum  and  C  is  the  Curie- Weiss  con¬ 
stant.  The  entropy  changes  A S  at  the  three  transition  points 
calculated  from  the  thermodynamic  potential  are  listed  in 
Table  IEL  The  experimental  data  from  Jona  and  Shirane11  are 
included  in  the  table  for  comparison.  Although  they  were  not 
used  in  the  fitting  of  the  potential,  the  calculated  and  experi¬ 
mentally  measured  values  agree  with  each  other  very  well. 


III.  PHASE  TRANSITIONS  UNDER  ELECTRICAL  FIELD 

Under  the  presence  of  an  external  electric  field  E 
=(ElJE2,E3)i  the  Landau-Devonshire  potential  is  rewritten 
as  follows  in  order  to  take  into  account  its  effect:1,21 

f=fo~ElPl-E2P2-E3P3.  (2) 

Based  on  Eq.  (2),  phase  diagrams  as  a  function  of  tempera¬ 
ture  and  external  applied  electric  field  were  constructed.  Fig¬ 
ures  3(a)-3(c)  show  the  results  for  the  cases  where  electric 
fields  2?=(0,0,£0)>  E=(Eq,0,E0),  and  E=(E0,E0,E0)  are 
applied,  respectively.  The  stable  ferroelectric  phase  corre¬ 
sponds  to  the  minimum  of  /  in  Eq.  (2)  for  a  given  tempera¬ 
ture  and  electric  field.  For  the  applied  electric  field  E 
=(0,0,£0)  our  results  are  in  good  agreement  with  the  experi¬ 
mental  results  from  Fesenko  and  Popov.22 

There  is  a  renewed  interest  in  the  piezoelectric  properties 
of  BaTi03  single  crystals  since  it  is  believed  that  a  signifi¬ 
cant  enhancement  could  be  achieved  through  domain  engi¬ 
neering  by  applied  electric  fields.19,20,23  We  calculated  the 
piezoelectric  coefficient  d3i(/=l,2,3)  through  the  deriva¬ 
tives  of  strains  with  respect  to  the  electric  field  for  a  BaTi03 
single  crystal  polarized  along  the  [001],  [101],  and  [111] 
crystallographic  axes,  respectively.  Under  a  stress-free  con¬ 
dition,  the  strains  in  the  single  crystal  are  the  spontaneous 
strains.  For  example,  the  spontaneous  strain  along  any 
arbitrary  direction  [/,m,/z]  with  respect  to  the  pseudocubic 
cell  is  calculated  through 


TABLE  II.  Properties  of  BaTi03  at  room  temperature  T=25  °C. 


Properties 

This  work 

Others 

Units 

Ps 

0.26 

0.25,b  0.26, c,d  0.265s 

C  m'2 

£33 

188 

168,'  150, b  130±5,f,8'h  193* 

ell 

3600 

2920,'  3600, b  4400±400,f'8'b  3300s 

^33 

86.3 

85.6,'  68.5/  90±5,e  125,'  128, 1 99.8k 

10"12  C  N"1 

^31 

-29.4 

-34.5,'  -33.4,f  -32.5±2,8  -40.8k 

10-|2CN-' 

“Reference  1. 

Reference  17. 

Reference  10. 

hReference  18. 

'Reference  13. 

’Reference  19. 

dReference  14. 

^Reference  20. 

'Reference  15. 
Reference  16. 

Reference  21. 
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Temperature  T  (°C) 

FIG,  I.  (Color  online)  Polarizations  vs  temperature  in  BaTi03  single  crystal 
under  stress-free  condition,  where  £=|£|,  P=(0,0t£,)  in  tetragonal  phase, 
P=(Po,0,Po)  in  orthorhombic  phase,  and  P=(Pr,Pr,Pr)  in  rhombohedral 
phase. 


e[l,m,n]  =  A \l2  +  e22m2  +  e33n 2  +  +  2^/n 

+  2^3  mny  (3) 

where  e?.  is  the  spontaneous  strain  or  transformation  strain 
along  pseudocubic  axes, 

e0n  =  QuP2i  +  QnPl  +  Qx2Pl 
e22  =  Gl2^1  +  Qll^2  +  Q  12^3* 

(4) 

e33  =  Ql2^2  +  Q  12^*2  +  Gll^3 

*12  =  044^ \P 2i  el3  “  Q44P 1 P 3>  ^23  =  044^ 2 P: 3* 

The  electrostrictive  coefficients  2u=0-10>  Gi2=_0*034,  and 
Q44=0’029(C_1  m2)  from  Yamada^4  were  used  for  calculat¬ 
ing  the  piezoelectric  coefficients.  Figures  4-6  show  the  pre¬ 
dicted  piezoelectric  coefficients  d3i(i=l  ,2,3)  at  room  tem¬ 
perature  as  a  function  of  applied  electric  field.  Wada  et  al. 
investigated  the  enhancement  of  piezoelectric  properties  by 
applying  an  electric  field  along  the  [111]  direction  and  found 
that  d33~203  pCN“1  when  £[m]<0.6  MV  nr1,  d33 
~  145  pC  N'1  when  4.0  MV  m_l<£[1U]<5.0  MV  m_1,  arid 
d33  had  its  maximum  of  295  pCN“!  when  1.6  MV  m"1 
<£[111]  <2.6  MV  m-1.19  The  experimental  data  agreed  well 


Temperature  T  (°C) 

FIG.  2.  (Color  online)  Dielectric  constants  in  the  corresponding  ferroelectric 
phases  where  e33  is  along  the  polarization  direction  and  eH  and  £22  arc  along 
the  two  directions  which  are  orthogonal  to  each  other  and  to  the  polariza¬ 
tion.  For  tetragonal  and  rhombohedral  phases,  £n =£22-  For  orthorhombic 
phase  of  P=(P0, <),/>„),  e22  is  along  the  [010]  direction. 


TABLE  m.  Entropy  changes  AS  of  BaTi03  for  the  three  ferroelectric  tran¬ 
sitions. 


This  work 


Others3 

Transition  [J  m-3  K"1]  [cal/(mole  K)]  [cal/(mole  K)] 


Cubic  to  tetragonal 

13535.1 

0.124 

0.12  — 

0.125 

Tetragonal  to  orthorhombic 

9939.8 

0.091 

0.054- 

-0.091 

Orthorhombic  to  rhombohedral 

6970.9 

0.064 

0.04- 

-0.07 

“Reference  11. 


with  our  prediction  (see  Fig.  4).  It  may  be  noted  that  the 
piezoelectric  coefficient  d33  is  the  highest  at  phase-transition 
temperatures.  A  similar  feature  is  found  in  Fig.  5  when  the 
applied  electric  filed  is  along  the  [101]  direction.  There  is  no 
significant  variation  in  the  piezoelectric  coefficients  d33  and 


FIG.  3.  (Color  online)  Phase  diagram  as  a  function  of  temperature  and 
applied  electric  field  (a)  E=(0,0,£0)  so  that  Efooij^o;  (b)  E=(£0,0,£0)  so 
that  \!2Eq]  (c)  E=(£0,£0,£0)  so  that  E^ny=\3EQ. 
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FIG.  4.  (Color  online)  Piezoelectric  coefficient  dv  vs  applied  electric  field 
£[1H]  at  room  temperature. 


FIG.  6.  (Color  online)  Piezoelectric  coefficients  d33  and  rf31  vs  applied  elec¬ 
tric  field  £[00 j}  direction  at  room  temperature. 


d3h  however,  with  applied  electric  field  £*[001]  (Fig-  6).  This 
is  consistent  with  the  fact  that  there  is  no  additional  ferro¬ 
electric  transition  induced  by  the  applied  electric  field  £[ooi]- 


IV.  SUMMARY 

An  eighth-order  polynomial  Landau-Devonshire  poten¬ 
tial  was  proposed  for  bulk  BaTi03  single  crystals.  It  repro¬ 
duces  the  ferroelectric  properties  of  BaTi03  single  crystals  as 
a  function  of  temperature.  The  dependence  of  the  ferroelec¬ 
tric  transition  temperatures  and  piezoelectric  properties  on 
applied  electric  field  calculated  from  the  potential  agrees 
well  with  reported  experimental  data.  It  is  found  that  the 
highest  piezoelectric  coefficient  d33  always  coincides  with  a 
ferroelectric  transition. 


FIG.  5.  (Color  online)  Piezoelectric  coefficient  d2i  vs  applied  electric  field 
Em]  at  room  temperature. 
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Polarization  switching  and  dielectric  temperature  behavior  have  been  investigated  in 
O.SSPKYbi^Nb^^Oj-O^SPbT^  (PYN-PT)  single  crystals.  These  compositions  are  of  interest 
owing  to  their  unusually  high  Curie  temperature  and  high  coercive  fields.  The  dielectric  permittivity 
was  determined  as  function  of  temperature  and  the  polarization  was  measured  as  a  function  of 
applied  field  for  different  drive  frequencies.  Analysis  of  the  experimental  results  with  a  random-field 
model  was  found  to  show  similar  trends  to  recent  reports  for  the  Pb(Mg1/3Nb2/3)03  relaxor 
materials.  The  frequency  dispersion  of  the  maximum  in  dielectric  spectra,  7^,  was  also  analyzed 
with  the  Vogel-Fulcher  relation  yielding  a  high  freezing  temperature  of  577  K.  Polarization 
behavior  was  found  to  follow  a  stretched  exponential  time  dependence  and  a  fractal  dimension  of 
—2.2  was  obtained.  It  is  proposed  that  the  polar  cluster  state  in  PYN-PT  single  ciystals  do  not  freeze 
in  random  orientations  but  rather  arrange  locally  in  preferred  configurations.  The  fraction  of  the 
clusters  decreases  under  increasing  external  field  drive,  resulting  in  the  formation  of  the  percolated 
network  at  higher  fields.  ©  2002  American  Institute  of  Physics,  [DOI:  10.1063/1.1467399] 


I.  INTRODUCTION 

The  lead  based  complex  perovskite  ferroelectric  single 
crystals  are  promising  candidates  for  piezoelectric  transduc¬ 
ers  because  of  their  excellent  dielectric  and  piezoelectric 
properties  close  to  the  moiphotropic  phase  boundary  (MPB) 
but  with  compositions  in  the  rhombohedral  phase  field.1'5 
Longitudinal  electromechanical  coupling  (£33)  values  as 
high  as  —94%  and  piezoelectric  coefficients  ( d 33) 
— 2500  pC/N  have  been  reported  for  (001)  oriented  crystals 
of  Pb(ZnJ/3Nb2/3)03-PbTi03  (PZN-PT)6  and 
PbCMg1/3Nb2/3)03-PbTi03  (PMN-PT),7  making  these  com¬ 
positions  attractive  for  both  actuator  and  ultrasonic  trans¬ 
ducer  applications.  The  magnitudes  of  piezoelectric  and  elec¬ 
tromechanical  properties  are  much  greater  than  those  of  the 
conventional  Pb(Zr,Ti)03  (PZT)  ceramics,  which  are  widely 
used  in  commercial  piezoelectric  devices.  The  implementa¬ 
tion  of  the  PZN-PT  and  PMN-PT  piezocrystals  in  practical 
devices  may  be  limited  by  their  low  Curie  temperature  ( Tc ) 
(150-170  °C)  near  the  MPB  compositions.  PZT  has  a  Curie 
temperature  ranging  from  200°C-350°C  depending  on  its 
doping,  however,  at  present  it  is  extremely  difficult  to  obtain 
PZT  single  crystals  large  enough  to  apply  to  piezoelectric 
devices.8  Naturally,  the  ideal  choice  would  be  a  single  crys¬ 
tal,  which  has  a  high  piezoelectric  and  electromechanical 
coupling  coefficient  like  that  of  PZN-PT  and  PMN-PT  and  a 
high  Tc  similar  to  that  of  PZT.  The  Curie  temperature  Tc  is 
desired  to  obtain  low  temperature  dependent  properties  near 
room  temperature.  Substantial  efforts  have  been  made  to  de¬ 
velop  ferroelectric  single  crystals  with  high  Tc  and  piezo¬ 
electric  properties.  PKIn^Nb^Cyj-PbTiC^  single  crystals 
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were  reported  to  have  a  Tc  of  about  260  °C  with  longitudinal 
electromechanical  coupling  coefficient  in  the  rectangular  bar 
mode  £33  of  about  77%.9,10  Recently,  the  growth  and  charac¬ 
terization  of  a  new  solid  solution  system  based  on 
Pb(Yb1/2Nb1/2)03-PbTi03  (PYN-PT)  has  been  reported.11’12 
The  Tc  near  MPB  (—0.46  PT)  in  this  system  is  around 
337  °C  and  piezoelectric  coefficients  as  high  as  2500  pC/N 
can  be  obtained. 

In  this  article,  we  extend  our  studies  to  explain  the  un¬ 
derlying  mechanisms  responsible  for  the  observed  properties. 
PYN-PT  crystals  close  to  MPB  with  45%  PT  were  chosen  in 
this  study.  Dielectric  measurements  as  a  function  of  tempera¬ 
ture  for  various  frequencies  shows  dispersion  around  Tc, 
typically  associated  with  ferroelectric  relaxorlike  behavior. 
Consequently,  an  in-depth  study  was  undertaken  to  charac¬ 
terize  this  relaxation  behavior.  Relaxors  are  believed  to  ex¬ 
hibit  highly  nonlinear  properties  due  to  a  dynamical  micro¬ 
scopic  polarization  existing  above  a  freezing  temperature. 
However,  the  polarization  can  be  put  in  a  metastable  state 
through  the  establishment  of  a  longer  range  cooperation  un¬ 
der  the  application  of  an  applied  electric  field.13"15  Recently, 
relaxor  behavior  has  been  suggested  to  be  a  consequence  of 
quenched  random  fields  due  to  the  defects,  chemical  inho¬ 
mogeneous,  and  ordering,  which  is  frozen-in  from  the  high 
temperature  processing.16,17  The  concentration  of  these  de¬ 
fects  is  sufficient  to  cause  breaks  in  the  translational  symme¬ 
try  of  the  lattice,  and  thereby  limit  the  development  of  long- 
range  polar  order  without  the  application  of  an  external  field 
to  override  those  internal  random  fields.18  The  random  inter¬ 
action  between  the  clusters  results  in  polarization  freezing 
and  glasslike  characteristics  on  reduction  of  the  thermal 
energy.16,17  Analysis  of  the  excess  polarization  obtained  from 
the  polarization-electric  field  (P-E)  hysteresis  loops  can  be 
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(») 


(b) 


FIG.  1.  Determination  of  the  freezing  temperature  in  the  system  PYN-PT 
for  x«0.45.  (a)  Dielectric  constant  and  dielectric  loss  as  a  function  of  tem¬ 
perature  at  2,  6,  10,  20,  40,  60,  and  100  kHz  with  heating  (solid  line)  and 
cooling  (dot  line)  processes,  (b)  as  a  function  of  frequency.  The  line  is 
analysis  of  this  data  with  Eq.  (1)  using  the  Levenberg-Marquadt  nonlinear 
analysis  program. 


used  to  establish  the  relaxation  rates  and  fractal  dimension¬ 
ality  of  the  polar  clusters.  It  is  noteworthy  that  in  an  inves¬ 
tigation  in  PZN-PT  crystals  by  Yu  and  Randall  direct  evi¬ 
dence  of  fractal  domain  behavior  was  observed  under  various 
electrical  field  applications.19 

The  objective  of  this  article  is  to  use  techniques  that 
have  been  applied  to  PZN-PT  and  PMN-PT  materials  to  con¬ 
trast  the  PYN-PT  relaxor  under  the  analysis  of  random-field 
models. 

II.  EXPERIMENT 

High  quality  rhombohedral  single  crystals  of  PYN-xPT 
(x=0.45)  were  grown  using  the  conventional  high  tempera¬ 
ture  solution  method.11  The  samples  used  in  this  study  were 
(001)  oriented  and  had  a  dimension  of  2.5 X  2.5 X  0.4  mm3. 
High  temperature  dielectric  behavior  was  determined,  from 
room  temperature  to  500  °C,  using  a  multifrequency  LCR 
meter  (HP  42 84 A).  The  data  was  recorded  using  a  General 
Purpose  Interface  Bus  (GP-IB)  interface.  Polarization- 
electric  field  hysteresis  loops  at  high  fields  were  measured  by 


FIG.  2.  Polarization-electric  field  hysteresis  loops  measured  at  frequencies 
0.1,  1,  10,  and  100  Hz,  at  the  maximum  applied  electric  field  level  of  50 
kV/cm. 


using  a  modified  Sawyer-Tower  circuit  driven  by  a  lock-in 
amplifier  (Stanford  Research  Systems,  Model  SR830).  High 
electric  fields  were  applied  via  a  Trek  609C-6  high  voltage 
dc  amplifier. 

III.  RESULTS  AND  DISCUSSION 

A.  Temperature  dependence  of  dielectric  constant 

Figure  1(a)  shows  the  dielectric  constant  and  dielectric 
loss  of  a  poled  (001)  oriented  PYN-0.45PT  single  crystal  as 
a  function  of  the  temperature.  The  data  in  the  figure  is  shown 
for  frequencies  2,  6,  10,  20,  40,  60,  and  100  kHz  with  in¬ 
creasing  and  subsequently  decreasing  temperature  (the  dot¬ 
ted  line  data  represents  cooling  and  the  solid  line  data  corre¬ 
sponds  to  the  heating).  Two  peaks  are  clearly  shown  in  the 
heating  process  around  —160  ( Tr )  and  — 325  °C  ( Tc ).  The 
first  dielectric  anomaly  is  believed  to  result  from  the  break¬ 
ing  of  the  metastable  macrodomain  structure  and  formation 
of  the  microdomains  on  the  scale  of  the  local  defect 
structure.20  This  transition  is  irreversible  in  the  absence  of  an 
external  electric  field.  This  behavior  indicates  that  in  the  ab¬ 
sence  of  an  external  electric  field,  the  polar  regions  are  dy¬ 
namic  and  contribute  to  the  dielectric  relaxation.  Corre¬ 
spondingly,  the  dielectric  loss  has  two  peaks  with  increasing 
temperature.  The  dielectric  constant  in  the  low  temperature 
region  of  25  0 C<  T<  100  °C  was  found  to  be  independent  of 
the  frequency,  and  the  magnitude  rises  from  1700  to  2400. 
[This  part  is  enlarged  and  inserted  in  Fig.  1(a)].  This  shows 
that  temperature  dependence  of  dielectric  behavior  is  rela¬ 
tively  low  for  ferroelectric  piezoelectrics,  fie/  fiT-9/  °C,  it  is 
significantly  low  compared  to  the  MPB  composition  of 
0.7PMN-0.3PT  (de/dT-445/°C  in  the  range  of  [25,rr]).21 
Thus,  in  PYN-0.45PT,  we  have  nearly  temperature  indepen¬ 
dent  properties  in  the  vicinity  of  the  room  temperature.  In  the 
temperature  range  of  Tr<T<Tc ,  frequency  dependent  per¬ 
mittivity  becomes  evident.  On  a  subsequent  cooling  process. 
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FIG.  3.  Plot  of  remanent  polarization  determined  from  the  hysteresis  loops  at  various  maximum  applied  electric  as  a  function  of  the  frequency  (a)  10,  (b)  20, 
(c)  30,  (d)  40,  and  (e)  50  kV/cm.  The  data  in  these  figures  was  analyzed  using  Eq.  (2)  and  the  fit  parameters  are  shown  in  the  plot. 


strong  ferroelectric  relaxor  characteristics  were  observed,  that  found  in  perovskite  based  relaxor  solid  solutions.  The 
with  no  secondary  transformation  at  lower  temperature  Tr .  A  frequency  dependence  of  the  dielectric  constant  is  weak  and 

notable  feature  in  Fig.  1(a)  is  the  dispersion  in  the  dielectric  the  shift  in  the  Tc  is  —4.5  °C  for  the  range  of  frequencies 
constant  around  Tc  and  a  diffused  phase  transition  similar  to  used  in  the  measurement 
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field. 


Previously,  it  has  been  established  that  a  relaxor  ferro¬ 
electric  is  a  polar-glassy  system  which  has  thermally  acti¬ 
vated  polarization  fluctuations  above  a  static  freezing 
temperature.16  The  Vogel-Fulcher  relationship  can  be  used  to 
determine  the  freezing  temperature  in  ferroelectric  relaxor 
materials.22  Figure  1(b)  shows  rmax  as  a  function  of  the  fre¬ 
quency.  The  data  in  this  figure  was  analyzed  using  the 
Vogel-Fulcher  relation 

f=fo  exp[  -Ea  IkiT^-Tf)^  (1) 

where  Tf  is  the  static  freezing  temperature,  Ea  is  the  activa¬ 
tion  eneigy,  and  /0  is  the  preexponential  onset  frequency 
(typically  — 1012  Hz),  yielding  a  freezing  temperature  (Tf)  of 
—577  K  and  the  activation  energy  Ea  of  0.065  eV.  The  freez¬ 
ing  temperature  of  PYN-0.45PT  is  much  higher  than  that  of 
the  PMN-0.3PT  single  crystals,  typically  found  to  be  —370 
K.21  The  close  fitting  of  the  data  in  Fig.  1(b)  suggests  that  the 
relaxor  behavior  in  the  PYN-PT  system  is  analogous  to  a 
dipolar  glass.  Note  the  ratio  Tf/T ^  can  be  related  to  extent 
of  interactions  between  the  clusters.22  Further  insight  into  the 
existence  of  the  glassy  state  can  be  obtained  from  the  time 
dependence  of  the  polarization,  and  is  discussed  below. 

B.  Dynamics  of  polarization  relaxation 

Polarization  switching  occurs  through  stabilization  and 
growth  of  micropolar  clusters  in  the  relaxorlike  materials.23 
The  process  follows  irreversible  thermodynamic  models, 
which  takes  into  account  the  thermostatic  and  thermoelastic 
contributions  of  the  first  and  second  order  of  lattice  defects. 
In  the  presence  of  quenched  disorder,  the  polarization 
switching  process  can  deviate  significantly  from  equilibrium 
conditions.  The  fundamentals  of  these  concepts  have  been 
discussed  by  Imry  and  Ma.24  Under  electrical  drive  condi¬ 
tions  greater  than  that  of  the  coercive  field,  a  breakdown  of 
the  long-range  ferroelectric  order  occurs,  which  results  in  the 
formation  of  the  polar  clusters.  The  functional  form  which 
has  been  used  to  analyze  the  relaxation  in  random-field  mag¬ 
netic  and  glassy  systems  is  a  stretched  exponential,  so  in  the 
ferroelectric  case,  we  expect  a  similar  form,  namely:17,24,25 

Pr(W)  =  P0exp[-(l/a>fc),~'’],  (2) 


FIG.  5.  Plot  of  the  characteristic  time  constant  as  a  function  of  the  applied 
electric  field. 


where  Pr  is  remanent  polarization,  to  is  a  drive  angular  fre¬ 
quency,  tc  is  a  characteristic  relaxation  time,  and  n  is  an 
exponent.  Figure  2  shows  typical  polarization  hysteresis 
loops  of  PYN-0.45PT  in  the  (001)  direction,  determined  at 
0.1,  1,  10,  and  100  Hz  at  the  maximum  electric  field  level  of 
50  kV/cm.  Similar  data  was  taken  for  the  maximum  applied 
electric  field  levels  of  10,  20,  30,  and  40  kV/cm.  The  rema¬ 
nent  polarization  determined  from  the  hysteresis  loops  at 
various  maximum  applied  electric  fields  was  plotted  as  func¬ 
tion  of  the  frequency  and  is  shown  in  Figs.  3(a)-3(e).  The 
data  was  analyzed  using  Eq.  (2)  and  the  fitted  parameters  are 
shown  in  the  respective  plots.  Good  correlations  between  the 
experimental  data  and  the  model  can  be  seen  in  these  figures. 


C.  Fractal  dimensionality 

The  exponent  n  in  Eq.  (2)  can  be  correlated  to  the  di¬ 
mensionality  of  the  fractal  cluster  order  using  the  relation¬ 
ship  proposed  by  Huse  et  al.26  as 

n  =  (d- 2)/(d-l).  (3) 

The  value  of  d  calculated  using  Eq.  (3)  from  the  data  in  Figs. 
3(a)-3(e)  is  —2  and  is  plotted  in  Fig.  4  as  function  of  the 
applied  E  field.  The  magnitude  of  d  shows  a  systematic  de¬ 
crease  with  the  E  field.  The  data  in  Fig.  4  was  fitted  using  an 
empirical  relation  of  the  form 

d=d0  +  dx  expf-(£-£jtf)/rf2],  (4) 

where  d0  is  the  initial  value  of  d  under  low  level  fields,  E ^ 
is  the  random  field,  dx  is  the  exponential  amplitude  prefactor, 
and  d2  is  the  exponential  relaxation  constant.  The  analysis 
yielded  values  of  dQ—2. 09,  10  kV/cm,  ^  =  0.15,  and 

^2  =  8.89  kV/cm.  The  empirical  fitting  to  an  exponential  de¬ 
caying  function  can  be  interpreted  as  that  fraction  of  the 
polar  clusters  and  can  be  reversed  in  various  external  drives. 
The  fractal  dimensionality  is  associated  with  the  growth  of 
the  clusters  and  a  network  formation  of  these  clusters  that 
can  undeigo  percolation  at  high  fields. 
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TABLE  L  Comparison  of  the  important  random-field  parameters  in  PYN-0.45PT  to  PMN-0.3PT  ceramic. 


Crystal 

n 

d 

W 

Eo 

(kV/cm) 

£rp 

(kV/cm) 

Tf 

do 

d> 

d2 

PYN-0.45PT 

(crystal) 

1.09  (40  kV/cm) 

2.2 

2.2  X  10"’ 

257 

10 

577  K 

2.09 

0.15 

8.89 

PMN-0.3PT 
(ceramic)  (Ref.  27) 

0.1554  (1  kV/cm) 

2 

2.5X10" 12 

138 

1.61 

366.1  K  (Ref.  21) 
(crystal) 

2.18 

0.0265 

0.305 

D.  Relaxation  time  constant 

A  characteristic  time  constant,  determined  from  the 
analysis  of  experimental  data  shown  in  Figs.  3(a)-3(e),  by 
using  Eq.  (2),  is  shown  in  Fig.  5.  The  data  was  analyzed  by 
using  the  modified  Arrhenius  equation  given  as 

tc=t0exp(-E0/E),  (5) 

where  tQ  is  the  time  required  for  cluster  nucleation  and  E0  is 
an  activation  field.  The  analysis  yielded  values  of  f0~2.2 
X  10~9  s  and  an  activation  field  of  257  kV/cm.  The  value  of 
E0  is  slightly  higher  than  that  reported  previously  by  analysis 
of  the  frequency  dependence  of  the  coercive  field.  In  general, 
a  value  for  E0  of  ~  100  kV/cm  is  obtained,  which  is  approxi¬ 
mately  half  of  the  magnitude  obtained  by  the  analysis  of  the 
polarization  dependence  for  E>EC  in  the  PYN-PT  system. 

A  comparison  of  the  random-field  parameters  in  a  PYN- 
0.45PT  single  crystal  to  PMN-0.3PT  ceramic27  is  given  in 
Table  I. 

IV.  CONCLUSIONS 

This  article  reported  the  presence  of  the  polar  glassy 
state  in  high  Tc  PYN-PT  single  crystals.  Polarization  reversal 
and  frequency  dependent  obeys  the  trends  of  the 
random-field  model.  A  correlation  between  the  experimental 
data  and  models  indicates  the  presence  of  a  high  concentra¬ 
tion  of  quenched  disorder  in  the  PYN-0.45PT  system.  High 
temperature  growth  of  the  crystals  followed  by  quenching 
results  in  freezing  of  the  high  temperature  thermodynamic 
equilibrium  defects.  Randomness  and  disorder  are  created  by 
mixed  interactions  and  competition  leading  to  cooperative 
freezing  process.  At  temperatures  around  Tf,  many  of  these 
random  polar  defects  arrange  themselves  in  locally  corre¬ 
lated  clusters  and  domains,  which  participate  in  the  switch¬ 
ing  process  as  a  whole.  The  results  of  this  article  demonstrate 
the  usefulness  of  the  random-field  approach  in  deterministi¬ 
cally  understanding  the  relaxation  behavior  of  the  high  Curie 
temperature  piezocrystals,  and  indicate  the  importance  to 
build  up  a  large  database  across  the  perovskite  relaxors  to  aid 
the  design  of  new  piezoelectric  crystals. 
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